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Abstract 

A  number  of  interstellar  precursor  missions  have  been  designed,  but  they  have  all 
made  use  of  unproven  future  technologies.  This  research  focuses  on  designing  an 
interstellar  precursor  mission  to  the  gravitational  lens  that  can  be  manufactured  solely 
with  today’s  technology.  The  design  includes  an  examination  of  trip  time  and  trajectory,  a 
system  trade  between  subsystems,  a  launch  vehicle  analysis  and  an  overall  spacecraft 
performance  analysis. The  trip  time  was  found  to  be  approximately  108  years  using  the 
NEXT  thruster  and  GPHS  power  source.  The  GPHS  provided  approximately  4  kW  of 
power  at  BOL.  The  communications  subsystem  relied  on  a  12  meter  high  gain  antenna, 
which  also  serves  as  the  primary  payload,  and  a  2  meter  medium  gain  antenna.  Finally,  the 
estimated  mission  cost  is  3-5  billion  USD 
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GRAVITATIONAL  LENS:  DEEP  SPACE  PROBE  DESIGN 


1  Introduction 


1.1  The  Focal  Space  Mission 

The  Focal  Space  Mission  (FSM)  revolves  around  the  concept  of  gravitational  lensing 
proposed  by  both  Albert  Einstein  and  Orest  Khvolson  [16]  .  Gravitational  lensing  occurs 
when  electromagnetic  waves  travel  in  the  vicinity  of  massive  objects.  The  gravitational 
field  of  the  object  bends,  or  lenses,  the  waves  and  causes  them  to  converge  to  a  focal  point 
where  they  can  be  resolved  [16].  Any  significantly  sized  astronomical  object  exhibits  this 
phenomenon,  which  was  confirmed  in  1979  by  Walsh,  Carswell,  and  Weymann  [17]. 
Furthermore,  this  lensing  effect  greatly  magnifies  objects  from  a  background  source. 
Utilizing  the  sun  as  a  gravitational  lens  provides  an  unmatched  opportunity  to  observe 
both  neighboring  and  distant  astronomical  phenomena  in  exceptional  detail. 

The  primary  mission  of  the  FSM  will  be  to  use  radio  telescopy  in  concert  with  the 
sun’s  gravitational  lens  to  image  nearby  solar  systems  and/or  other  pertinent  astronomical 
objects,  e.g.  the  Cosmic  Microwave  Background  (CMB)  radiation. 

1.2  Thesis  Objective 

In  Deep  Space  Flight  and  Communications  Claudio  Maccone  uses  the  mass,  radius, 
mean  density,  and  Schwarzschild  radius  of  the  sun  to  find  its  minimum  focal  distance.  His 
result  was  548  Astronomical  Units  (AU)  or  13.86  times  the  Sun-to-Pluto  distance  [18]. 
With  that.  Dr.  Maccone  gives  an  in-depth  review  of  the  physics  necessary  to  complete  this 
mission;  however,  there  has  not  been  a  fully  developed  system  design  of  a  spacecraft  to 
move  this  idea  forward.  Therefore,  the  two  aims  of  this  thesis  are  to  provide  a  FSM 
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mission  analysis  and  FSM  spacecraft  design.  However,  unlike  other  notable  innovative 
interstellar  precursor  designs,  e.g.  Project  Prometheus,  the  FSM  design  will  consist  solely 
of  commercially  available  technology  (COTS).  This  will  allow  current  and  future 
investigators  to  understand  the  limits  of  current  deep  space  capabilities  so  that 
recommendations  can  be  made  for  going  forward.  Finally,  it  is  important  to  note  the 
design  done  in  this  study  is  limited  in  scope  owing  to  the  fact  that  it  is  a  masters  thesis  and 
does  not  have  the  requisite  resources  for  a  full-up  spacecraft  design.  Consequently,  there 
is  no  fault-tree  analysis  and  certain  subsystems  were  explored  in  more  depth  than  others. 

The  mission  analysis  will  explore  the  feasibility  of  traveling  to  the  focal  point  of  the 
gravitational  lens  and  provide  estimates  for  what  is  required  to  perform  the  mission.  This 
includes  examining  the  AC  required,  the  trip  time,  and  the  trajectory  of  the  spacecraft.  On 
the  other  hand,  the  spacecraft  design  will  include  a  system  trade  between  subsystems,  a 
spacecraft  model,  and  performance  measures.  The  spacecraft  design  will  include  a  system 
trade  between  subsystems,  a  spacecraft  model,  and  performance  measures.  Therefore  the 
five  deliverables  are  a  trajectory  analysis,  subsystem  trade,  spacecraft  model,  spacecraft 
performance  analysis  and  launch  vehicle  analysis. 

More  specifically,  the  trajectory  analysis  will  look  into  the  AV  required,  trip  time, 
and  trajectory.  The  subsystem  trade  will  look  at  the  attitude  control  system,  propulsion 
system,  electrical  power  system,  environmental  control  system,  and  communications 
system,  and  of  these  subsystems  the  power,  propulsion  and  communications  subsystems 
will  receive  the  most  scrutiny  because  they  are  the  greatest  challenges  to  mission  success. 
The  two  remaining  systems  will  loosely  adopt  solutions  used  on  other  comparable 
missions.  Next,  the  spacecraft  model  will  display  both  the  internal  spacecraft  layout  and 
the  external  spacecraft  layout.  The  external  spacecraft  layout  will  include  models  of  the 
spacecraft  in  the  stored  and  deployed  configuration.  Following  that,  the  spacecraft 
performance  analysis  will  address  whether  the  mission  can  be  completed  with  the 
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assembled  spacecraft.  The  launch  vehicle  analysis  will  explore  things  from  mass  and 
volume  constraints  to  the  AV  it  provides.  The  final  and  unifying  goal  of  the  study  is  to 
understand  whether  there  is  sufficient  technological  maturity  to  complete  an  interstellar 
precursor  mission  with  COTS.  Concluding  this  section,  Table  1.1  outlines  the  results  and 
where  each  can  be  found  in  the  document  while  Table  4.1  outlines  the  mission  payloads. 


Table  1.1:  Thesis  Objectives  Summary 


Thesis  Objectives 

Objective  Breakout 

Result 

Section  Number 

AV 

51.79  km/sec 

4.5. 1.4 

Trajectory  Analysis 

Trip  Time 

108.7years 

4.5. 1.4 

Trajectory 

Spiral 

4.5.1 

Star  Tracker 

Attitude  Control  System 

IMU 

Att.  Thruster 

4.3.5 

Propulsion  System 

NEXT 

4.3.4 

Subsystem  Trade 

Power  System 

GPHS 

4.3.3 

Environmental  Control 
System 

MU 

RHU 

Radiators /w  Louvers 

4.3.6 

Communication  System 

12-m  HGA 

2-m  MGA 

4.3.2 

Model 

External 

4.4 

Internal 

4.4 

Data  Rate 

Variable 

4.5.4 

Dry  Mass 

2345.93  kg 

4.5.2 

Performance  Analysis 

Power  Consumption 

188.75  W 

4.5.3 

Cost 

$3-5  Billion  (FY11  US$) 

4.5.5 
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2  Background  &  Theory 


2.1  Overview 

The  first  step  in  a  novel  scientific  study  is  a  comprehensive  literature  review.  Thus, 
this  chapter  is  present  to  positively  assert  the  unique  nature  of  the  FSM  mission  by 
analyzing  NASA’s  previous  deep  space  missions,  where  deep  space  is  defined  to  be 
greater  than  5  AU,  and  to  present  a  basis  from  which  to  design  the  FSM.  Accompanying 
the  analysis  of  past  space  missions  is  a  rigorous  treatment  of  the  science  and  history  of 
gravitational  lensing,  since  it  is  of  critical  importance  to  the  subsequent  work  presented, 
along  with  a  discussion  of  secondary  missions  and  payloads. 

2.2  Gravitational  Lens 

It  is  only  recently  that  the  phenomenon  of  gravitational  lensing  has  come  to  be 
understood.  As  stated  earlier,  it  was  posited  in  the  early  20,/'  century  and  not  verified  until 
1979.  Einstein  introduced  the  world  to  the  physics  of  the  gravitational  lens  in  his  seminal 
Science  journal  article  titled  Lens-like  Action  of  a  Star  by  the  Deviation  of  Light  in  the 
Gravitational  Field  [16].  The  geometry  of  the  sun’s  gravitational  lens  is  displayed  in 
Figure  2.1,  which  is  not  to  scale.  It  must  be  noted  that  A  does  not  need  to  be  a  star  and 
indeed  it  can  be  any  observable  phenomenon.  In  the  figure,  ao  corresponds  to  the 
deflection  angle,  x  to  the  distance  from  the  line  that  intersects  stars  A  and  B,  R0  to  the 
radius  of  the  sun,  and  D  corresponds  to  the  radial  distance  from  the  sun  to  its  focus. 

In  his  article,  Einstein  put  forth  that  at  the  focal  point  D  ,  of  star  B,  the  star  A  does  not 
appear  to  be  a  star  at  all  but  a  “luminous  circle”  with  angular  radius  (L  where 


P  = 


[16i 


(2.1) 
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Figure  2.1:  Geometry  of  the  Sun’s  Gravitational  Lens.  Note:  figure  is  not  drawn  to  scale. 


and 


«o 


4  GMS, 
c2r 


[18] 


(2.2) 


With  the  introduction  of  the  first  equations,  it  is  important  to  note  that  all  equations 
listed  in  this  thesis  use  base  International  System  units  unless  otherwise  noted. 

In  the  Equation  2.2 ,  G  represents  the  universal  gravitational  constant,  Msun 
represents  the  mass  of  the  sun,  c  represents  the  speed  of  light,  and  r  represents  the  radial 
distance  at  which  the  light  passes  over  the  sun.  Furthermore,  anything  in  this  luminous 
circle  will  have  an  increased  apparent  brightness  due  to  the  gravitational  field  of  B. 
However,  this  is  limited  by  x,  the  distance  from  the  line  that  intersects  the  center  of  star  A 
and  B,  because  as  the  observer  moves  farther  away  from  the  line  of  intersection  the 
magnification  decreases  [16].  This  ratio  of  magnification,  q,  is  dictated  by  the  expression 


q  = 


l  1  +  2P 


[16] 


(2.3) 


x2 

V  +  4P 


where 


/  =  yja^DRo  [16] 


(2.4) 


Looking  at  Equation  2.3  the  reader  will  immediately  notice  that  to  have  a 
considerable  level  of  magnification,  q  »  1,  x  must  be  very  small  compared  to  /.  Thus,  the 
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value  j2  can  be  excluded  from  the  calculation  for  the  cases  germane  to  this  thesis  yielding 
the  equation  [16] 

q=-  [16]  (2.5) 

x 

To  maximize  the  magnification  x  must  be  minimized  and  l  must  be  maximized.  In  his 
article,  Einstein  noted  the  interesting  quality  that  as  the  displacement  from  the  line  of 
intersection  approaches  zero  the  magnification  goes  to  infinity  [16].  The  first  quantity  x 
can  only  be  affected  by  the  attitude  and  trajectory  of  the  spacecraft.  The  second  quantity  l 
is  a  function  of  two  variables,  one  being  the  distance  from  the  focusing  star  D  and  the 
second  being  the  deflection  angle  a().  As  the  distance  along  the  line  of  intersection 
increases  the  magnification  increases  proportionally  to  Vd.  This  means  that  the 
gravitational  focus  is  a  line  extending  outward  to  infinity  from  the  minimum  focal  point 
[18].  So  maximizing  this  value  is  also  dependent  upon  the  trajectory  of  the  spacecraft.  On 
the  other  hand,  the  deflection  angle  ao  has  a  definite  maximum.  It  occurs  at  the  minimum 
possible  radius.  In  the  case  of  the  Sun  this  minimal  radius  is  just  the  radius  of  the  sun  R0. 
Therefore,  electromagnetic  waves  deflected  at  Rq  will  maximize  ao,  /,  and  as  a  result  q. 
However,  at  this  point  and  with  the  information  provided  thus  far,  utilizing  the 
trigonometry  of  Figure  2.1  the  angle  a  can  be  related  to  the  distance  from  the  focusing 
star,  D  and  the  radius  of  the  sun  Rq 


tan  a  -  —  [18] 

D 

which  through  the  small  angle  approximation  and  substitution  becomes 

D= ^  [18] 


(2.6) 


(2.7) 


where  the  Schwarzschild  Radius  r„  is 


2  GMsun 

rs  = - o —  [!8] 


(2.8) 


6 


This  yields  a  minimum  focal  distance  of  542  AU  [18].  Yet,  caution  should  be  exercised 
with  this  result  because  due  to  the  sun’s  chaotic  corona  the  focal  distance  for  specific 
frequencies  can  be  pushed  out  hundreds  of  astronomical  units  or  even  completely 
eliminated. 

The  sun’s  corona,  according  to  the  Baumbach- Allen  Model,  is  composed  primarily  of 
three  sources:  K-corona,  F-corona,  and  E-corona  [18].  Utilizing  this,  the  most  accurate 
model  of  the  corona,  Anderson  and  Turyshev  produced  Equation  2.9  for  the  deflection 
angle  due  to  the  corona. 


ar 


6.36  MHz) 


)2H£f 


+  228 | ^ 
b 


+  u[£ 

b 


[19] 


(2.9) 


where  b  is  the  impact  parameter,  or  radius  from  the  center  of  the  sun  to  the  point  of 
deflection,  and  v  is  the  frequency  of  the  wave  to  be  observed.  So  the  total  deflection 
becomes 

&  total  ®  ® corona  M  8 1  (2.10) 


This  implies  that  the  total  deflection  angle  is  a  function  of  the  impact  parameter  as  well  as 
frequency.  This  brings  about  a  very  important  result-the  critical  frequency.  If  the 
frequency  of  the  wave  is  below  this  critical  frequency  for  a  given  impact  parameter  then 
the  focus  disappears  [18].  The  critical  frequency  is  given  by  the  equation 
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Plotting  this  equation  (Eq.  2.11),  by  varying  the  impact  parameter  b  from  1  to  15, 
produces  Figure  2.2.  Also,  with  the  equation  for  the  critical  frequency  a  new  equation  for 


the  location  of  the  gravitational  focus  emerges.  Equation  2.12  clearly  demonstrates  that  if 
the  frequency  is  less  than  the  critical  frequency  of  its  impact  parameter  then  the  focus 
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Figure  2.2:  Critical  Frequency  as  a  Function  of  Impact  Parameter 

disappears  because  its  point  of  intersection  with  the  focal  axis  is  negative. 

542  ( — )2 

Dtotai  » - [18]  (2.12) 

v2cntuJb) 

1 

Noting  that  the  focus  existence  is  dependent  upon  the  impact  parameter  and  frequency  Dr. 
Maccone  developed  focusing  conditions  based  on  Baumbach- Allen  model.  To  do  this  Dr. 
Maccone  split  the  Baumbach- Allen  model  into  its  three  constituent  pieces  and  found 
focusing  conditions  for  each.  The  E-Corona  corresponds  to  impact  parameters  in  the 
range  of  1  to  1.3  solar  radii  [18].  The  K-Corona  corresponds  to  impact  parameters  in  the 
range  of  1.3  to  3  solar  radii  [18].  Finally,  the  F-Corona  corresponds  to  impact  parameters 
from  three  solar  radii  to  infinity  [18].  Looking  at  the  equations  for  coronal  deflection 
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Figure  2.3:  Focusing  Conditions 


(Eq.  2.9)  the  E-Corona,  K-Corona,  and  F-Corona  are  represented  by  the  three  terms  in 
brackets  and  in  that  order.  After  this  separation  Dr.  Maccone  noted  that  the  observing 
frequency,  which  is  a  function  of  the  impact  parameter,  for  each  piece  of  the  corona  must 
be  greater  than  its  respective  critical  frequency  [18].  By  plotting  the  equations  for  the 
observing  frequency  of  each  of  the  coronal  constituents  versus  the  critical  frequency  the 
focusing  conditions  can  be  found.  This  is  shown  in  Figure  2.3. 

Each  focusing  condition  is  located  by  finding  the  intersection  between  the  critical 
frequency  and  coronal  constituent.  For  the  E-Corona  the  intersection  point  occurs  at 
(1.588,1 1.07).  This  implies  that  focusing  occurs  for  impact  parameters  less  than  1.588 


9 


solar  radii  and  frequencies  above  11.07  GHz  in  this  zone.  The  K-Corona’s  intersection 
point  is  located  at  (10,  .75)  .  Finally,  the  F-Corona  has  focusing  for  all  observing 
frequencies  since  there  is  no  intersection. 

Aggregating  all  of  the  information  presented  thus  far  yields  the  most  important 
equation  of  the  section-Equation  2.13.  From  Einstein  it  was  apparent  that  the 
gravitational  focus  is  not  a  point  but  a  line  extending  to  infinity  from  the  minimum  focal 
distance  [16].  Referring  back  to  Equation  2.7  the  sun  sans  corona,  or  naked  sun,  focal 
distance  minimum  occurs  when  the  waves  graze  the  naked  sun  at  R0.  The  introduction  of 
the  Baumbach- Allen  corona  model  resulted  in  Equation  2.12.  This  new  equation 
necessitated  the  exploration  of  focusing  conditions.  To  explore  these  conditions  the 
corona  model  was  broken  down  into  its  three  constituent  pieces.  From  this  analysis  it  was 
clear  that  even  with  the  introduction  of  the  corona  the  logic  gleaned  from  Equation  2.7 
still  holds  meaning  the  minimum  focal  distance  will  occur  for  waves  travelling  closest  to 
the  naked  sun  in  the  E-Corona.  Thus  by  going  back  and  isolating  the  E-Corona 
component,  producing  the  equation  for  the  straight  ray  path,  minimizing  that  equation  and 
finally  simplifying  Equation  2.13  is  produced.  Equation  2.13  is  the  locus  of  focal  distance 
minima.  Figure  2.4  shows  the  simple  parabolic  nature  of  the  locus.  The  E-Corona 
focusing  conditions  state  that  the  waves  must  travel  less  than  1.58  solar  radii  and  that  is 
why  Figure  2.4  terminates  at  that  point.  With  that  bmin  is  plotted  as  a  function  of 
frequency  in  Figure  2.5  which  has  a  lower  limit  defined  by  the  E-corona  focusing 
condition  and  an  upper  limit  defined  by  the  minimum  impact  parameter  of  one. 


where 
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Figure  2.4:  Locus  of  Focal  Minima 

With  these  two  figures  the  absolute  minimum  focal  distance  is  found  to  reside  at 
-625  AU  and  at  a  frequency  of  -355  GHz  .  To  find  other  focal  distances  the  reader  may 
choose  a  desired  frequency  in  the  range  of  1 1 .07  to  455  GHz  and  use  Figure  2.5  to  find 
the  corresponding  impact  parameter.  Then  using  that  impact  parameter  go  to  Figure  2.4 
and  locate  the  corresponding  minimum  focal  distance. 

To  utilize  the  gravitational  lens  Dr.  Maccone  proposed  the  use  of  a  12  meter 
parabolic  radio  telescope,  which  will  serve  as  the  primary  payload  in  this  design  [18]. 
This  payload  will  allow  the  spacecraft  to  measure  the  variations  of  the  sun’s  focus  to 
characterize  it  and  to  take  images  of  opportunity  using  radio  telescopy.  Though  optical 
telescopes  have  been  used  to  investigate  gravitational  lensing  on  Earth  optical  telescopes 
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Figure  2.5:  Minimum  Impact  Parameter  as  a  Function  of  Frequency 


were  not  investigated  in  this  study  because  optical  frequencies,  in  the  range  of  400  -  750 
THz,  are  not  present  in  the  close  sun  approximation.  As  a  result,  their  focusing  distance 
would  be  beyond  1000  AU. 

The  resultant  gain  from  the  radio  telescope  plus  the  sun  would  be 


16  n4GMSimrl 


r  Total 


[18] 


(2.15) 


where  G  is  the  gravitational  constant,  MSm  is  the  mass  of  the  sun,  rantenna  is  the  radius  of 
the  antenna,  c  is  the  speed  of  light,  and  v  is  the  frequency  being  observed.  This  equation 
motivates  the  mission. 
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It  is  clear  from  Equation  2.15  that  increasing  the  observed  frequency  has  a  dramatic 
effect  on  the  total  gain.  However,  the  opposite  is  true  for  the  image  size.  The  image  size  is 
defined  as  the  distance  x  in  Figure  2.1  where  the  gain  falls  by  6  dB  [18].  Mathematically, 
the  image  size  is  illustrated  by  the  equation 


V ImageS  ize  ~ 


c2Vd 


[18] 


(2.16) 


In1  V GMsunv 

where  D  is  the  distance  of  the  aperture  from  the  sun  and  all  other  variables  are  the  same  as 
they  appear  in  Equation  2.15.  Though  the  frequency  makes  the  image  smaller,  as  the 
spacecraft  travels  away  from  the  sun  the  radius  of  the  image  will  get  larger  and  the 
magnification/resolution  will  also  improve. 

The  physics  of  gravitational  lensing  is  well  established  and  has  been  verified. 
Nonetheless,  the  exact  location  of  the  sun’s  gravitational  focus  is  still  merely  an 
approximation  due  to  the  calculation’s  reliance  on  computational  and  theoretical  models. 
But  what  is  known  is  that  the  focus  will  lie  somewhere  between  550  and  800  AU. 
Consequently,  that  will  be  the  primary  mission  of  the  spacecraft-determine  the  exact 
location  gravitational  lens  focus,  utilize  it,  improve  the  model  of  the  lens,  and  study 
astronomical  phenomena. 


2.3  Secondary  Missions 

Since  utilizing  the  gravitational  lens  requires  the  spacecraft  to  travel  out  to  such  great 
distances,  it  provides  the  opportunity  for  the  spacecraft  to  explore  the  interstellar  medium 
in-situ.  As  a  result,  exploring  the  interstellar  medium  and  heliosphere  constitute  important 
secondary  missions. 

NASA’s  Thousand  Astronomical  Unit  (TAU)  has  already  examined  what  would  be 
important  and  scientifically  viable  to  study  in-situ.  Based  on  the  conclusions  in  the  TAU 
NASA  study,  three  secondary  missions  were  derived. 


13 


The  first  mission  is  to  examine  dust  particles  encountered  in  the  heliosphere  and 
compare  them  to  those  in  the  interstellar  medium.  The  mass,  speed,  direction,  and 
composition  are  the  principal  characteristics  of  concern.  These  values  will  allow  for 
improved  models  of  the  distribution  of  dust  particles  and  inform  scientists  on  the  origin  of 
the  particles  (interstellar  or  from  within  the  solar  system)  [20].  Six  spacecraft  have 
previously  flown  dust  detectors  within  the  solar  system  [21].  The  most  recent  to  do  so  was 
New  Horizons.  The  Venetia  Burney  Student  Dust  Counter  (VBSDC),  which  is  shown  in 
Figure  2.6,  was  flown  on  New  Horizons  to  measure  dust  levels  in  the  inner  solar  system 
and  Kuiper  Belt.  It  can  measure  dust  particles  with  masses  from  10  12  to  10-9  kg  and  radii 
in  the  range  of  1  to  10  pint  [21].  Since  the  VBSDC  would  also  achieve  the  FSM  secondary 
mission  objectives  its  specifications  will  be  used  for  the  design.The  VBSDC  has  a  mass  of 
1.6  kg,  dimensions  of  45.72  cm  x  30.48  cm,  and  an  average  power  of  5.1  W  [21]  . 


Figure  2.6:  Venetia  Bumey  Student  Dust  Counter  [1] 
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The  second  mission  is  to  map  the  magnetic  field  as  a  function  of  heliocentric  distance 
[20].  This  will  help  physicists  create  more  accurate  models  for  the  magnetic  field  of  the 
heliosphere.  It  will  also  allow  the  magnetic  instabilities  at  the  heliopause  to  be  studied 
in-situ  [20].  Concurrently,  the  magnetic  field  mapping  outside  of  the  heliosphere,  in  the 
interstellar  medium,  will  provide  useful  information  that  will  enable  physicists  to 
construct  models  for  the  origin  and  generation  of  the  galactic  magnetic  field  [20] . 
Recently,  Juno  has  flown  magnetometers  to  map  the  magnetic  field  of  Jupiter.  These  state 
of  the  art  magnetometers  could  also  map  the  heliosphere’s  magnetic  field  as  a  function  of 
heliocentric  distance. 

Figure  2.7  displays  the  two  magnetometers  flown  on  the  Juno  mission  while  Table 
2. 1  displays  their  performance. 


Figure  2.7:  Juno’s  Magnetometer  Locations  [2] 


Table  2.1:  Juno  Magnetometer  Table  [2] 


Juno  Magnetometers  j 

Flux  Gate  Magnetometer 

Scalar  Helium  Magnetometer 

Range 

(Gauss) 

2  x  10-6  -  12 

.1—12 

Accuracy 

(%) 

.05 

.002% 

Power 

(W) 

11.3 

6.5 

Mass 

ate 

15.25 

9.08 

15 


Notice  where  the  magnetometers  are  placed  in  Figure  2.7.  In  missions  where  the 
objective  is  to  use  magnetometers  for  characterizing  the  magnetic  field,  like  Juno  and 
MagSat  (which  characterized  the  magnetic  field  of  the  Earth),  the  magnetometers  are 
always  placed  on  booms  to  avoid  registering  interference  from  the  rest  of  the  spacecraft 
bus  [22].  Also,  it  is  important  to  note  that  while  the  power  and  mass  of  the  Flux  Gate 
Magnetometer  (FGM)  are  low  they  are,  in  actuality,  even  lower  because  Juno’s  FGM 
possessed  an  advanced  star  compass  that  will  not  be  necessary  for  the  heliospheric 
mapping.  Nonetheless,  the  numbers  for  the  Flux  Gate  Magnetometer  and  the  Scalar 
Helium  Magnetometer  in  Figure  2.1  will  be  used  for  the  FSM  design. 

The  final  mission  is  to  detect  neutral  and  ionized  particles  in  the  interstellar  medium 
and  heliopause.  This  will  allow  scientists  to  determine  the  initial  energy  distribution  of  the 
interstellar  medium  as  well  as  understand  how  the  solar  system  interacts  with  the  galactic 
environment  in  the  outer  heliosphere  [20].  The  Pluto  Energetic  Particle  Spectrometer 
Science  Investigation  (PEPSSI)  on  New  Horizons  could  accomplish  this  mission.  The 
PEPSSI  can  measure  ions  in  the  range  of  1  keV  to  1  MeV  in  a  120  degree  by  12  degree 
beam  [23].  It  can  determine  the  composition  and  spectrum  of  neutral  particles,  ions,  and 
electrons  [23].  The  PEPSSI  has  a  mass  of  1.5  kg  and  a  power  requirement  of  2.5  W. 

2.4  Space  Missions 

In  1958  the  newly  formed  National  Aeronautics  and  Space  Administration,  NASA, 
set  out  on  the  ambitious  mission  of  placing  a  satellite  into  a  heliocentric  orbit.  Two 
months  after  the  Soviet  Union  accomplished  this  lofty  goal  with  Funa-1  NASA  achieved 
Earth  escape  with  Pioneer  IV.  With  their  newfound  understanding  of  how  to  achieve  a 
heliocentric  orbit  NASA  pushed  their  goals  higher.  One  of  the  results  of  this  was  Pioneer 
X  and  XI  which  would  explore  Jupiter  and  Saturn. 
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Figure  2.8:  Pioneer  X  and  XI  Configuration  [3] 

From  the  Pioneer  missions  of  the  1960s  to  the  Juno  mission  launched  in  201 1  bold 
interplanetary  probes  have  laid  the  foundation  for  the  FSM.  Due  to  the  diligent  work  of 
the  mission  engineers  invaluable  lessons  were  learned  on  how  to  develop  and  implement 
all  of  the  primary  subsystems  including  attitude  control  and  propulsion,  communications 
and  data-handling,  electrical  power,  and  environmental  control.  This  section  is  meant  to 
be  an  overview  of  current  deep  space  capabilities.  For  this  overview  I  centered  my 
analysis  on  the  most  recent  deep  space  probes  and  touched  on  other  relevant  missions 
where  appropriate. 

2.4.1  Attitude  Control  and  Propulsion.  Like  all  spacecraft  systems  the  attitude 
control  subsystem  is  driven  by  mission  requirements.  Since  every  mission  varies  the 
sensor  and  actuator  complements  vary  also. 
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Table  2.2:  Attitude  Control  History 


Spacecraft 

Mission 

Mode 

Sensor 

Actuator 

Pioneer 

Explore  interplanetary  medium 

• 

Spin  Stabilized 

• 

Sun  Sensors 

Six  4.5  N  MMH 

• 

• 

Investigate  the  asteroid  belt 
Examine  Jupiter 

• 

Star  Tracker 

Thrusters 

Voyager 

• 

ExploreJupiterand  Saturn 

• 

Three  Axis 

• 

Sun  Sensor 

Eight  MMH 

Stabilized 

• 

Star  Tracker 

Thrusters 

Galileo 

• 

Determine  the  composition  of 

• 

Dual  Spin 

• 

Star  Trackers 

Twelve  10  N 

Jupiter's  atmosphere 

• 

Gyroscope 

MMH/NiOa  Thruster 

Investigate  the  composition  of 

• 

Accelerometers 

400  N  MMH/N;04 

• 

Jovian  Satellites 

Determine  the  structure  of  the 
Jovian  Magnetosphere 

Thruster 

Ulysses 

. 

Investigate  the  solar  wind 

• 

Spin  Stabilized 

• 

Sun  Sensor 

Four  Hydrazine 

• 

Study  the  heilospheric  magnetic 

• 

CONSCAN  system 

Thrusters 

field  and  radiation 

Three  Nutation 

• 

Examine  interplanetary  gas  and 
dust 

Dampers 

Cassini 

• 

Explore  the  Saturn  system 

• 

Three  Axis 

• 

Sun  Sensors 

Reaction  Wheel 

Stabilized 

• 

Star  Trackers 

Assembly 

• 

Inertial  Reference 

16  MMH/NjOj 

Units 

thrusters 

• 

Accelerometers 

445  Nthruster- 
MMH/N2Oj 

New 

• 

Explore  the  composition  and 

• 

Spin  Stabilized 

• 

Inertial  reference 

Twelve  .8  N 

Horizons 

physical  environment  of  Pluto 

• 

Three  Axis 

Units 

MMH/N204  thrusters 

Stabilized 

• 

Star  Trackers 

Four 4.4  N 

• 

Sun  Sensors 

MMH/N2Oa  thrusters 

Juno 

• 

Examine  abundance  of  water  on 

• 

Spin  Stabilized 

• 

Sun  Sensor 

Twelve  MMH/N;0, 

Jupiter 

• 

Inertial  reference 

Thrusters 

• 

Examine  Jupiter's  core  structure 

units 

645  N  MMH/N204 

• 

Characterize  Jupiter's 

• 

Star  Trackers 

thruster 

_ 

magnetosphere 

Nutation  Damper 

Table  2.2  shows  a  complete  history  of  sensor  and  actuator  complements  for  deep 
space  probes.  Furthermore,  it  demonstrates  the  differing  approaches  to  attitude  control. 
For  example,  looking  at  the  most  recent  missions-Juno  and  New  Horizons- Juno  is  spin 
stabilized  while  New  Horizons  has  multiple  attitude  modes.  Juno’s  primary  missions 
involve  surveying  characteristics  of  Jupiter  and  since  it  will  be  in  a  Jovian  orbit  spin 
stabilization  is  adequate  for  the  duration.  On  the  other  hand,  New  Horizons  will  not  enter 
into  orbit  around  Pluto  so  it  requires  three  axis  stabilization  to  acquire  specific  targets  on 
Pluto,  and  spin  stabilization  for  other  portions  of  its  mission. 

Though  the  modes  vary,  there  is  a  commonality  between  how  stabilization  is 
maintained.  There  is  a  clear  trend  in  both  the  sensor  and  actuator  suites  of  the  deep  space 
probes  listed  in  Table  2.2.  In  regard  to  the  sensor  complement,  the  standard  approach  is  to 
use  star  trackers,  sun  sensors,  and  inertial  reference  units  to  gauge  the  attitude  of  the 
spacecraft.  The  aptly  named  star  trackers  and  sun  sensors  gauge  the  position  of  the 
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spacecraft  based  on  the  location  of  a  star  and  the  sun  respectively.  The  inertial  reference 
units  are  made  up  of  four  gyroscopes.  Three  of  the  gyroscopes  are  orthogonal  and  the 
fourth  is  skewed  so  that  it  can  compensate  for  any  gyroscope  that  malfunctions  [24] .  The 
inertial  reference  units,  sometimes  called  inertial  measurement  units,  measure  the 
spacecraft’s  orientation  as  well  as  velocity 

To  maintain  and  correct  to  the  proper  attitude  spacecraft  primarily  employ 
MMH/N2O4  thrusters.  This  is  due  to  the  high  efficiency,  moderate  cost,  and  storability  of 
the  fuel.  On  top  of  their  attitude  control  thrusters,  both  Juno  and  Cassini  boast  relatively 
high  thrust  propulsion  for  conventional  large  AV  changes  like  orbit  corrections  [25] [26]. 
While  they  have  yet  to  be  flown  on  a  deep  space  probe  ion  thrusters  are  a  burgeoning  field 
that  may  revolutionize  deep  space  flight.  As  a  result,  it  is  appropriate  to  touch  on  a  few 
outer  space  missions  that  have  used  ion  thrusters  for  propulsion. 

The  first  mission  to  utilize  an  ion  thruster  was  SERT-1  in  1964  [27].  It  validated  the 
concept  of  using  ion  propulsion  in  space.  This  led  to  the  development  and  launch  of  Deep 
Space- 1.  Deep  Space- 1  flew  the  NASA  NSTAR  thruster,  demonstrated  the  viability  of 
using  an  ion  thruster  for  interplanetary  missions,  and  paved  the  way  for  the  Dawn  mission 
[28].  The  Dawn  mission  objective  is  to  explore  Vesta  and  Ceres.  This  mission  required  an 
innovative  propulsion  system  to  supply  the  required  AV.  Dawn  used  ion  propulsion  along 
with  a  gravity  assist  to  reach  the  orbits  of  Vesta  and  Ceres.  Figure  2.9  displays  Dawn’s 
exact  orbit  and  the  characteristic  spiraling  of  ion  engines.  In  total  Dawn  will  be  thrusting 
for  1100  days  [29]. 

2.4.2  Communications  and  Data-handling.  Due  to  the  high  speed  of  advancement 
in  the  field  of  electronics  for  the  past  two  decades,  the  review  on  the  communications  and 
data  handling  systems  will  be  restricted  to  the  most  recent  spacecraft,  New  Horizons. 
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Figure  2.9:  Dawn  Trajectory  [4] 


New  Horizons  is  capable  of  storing  64  Gbits  in  non-volatile  memory  [5].  These  64 
Gbits  are  divvied  into  the  16  independently  addressable  segments  which  make  up  the  solid 
state  recorder  [5].  The  maximum  rate  that  New  Horizons  can  record  new  data  is  13 
Mbits/s  [5].  After  the  original  recording  is  complete,  data  can  processed  off  of  the  solid 
state  recorder  and  compressed  via  either  lossless  or  lossy  compression.  To  further 
compress  the  data  New  Horizons  is  also  capable  of  subframing  whereby  the  image  data 
can  be  separated  into  eight  subframes  and  compressed  individually.  Paired  with  New 
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Horizons’  software  that  can  detect  which  subframe  the  observation  target  is  on  the  other 
subframes  can  be  deleted  thus  conserving  valuable  space.  Following  compression  the  data 
can  be  transmitted  back  to  Earth  or  rewritten  to  the  solid  state  recorder.  To  help  organize 
and  control  the  flow  of  data  New  Horizons  designates  a  data  type  for  each  recording. 

There  are  51  data  types  [5].  An  example  of  the  data  types  are  compressed  science  data, 
instrument  data,  and  spacecraft  housekeeping  data  [5].  Finally,  New  Horizons  can 
bookmark  data.  This  allows  for  quick  access  to  commands  and  data. 

The  communications  subsystem  has  two  primary  tasks:  command  uplink  and 
telemetry  downlink.  To  accomplish  these  tasks  New  Horizons  employs  an  antenna 
assembly,  traveling  wave  tube  amplifiers,  an  ultra-stable  oscillator,  an  uplink  command 
receiver,  and  a  regenerative  ranging  circuit  [5]. 

The  Antenna  assembly  is  composed  of  a  forward  and  underside  antenna  system.  The 
forward  antenna  system  possesses  a  hemispherical  coverage  low  gain  antenna  (LGA),  a 
high  gain  antenna  (HGA),  and  a  medium  gain  antenna  (MGA)  [5].  The  LGA  is  capable  of 
providing  close  Earth  communication.  It  was  able  to  communicate  with  Earth  until  it  was 
1  AU  away  [5].  The  HGA  measures  2.1  meters  and  was  designed  to  meet  the  requirement 
of  a  600  bit/s  data  rate  at  its  mission  distance  [5].  The  mission  distance  is  36  AU.  The 
HGA  provides  a  42  dB  downlink  gain  when  it  is  pointed  within  .3  degrees  of  the  Earth  [5]. 
This  means  the  spacecraft  will  transmit  all  5  Gbits  of  mission  data  in  172  days  if  8  hour 
passes  are  executed  every  day  [5].  This  calculation  includes  a  2  dB  margin  [5]  .  Finally, 
the  MGA  allows  the  spacecraft  to  transmit/acquire  Earth  at  angles  up  to  4  degrees  [5].  The 
underside  antenna  system  is  comprised  of  a  single  hemispherical  coverage  low  gain 
antenna  that  serves  as  a  redundancy  for  its  twin  on  the  forward  antenna  system  [5]. 

The  ultra-stable  oscillator  (USO)  maintains  the  spacecraft’s  time  base  and  provides 
frequency  stability  for  the  uplink  and  downlink  [5].  The  USO  is  measured  to  have  an 
Allan  Deviation  of  3  x  10-13  ( unitless )  for  a  1  second  interval  and  2  x  10-13  for  an  interval 
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of  10  seconds.  Allan  Deviation  is  the  square  root  of  the  Allan  Variance  where  the  Allan 
Variance  is  “one  half  the  time  average  of  the  squares  of  the  differences  between  successive 
readings  of  the  frequency  deviation  sampled  over  the  sampling  period”  [30].  Thus,  a  low 
Allan  Deviation  is  desired  because  it  implies  good  frequency  stability  over  the  period. 

The  uplink  command  receiver  on  New  Horizons  is  innovative.  It  performs  the  tasks 
of  command  decoding,  ranging  tone  demodulation,  and  carrier  tracking  and,  unlike  its 
predecessors  utilizes  a  very  low-power  digital  design  [5].  Previously  uplink  command 
receivers  used  approximately  12  W  of  power,  but  the  uplink  command  receiver  on  New 
Horizons  consumes  only  4  W  of  power  [5].  The  8  Watts  saved  here  is  of  great  importance 
for  a  deep  space  mission  and  accounts  for  almost  5%  of  the  entire  power  budget  [5]. 

Deep  space  probes  are  tracked  using  phase  modulated  tones  sent  from  the  deep  space 
network  [5].  The  engineers  of  New  Horizons  improved  this  method  of  tracking  by 
integrating  a  regenerative  ranging  circuit  into  the  spacecraft.  The  regenerative  ranging 
circuit  works  by  using  a  delay-locked  loop  to  replicate  the  uplink  signal  and  adjust  the 
spacecraft  timing  accordingly.  This  eliminates  wide  band  uplink  noise,  a  major 
contributor  to  tracking  error,  from  the  signal  the  spacecraft  emits.  This  means  New 
Horizons  can  be  tracked  to  within  10  meters  until  it  passes  50  AU  [5]. 

The  last  pieces  of  the  communications  system  are  the  travelling  wave  tube  amplifiers 
(TWTA).  The  TWTAs  amplify  the  downlink  signal  to  provide  an  improved  data  rate. 
There  are  two  TWTAs  one  connected  to  the  MGA  and  one  connected  to  the  HGA  [5]. 

Due  the  the  presence  of  a  hybrid  coupler  both  TWTAs  can  be  connected  to  the  HGA  [5]. 
In  this  configuration  one  TWTA  will  transmit  right  hand  circular  and  the  other  will 
transmit  left  hand  circular  [5].  The  deep  space  network  architecture  has  the  capability  to 
combine  the  two  signals  on  the  ground.  Furthermore,  operating  the  TWTAs  in  this 
configuration  will  increase  the  date  rate  by  almost  2  times  as  shown  by  Figure  2.10  [5]. 
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Figure  2.10:  New  Horizons  Data  Rate  [5] 


An  important  final  note,  every  deep  space  probe  currently  communicates  with  Earth 
through  the  Deep  Space  Network  (DSN).  The  operating  frequencies  of  the  Deep  Space 
Network  are  shown  in  Table  2.3. 


Table  2.3:  DSN  Allocated  Frequency  Band  (MHz)  [14] 


Band 

Designation 

Deep  Space  Bands 
(for  spacecraft  greater  than 

2  million  km  from  Earth) 

Near  Earth  Bands 
(for  spacecraft  less  than 

2  million  km  from  Earth) 

Uplink 
(Earth  to 
space) 

Downlink 
(space  to 
Earth) 

Uplink 
(Earth  to 
space) 

Downlink 
(space  to 
Earth) 

S-band 

2110-2120 

2290-2300 

2025-2110 

2200-2290 

X-band 

7145-7190 

8400-8450 

7190-7235 

8450-8500 

K-band 

• 

25500-27000 

Ka-band 

34200-34700 

31800-32300 

• 

• 

*  No  allocation  or  not  supported  by  the  DSN. 


2.4.3  Electrical  Power.  Travelling  to  deep  space  puts  serious  restrictions  on  power 
because  it  makes  the  use  of  solar  panels  infeasible.  Even  a  theoretical,  100  percent 
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efficient  solar  panel  array  would  output  a  very  small  amount  of  power  due  to  the  rapidly 
declining  solar  irradiation.  Figure  2.1 1  clearly  demonstrates  this  phenomenon. 


1200 


1000 


E 

800 

0) 

u 

c 

(TJ 

2  600 


o 

to 


400 


200 


Distance  (AU) 

Figure  2.11:  Solar  Irradiance 


Mathematically,  the  relationship  between  the  solar  irradiance  and  distance  is  the  solar 
irradiance  is  inversely  proportional  to  the  square  of  the  distance  from  the  sun.  At  1  AU  the 
solar  irradiance  is  1358  W/m2  however  once  a  spacecraft  reaches  Jupiter,  at  approximately 
5  AU,  the  solar  irradiance  dwindles  to  a  paltry  270  W/m 2  and  again  it  drops  precipitously 
to  135  W/m2  at  10  AU.  As  a  result,  historically,  deep  space  missions  have  not  used  solar 
panels  as  their  primary  source  of  power. 
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Nonetheless,  Juno,  NASA’s  newest  deep  space  satellite,  uses  solar  panels.  Figure 
2.12  shows  the  gargantuan  size  of  Juno’s  solar  arrays.  Since  Juno’s  mission  only  involves 
traveling  to  and  orbiting  Jupiter  it  is  possible  to  use  efficient,  sizable  solar  arrays. 


Figure  2.12:  Juno  Solar  Panels  [6] 


More  specifically,  each  of  Juno’s  solar  panel  arrays  measure  23.85  nr  making  the 
total  array  size  71.55  m 2  [6].  The  panels  have  an  efficiency  of  28.3  percent  [31].  This 
provides  a  total  power  output  of  approximately  14  kW  at  Earth  and  450  W  at  Jupiter  [6]. 
However,  as  the  mission  progresses  the  solar  panel  efficiency  will  degrade  at  a  rate 
dependent  upon  the  harshness  of  the  environment.  An  electron  fluence  of  5  x  1014  at  1 
MeV  will  result  in  an  eight  percent  loss  of  beginning  of  life  power  [32].  Ultraviolet  rays 
and  impacts  (e.g.  micrometeors)  can  cause  a  1.7  %  and  1%  loss  in  beginning  of  life  power 
respectively  over  a  twenty  year  span[32]. 
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While  solar  panel  arrays  have  been  utilized,  most  of  other  deeps  space  missions  have 
opted  to  use  radioisotope  thermoelectric  generators  ,  or  RTGs  due  to  the  inefficiency  of 
solar  panels  in  deep  space.  Radioisotope  thermoelectric  generators  work  by  converting  the 
heat  produced  by  radioisotopes  into  electric  power.  Table  2.4  shows  the  progression  of  the 
RTG. 


Table  2.4:  RTG  Evolution 


RTG 

Mission  (Yr) 

BOL  Power  Output  (W) 

Power  Density  (W/kg) 

RTG  Mass  (Kg) 

Dimensions 

SNAP-3B 

Transit  4A  (1961) 

Transit  4B  (1961) 

2.7 

1.28 

2.1 

Diameter  - 12.1  cm 
Height  - 14  cm 

SNAP-9A 

Transit  5BN-1  (1963) 
Transit  5BN-2  (1963) 

26.8 

2.18 

12.3 

Diameter  -  50.8  cm 
Height  -  26.7  cm 

SNAP-19B 

Nimbus-3  (1969) 

28.2 

2.10 

13.4 

Diameter  -  53.8  cm 
Height  -  26.7  cm 

SNAP-27 

Apollo  12  (1969) 

Apollo  14  (1971) 

Apollo  15  (1971) 

Apollo  16  (1972) 

Apollo  17  (1972) 

63.5 

3.24 

19.63 

Diameter  -  40  cm 
Height  -  46  cm 

SNAP-19 

Pioneer  X  (1972) 

Pioneer  XI  (1973) 

Viking  (1975) 

40.3 

2.96 

13.6 

Diameter  -  50.8  cm 
Height  -  28.2  cm 

MHW-RTG 

Voyager  (1977) 

LES  8  (1976) 

LES  9  (1976) 

158 

4.19 

37.69 

Diameter  -  39.73  cm 
Height  -  58.31  cm 

GPHS-RTG 

Cassini  (1997) 

New  Horizons  (2006) 
Galileo  (1989) 

Ulysses  (1990) 

245.7 

4.25 

57.8 

Diameter  -  42.2  cm 
Height  - 114  cm 

MMRTG 

Mars  Science  Lab  (2011) 

120 

2.79 

43 

Diameter  -  64  cm 
Height  -  66  cm 

ASRG 

N/A 

143 

7 

20.2 

Length  -  72.5  cm 
Width -29.3  cm 
Height  -  41.0  cm 

Every  one  of  the  RTGs  listed  in  the  table  used  or  uses  plutonium-238  as  its 
radioisotope.  Plutonium-238  has  a  half-life  of  87.7  years  and  an  average  power  density  of 
570  W/kg  [33].  Figure  2.13  displays  the  exponential  decay  of  plutonium  with  time. 
Therefore,  the  initial  radioisotope  mass  must  be  chosen  so  that  it  provides  sufficient  power 
at  the  end  of  life  (EOL). 

Table  2.4  also  shows  that  great  strides  have  been  made  in  improving  the  power 
density  over  the  last  fifty  years.  The  Advanced  Stirling  Radioisotope  Generator  (ASRG) 
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Figure  2.13:  Theoretical  Decay  of  Plutonium-238 


constitutes  a  65%  increase  in  power  density  from  it  predecessor  the  General  Purpose  Heat 
Source  (GPHS).  Interestingly,  the  Multi-Mission  RTG  (MMRTG)  deviates  from  the  trend 
of  increasing  power  density.  The  goal  of  the  MMRTG,  unlike  the  ASRG,  was  not  to 
maximize  power  density  but  to  minimize  the  amount  of  plutonium-238  used  while 
maximizing  the  output  power.  In  other  words  it  was  engineered  to  maximize  the 
conversion  efficiency.  Therefore,  it  does  not  follow  the  general  trend. 

Over  the  years  deep  space  electrical  power  sources  have  evolved  from  rudimentary 
RTGs  to  highly  efficient  soiar  paneis  at  ciose  distances  and  high  power  density  RTGs  in 
the  deepest  reaches  of  space. 

2.4.4  Environmental  Control.  The  harsh  environment  of  space  (e.g.  charged 
particies,  the  vacuum  of  space,  radiation)  requires  the  spacecraft  to  be  hardened  against  a 
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multitude  of  environmental  dangers.  Yet,  out  of  all  of  the  subsystems  the  environmental 
control  subsystem  has  shown  the  least  evolution  over  time.  The  techniques  that  were 
utilized  on  the  Pioneer  missions  are  still  the  primary  methods  of  environmental  control. 
This  is  due,  in  part,  to  the  fact  that  what  was  optimal  in  the  1960s  is  still  optimal  today.  As 
a  result,  the  vacuum  of  space  presents  many  unique  challenges.  The  most  notable  among 
them  is  the  challenge  of  heat  transfer. 

In  space,  heat  can  be  transferred  away  from  the  spacecraft  only  by  radiation.  Which 
is  expressed  as, 

q  =  6<tA(T*/c  -  T4J  [34]  (2.17) 

where  €  is  the  emissivity,  cr  is  the  Stefan-Boltzmann  constant,  A  is  the  area  of  the  body, 
Ts/C  is  the  temperature  of  the  spacecraft  and  TIX  is  the  temperature  of  the  surroundings. 
However,  heat  can  be  transferred  throughout  the  spacecraft  through  all  three 
modes-radiation,  conduction,  and  convection.  So,  the  methods  of  ensuring  that  the 
spacecraft  does  not  overheat  or  freeze  are  limited. 

The  most  difficult  challenge  in  temperature  control  is  preparing  for  large  temperature 
swings.  The  temperature  instrumentation  on  the  Ulysses  will  be  examined  because  the 
spacecraft  travels  out  to  5  AU  and  close  to  the  sun  so  it  was  engineered  to  radiate  heat  as 
well  as  be  insulated  depending  on  the  varying  heat  flux.  Ulysses’  instrumentation  required 
that  the  temperature  remain  between  5  and  25  degrees  Celsius  to  guarantee  that  its  fuel 
source  did  not  freeze  and  its  electronics  did  not  overheat  [24].  To  accomplish  this  its 
engineers  used  standard  techniques. 

The  first  technique,  insuiation,  reiied  on  muiti-iayer  insuiation  biankets,  or  ML1 
biankets,  which  are  shown  in  Figure  2.14  as  the  golden  colored  covering.  These  blankets 
are  layered  to  help  insulate  and  maintain  thermal  balance.  With  that,  their  outer  layer  is 
reflective  to  decrease  the  emissivity  and  absorptivity  of  the  spacecraft.  Mathematically, 
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Figure  2.14:  Ulysses  Spacecraft  [7] 


the  easiest  way  to  understand  their  usefulness  is  to  examine  the  one-dimensional 
composite  heat  transfer  equations. 


q  =  UAAT  [34] 


(2.18) 


The  blankets  use  the  principles  of  conduction  to  limit  the  flow  of  heat  and  maintain 
an  internal  balance  in  the  spacecraft.  Equation  2.19  shows  how  the  overall  heat  transfer 
coefficient  U  effects  the  heat  transfer  rate.  In  Equation  2.21,  A  is  the  surface  area  and  A T 
is  the  temperature  difference  between  the  outside  and  inside  of  the  composite  wall.  The 
overall  heat  transfer  coefficient  is  dictated  by 


(2.19) 


and 


Rtot  =  'ZRt  [34] 


(2.20) 
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where  Rt  is  the  thermal  resistance  of  a  single  layer  of  the  blanket.  The  thermal  resistance 


equation  for  radiation  is 


Ts/C  -  T0 o 

eA°iT*c  -  Tt) 


[34] 


(2.21) 


Consequently,  there  is  direct  relationship  between  the  number  of  layers,  their  thermal 
resistances,  thermal  emissivities,  and  the  heat  transfer  rate.  An  increased  number  of 
layers,  decreased  emissivities  and  higher  thermal  resistances  correspond  to  a  lower  heat 
transfer  rate;  thus,  keeping  the  satellite  warm  or  cool.  Ulysses’  MLI  blankets  were 
composed  of  20  layers  of  aluminized  mylar  that  possess  emissivities  in  the  range  of  0.06 
to  0.08  [24].  On  top  of  the  MLI  blankets,  Ulysses  also  possesses  a  commandable  thermal 
radiator  which  radiates  heat  away  from  the  spacecraft. 

The  second  technique  for  thermal  control  is  optimizing  the  layout  of  the  spacecraft. 
Distance  provides  a  barrier  to  overheating.  Using  the  same  principles  as  the  blankets,  the 
more  distance  between  electronic  devices  the  more  surfaces  the  heat  must  be  conducted 
through.  By  spacing  the  electronics  in  this  manner,  the  entire  spacecraft  can  maintain  a 
steady,  uniform  temperature. 

Finally,  though  Ulysses  did  not  possess  them,  radioisotope  heater  units  (RHU)  have 
been  used  in  the  past  to  keep  spacecraft  warm  in  the  deepest  reaches  of  space.  Similar  to 
the  aforementioned  RTGs,  RHUs  utilize  radioisotopes,  but  they  radiate  the  heat  produced 
by  the  radioisotopes  and  do  not  convert  the  heat  to  electric  power.  As  a  result,  they  are 
much  smaller  and  very  efficient  as  a  long  term  heat  source.  Radioisotope  heater  units  also 
use  plutonium-238  [33].  They  weigh  only  39.7  grams  and  measure  2.54  cm  in  diameter 
and  3.3  cm  in  height  [33].  They  output  1  W  of  heat  and  when  these  units  are  placed 
properly  they  can  adequately  heat  a  small  area  [33]. 

Charged  particles  present  another  unique  challenge.  As  spacecraft  travel  on  their 
journey  they  are  bombarded  with  charged  particles  which  causes  the  phenomenon  known 
as  spacecraft  charging.  Understanding  Space  defines  spacecraft  charging  as  a  build  up  of 
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charge  on  different  parts  of  a  spacecraft  [35].  What  is  more,  this  build  up  of  charge  is 
dangerous  because  it  may  eventually  result  in  a  discharge,  which  is  the  natural  result  when 
there  is  a  concentration  of  charge  in  one  place  and  an  absence  in  another.  These 
discharges  can  damage  the  electronics,  solar  panes,  and  surface  coatings  [35].  To  combat 
this,  spacecraft  are  covered  with  highly  conductive  materials  like  Indium  Tin  Oxide  [24]. 
This  conductive  coating  provides  a  passive  defense  against  charging  by  distributing  any 
charge  build  up  so  that  a  discharge  will  not  occur. 

Charged  particle  collisions  also  result  in  sputtering.  To  borrow  the  analogy  from 
Sellers’  text  sputtering  is  like  “sand  blasting”  the  spacecraft  [35].  A  multitude  of  small 
collisions  between  charged  particles  and  the  spacecraft  cause  the  spacecraft’s  exterior  to 
ablate.  The  materials  ablated  from  the  surface  often  return  due  to  lighter  pressure  from 
radiative  heat  transfer  to  another-usually  cold-part  of  the  spacecraft  and  cause  damage. 
This  in  turn  results  in  lower  performance  from  the  spacecraft’s  external  sensors  and  a  slow 
degradation  in  efficiency  of  the  spacecraft’s  solar  panels  [35].  To  overcome  this, 
astronautical  engineers  utilize  specific  space  qualified  materials  when  manufacturing 
parts.  Although  the  controlled  manufacturing  processes  do  not  eliminate  sputtering,  they 
do  minimize  it.  These  specialized  parts  are  also  selected  to  prevent  out-gassing  and  cold 
welding. 

Out-gassing  is  the  expulsion  of  gasses  from  material.  These  gasses,  e.g.  H20,  are 
trapped  in  the  material  on  Earth  because  of  the  pressure  exerted  by  the  atmosphere  [35]. 
However,  once  these  materials  are  introduced  into  the  vacuum  environment  of  space  the 
gasses  are  expelled.  This  expulsion  can  cause  arcing  and  damage  the  electronics  and 
sensor  payload,  and  is  again  combated  by  using  space  qualified  parts  [35]. 

Cold  welding  is  when  two  parts  on  the  spacecraft  fuse  together  due  to  having  little 
separation  between  them.  On  Earth  there  may  be  a  free  flow  of  air  molecules  between  the 
two  parts  to  prevent  this  phenomenon  but  once  in  space  the  air  flow  stops  and  the  parts  are 
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welded  together  [35].  Cold  welding  is  combated  by  using  material  combinations.  Two 
different  materials  when  placed  close  together  are  less  likely  to  fuse  together  than  two 
similar  materials. 

The  final  major  challenge  of  the  space  environment  is  radiation.  In  deep  space,  the 
primary  source  of  radiation  will  come  from  galactic  cosmic  radiation.  Galactic  cosmic 
radiation  (GCR)  is  made  up  of  high  energy  nuclei  in  the  MeV  to  GeV  range  [36].  The 
maximum  yearly  exposure  to  GCR  for  a  spacecraft  is  approximately  1.3  x  108 
protons  I  cm2  which  is  equivalent  to  10  rads/yr  [36].  The  high  kinetic  energy  of  the  GCR 
can  cause  errors  and  permanent  damage  to  solid  state  electronics  [36].  To  guard  the 
electronics  from  the  GCR  electronics  are  space  hardened  or  put  in  an  electronics  vault  like 
on  the  Juno  mission,  which  is  made  of  tantalum  [37].  The  electronics  vault  on  the  Juno 
mission,  which  expects  a  radiation  dose  greater  than  100  Mrad,  decreases  the  exposure  the 
electronics  see  by  almost  four  orders  of  magnitude  to  25  Krad  [37]. 

2.5  Conclusion 

This  chapter  presented  what  has  been  done  thus  far  in  satellite  systems  and 
introduced  the  physics  of  the  gravitational  lens.  By  understanding  what  has  been  done  in 
satellite  systems  an  approach  was  developed  for  the  design  of  the  FSM  spacecraft,  which 
is  to  follow  in  Chapter  3. 
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3  Analysis 


3.1  Introduction 

After  extensive  study  on  the  concept  of  gravitational  lensing  between  1987  and  1993, 
Dr.  Maccone,  the  originator  of  the  idea,  submitted  a  medium-class  mission  (M3)  proposal 
to  the  European  Space  Agency  for  the  design  of  the  Focal  Space  Mission.  This  proposal 
addressed  a  number  of  things  ranging  from  the  impetus  for  the  mission  to  the  payload 
requirements.  Following  his  proposal,  Dr.  Maccone  continued  to  expand  his  ideas,  and  in 
Deep  Space  Flight  and  Communications  he  compiled  the  most  recent  information  on  the 
topic.  In  the  book,  he  touches  on  the  physics  of  gravitational  lensing,  its  enabling 
capabilities,  and  the  myriad  of  ways  to  exit  the  solar  system.  However,  what  is  not 
discussed  in  real  terms  are  the  communication  challenges,  the  specific  subsystems 
requirements,  the  subsystems  themselves,  an  in-depth  AV  analysis,  trajectory  analysis, 
and  possible  launch  vehicles.  Dr.  Maccone  presents  a  very  thorough  scientific  analysis  but 
does  not  delve  into  the  engineering  required  and  that  is  the  heart  of  the  next  two  chapters. 

3.2  Launch  Vehicle 

The  launch  vehicle  analysis  for  this  mission  comes  down  to  one  thing-C3.  C3  is  the 
term  used  to  describe  the  excess  energy  provided  by  the  launch  vehicle.  It  is  defined  by 
the  equation 

C3  =  V^  [15]  (3.1) 

where  VA  is  the  hyperbolic  excess  velocity.  The  hyperbolic  excess  velocity  is  the  velocity 
difference  between  the  escape  velocity  and  the  velocity  of  the  spacecraft.  The  Earth 
escape  velocity  for  a  circular  orbit  is  given  by 


where  G  is  the  gravitational  constant,  ME  is  the  mass  of  the  Earth,  and  r  is  the  orbit  radius. 
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When  the  same  launch  vehicle  configuration  is  considered  a  higher  C3  implies  a 
higher  initial  mission  velocity  but  a  lower  delivered  mass.  The  remainder  of  this  section 
will  explore  the  optimization  of  A V  as  a  function  of  C3.  In  more  specific  terms,  the  goal  is 
to  find  the  launch  C3  that  maximizes  the  mission  AV. 

The  first  step  in  this  examination  was  to  look  at  the  ideal  rocket  equation 

AV  =  Ispgo  ln|~j  P8]  (3.3) 

where  Isp  represents  the  specific  impulse,  go  is  the  acceleration  due  to  gravity  at  the 
Earth’s  surface,  m0  is  the  initial  mass,  and  nif  is  the  final  spacecraft  mass.  The  initial  mass 
is  equal  to  the  sum  of  the  propellant  mass  mp  and  the  final  spacecraft  mass.  Thus  clearly 
an  increase  in  propellant  mass  with  no  change  to  the  final  mass  will  result  in  a  logarithmic 
increase  in  AV.  Also,  an  increase  in  specific  impulse  will  result  in  a  linear  increase  in  AV . 
With  this  knowledge  it  is  now  appropriate  to  look  at  trends  in  the  excess  energy. 

Quickly  looking  at  Figure  3.1  the  trend  is  unmistakable.  As  mentioned  earlier,  for  a 
given  launch  vehicle  configuration  the  larger  the  payload  mass  the  lower  the  excess 
energy.  Therefore,  if  the  final  spacecraft  mass  is  assumed  to  be  the  minimum  mass 
required  to  accomplish  the  mission  then  the  only  mass  that  is  varying  is  that  of  the 
propellant.  So  an  increase  in  propellant  mass  corresponds  to  a  decrease  in  excess  energy. 

Exploring  the  concept  in  more  depth,  the  launch  vehicle  configuration,  similar  to  the 
spacecraft,  operates  under  the  constraints  of  the  ideal  rocket  equation.  This  means  the  AV 
of  the  launch  vehicle  is  also  a  function  of  the  specific  impulses,  final  masses  and 
propellant  masses  of  the  launch  vehicle  configuration.  Therefore,  the  optimization  of  the 
AV  as  a  function  of  the  C3  is  dependent  only  upon  these  values  for  the  launch  vehicle  and 
the  spacecraft.  With  that,  those  values  are  constants  for  a  given  launch  vehicle 
configuration  and  are  what  produce  the  curves  in  Figure  3.1.  As  a  result,  the  optimization 
is  constrained  to  the  curve  meaning  it  is  now  dependent  solely  on  the  spacecraft’s  values 
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Figure  3.1:  Selected  Launch  Vehicle  Performance  [8] 


for  specific  impulse,  propellant  mass,  and  final  mass.  Interestingly,  the  spacecraft’s 
propellant  mass  and  final  mass  are  also  constrained  by  the  launch  vehicle  curves  because 
the  curves  specify  the  delivered  mass,  which  is  the  sum  of  the  propellant  mass  and  final 
mass.  This  leaves  the  specific  impulse  of  the  spacecraft  as  the  only  free  variable. 

Figure  3.2  shows  the  scaling  effect  that  the  specific  impulse  has  on  the  AV.  An 

increased  specific  impulse  results  in  a  linear  increase  in  the  AV  for  a  given  mass  ratio 
which  is  defined  as 


(3.4) 
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Figure  3.2:  Effect  of  Specific  Impulse  on  AE 

Due  to  this  fact,  if  the  spacecraft  has  a  higher  specific  impulse  than  the  launch  vehicle 
upper  stage  then  optimizing  the  AE  by  varying  the  C 3  requires  the  C3  to  be  minimized 
because  the  spacecraft  will  use  the  mass  more  efficiently.  Upper  stages  have  low  specific 
impulses  in  the  hundreds,  e.g  the  Star  48  has  a  specific  impulse  of  289.9  seconds, 
compared  to  electric  propulsion  options  with  thousands  of  seconds  of  specific  impulse 
[39]  .  Hence,  the  presence  of  an  electric  thruster,  e.g.  the  NSTAR  thruster  mentioned  in 
the  literature  review,  in  the  range  of  thousands  of  seconds  would  imply  that  to  optimize 
launch  AE  ,  C3  should  be  minimized.  However,  if  a  thruster  is  not  present  or  it  will  only 
be  used  in  orbit  corrections/gravity  assists  the  opposite  is  true.  The  optimum  C3  is  the 
maximum  C3. 
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The  optimum  launch  vehicle  for  the  first  case  (electric  thrusters  are  present)  will  be 
the  one  that  provides  the  maximum  delivered  mass  to  the  minimum  C3.  In  regard  to  this 
mission,  the  minimum  C 3  is  the  parabolic  excess  energy  of  0  km2 Is1.  This  is  taken  as  the 
minimum  to  avoid  spiraling  around  the  Earth.  For  the  second  case  (thrusters  only  used  for 
orbit  corrections  and/or  gravity  assist),  the  optimum  launch  vehicle  would  be  the  one  that 
delivers  the  spacecraft  mission  mass  to  the  highest  C3.  So  based  on  these  principles  after  a 
propulsion  system  is  chosen  a  launch  vehicle  may  be  chosen. 


3.3  Communications  and  Data  Handling 


The  communications  and  data  handling  subsystem  gathers  data  from  the  spacecraft, 
processes  the  data  and  transmits  the  data  to  the  ground.  The  communications  and  data 
handling  subsystem  also  receives  commands  from  the  ground  station.  The  transmission 
and  receipt  of  data  is  analyzed  through  the  link  budget. 

The  link  budget  is  the  relationship  between  data  rate,  transmission  power, 
transmission  path  and  antenna  size.  It  is  best  illustrated  by  the  energy  per  bit  to  noise 
power  spectral  density  ratio,  which  is  given  in  the  following  equation  where  P  is  the 
transmitter  power,  L/  is  the  transmitter-to-antenna  line  loss,  G,  is  the  transmit  antenna 
gain,  Ls  is  the  space  loss,  La  is  the  transmission  path  loss,  G,  is  the  receiver  antenna  gain, 
k  is  Boltzmann’s  constant,  Ts  is  the  system  noise  temperature,  and  R  is  the  data  rate  in  bits 
per  second. 


Eb  _  PE]G,LsLaGr 

WG  ~  Wji 


(3.5) 


The  transmitter  power  is  an  input  parameter  with  units  of  W.  The  amount  of  transmitter 
power  available  is  dependent  upon  the  choice  of  the  power  system.  However,  the  power 
system  has  yet  to  be  determined  at  this  point  in  the  analysis.  As  a  result,  a  range  of  1  to 
300  W  was  examined  throughout  this  section. 
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The  transmitter-to-antenna  line  loss  represents  the  communication  losses  that  occur 
when  data  is  transferred  from  transmitter  on  the  spacecraft  to  the  spacecraft  antenna.  On 
average,  this  loss  is  .8,  a  unitless  value  [15]. 

The  transmitter  and  receiver  antenna  gains  are  given  by  the  equation 


(3.6) 


where  tt  is  the  familiar  mathematical  constant,  D  is  the  diameter  of  the  antenna,  77  is  the 
antenna  efficiency,  /  is  the  frequency,  and  c  is  the  speed  of  light.  Thus  an  increase  in 
antenna  size  will  lead  to  a  parabolic  increase  in  gain,  as  will  an  increase  in  frequency.  The 
antenna  efficiency  is  a  function  of  the  antenna  shape.  Three  of  the  most  common  antenna 
types  are  the  paraboiic  reflector,  helix,  and  horn  which  have  efficiencies  of  .55,  .7,  and  .52 
respectiveiy.  Finaiiy,  the  frequency  of  the  transmission  is  dependent  upon  the  channeis 
avaiiabie  in  the  DSN.  Each  channel  has  a  bandwidth  of  approximately  370  kHz  and  falls 
within  the  frequency  ranges  shown  in  Table  2.3.  Though  the  size  of  the  spacecraft 
antennas  have  yet  to  be  determined  in  this  analysis  the  DSN  receiver  antennas  are 
parabolic  and  70  meters  in  diameter  [14]. 

For  deep  space  missions  the  space  loss  has  a  gargantuan  effect  on  the  signal  to  noise 
ratio.  Given  by  the  equation 


the  space  loss  represents  the  attenuation  of  the  transmitted  signal  due  to  travelling  through 
free  space.  In  Equation  3.7,  c  is  again  the  speed  of  light  in  meters  per  second,  n  is  the 
familiar  constant,  S  is  the  distance  from  the  antenna  to  the  ground  receiver  and  /  is  the 
frequency  of  the  signal.  Therefore,  the  further  the  spacecraft  travels  from  Earth  the  greater 
the  loss  as  shown  in  Figure  3.3  .  Notice  that  Figure  3.3  is  plotted  on  a  logarithmic  scale  in 
the  vertical  axis.  Finally,  for  the  figure,  the  frequency  was  held  constant  at  2291.66  MHz, 
one  of  the  DSN  frequencies.  Furthermore,  this  will  be  the  assumed  FSM  frequency. 
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Figure  3.3:  Space  Loss 

The  transmission  path  loss  is  variable,  and  it  is  based  on  the  conditions  of  the  signal 
medium.  Space  Mission  Analysis  and  Design  by  Larson  and  Wertz  outlines  the  sources  of 
transmission  path  loss:  man-made  noise,  solar  noise,  galactic  noise,  atmospheric  noise, 
quiet  sun,  sky  noise,  and  black  body  radiation  [15]. The  sources  that  have  the  greatest 
affect  on  the  assumed  frequency  of  2291.66  MHz  are  black  body,  galactic  noise,  and  quiet 
sun.  The  total  effect  of  all  sources  results  in  a  transmission  path  loss  of  ten  percent 
meaning  La  is  equal  to  .9. 

The  data  rate,  which  has  units  of  bits  per  second,  is  the  rate  information  is  transferred 
over  a  communications  link.  This  rate  is  determined  based  on  the  modulation  which  will 
be  discussed  later  in  this  section. 
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The  final  piece  of  the  Eb/N0  ratio  is  the  noise  spectral  density.  The  noise  spectral 
density  is  the  Boltzmann  Constant,  1.38  x  10_23!^,  multiplied  by  the  system  noise 
temperature.  The  system  noise  temperature  is  dependent  upon  the  frequency.  For  the  DSN 
ranges  the  uplink  and  downlink  system  noise  temperatures  are  614  K  and  135  K 
respectively. 


3.3.1  Shannon  Limit.  In  the  Eb/Na  ratio  the  data  rate  is  limited  by  the 
Shannon-Hartley  Theorem.  The  theorem  specifies  that  there  exists  a  maximum  rate  of 
information  transfer  for  a  given  bandwidth,  which  is  dictated  by  the  expression 


Rmax  ~  B  lo§2  I  1  H" 


C ElRP)LsLaG , 
kTsB 


[15] 


(3.8) 


where  the  E1RP  is  the  effective  isotropic  radiated  power  in  watts;  Ls,  La,  k,  Ts,  and  G,  are 
the  same  as  in  Equation  3.5;  and  B  is  the  bandwidth.  The  Shannon  Limit  is  being  explored 
because,  with  the  enormous  space  loss  term  at  the  mission  distance,  it  is  important  to 
ensure  that  the  maximum  data  rate  is  large  enough  that  data  can  be  transferred  in  a 
realistic  amount  of  time.  If  it  is  not,  then  the  mission  is  not  viable. 

To  calculate  the  maximum  data  rate  an  E1RP  of  10000  W  was  assumed.This 
corresponds  to  an  approximate  antenna  diameter  of  3  m,  transmitter  power  of  4  W, 
frequency  of  2291.66  MHz,  transmitter  efficiency  of  .55  and  line  loss  of  .8.  The  equation 
for  the  E1RP  is 

EIRP  =  PLjG,  [15]  (3.9) 


where  the  definition  of  the  variables  have  been  previously  provided.  Also,  since  the 
relevant  case  is  the  worst  case  the  mission  distance  was  taken  to  be  1000  AU,  the  EOL 
distance.  The  Shannon  Limit  data  rate  at  1000  AU  given  the  values  stipulated  earlier  in 
this  section,  a  bandwidth  of  100  kHz  and  using  Equation  3.8  is  44.65  kB/ sec.  This  rate  is 
sufficient.  The  EIRP  and  mission  distance  were  varied  in  Figure  3.4  to  produce  a  more 
detailed  picture  of  the  Shannon  Limit .  With  that,  the  figure  displays  the  positive 
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Figure  3.4:  Shannon  Limit 


correlation  between  EIRP  and  the  Shannon  Limit  as  well  as  the  negative  correlation 
between  the  mission  distance  and  the  Shannon  Limit.  Moreover,  Figure  3.5  shows  the 
constant  EIRP  curves.  These  curves  demonstrate  how  each  EIRP  in  Figure  3.4  can  be 
achieved  by  varying  the  transmitter  power  and  antenna  diameter.  While  the  graph 
describes  configurations  that  would  be  undesirable  such  as  requiring  197  W  for  a  1 -meter 
antenna  to  achieve  an  EIRP  of  50000  W,  it  is  also  provides  reasonable  options  such  25  W 
to  an  8-meter  antenna  to  produce  the  same  EIRP. 

Concluding  this  subsection  on  the  Shannon  Limit,  it  is  important  to  note  that  the 
Shannon  Limit  is  a  theoretical  limit  that  can  be  approached  but  not  achieved.  It  is 
approached  through  the  use  of  modulation  and  coding  schemes.  The  modulation  and 
coding  scheme  chosen  dictates  the  actual  data  rate  and  bandwidth. 
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Figure  3.5:  Constant  EIRP  Curves  (W) 


3.3.2  Signal  Modulation.  Modulation  allows  information  to  be  conveyed  through 
radio  frequencies.  The  four  types  of  modulation  are  polarization,  frequency,  phase,  and 
amplitude  [15].  Spacecraft  primarily  use  frequency  and  phase  modulation  due  to 
improved  power  efficiency  [15].  Table  3.1  shows  the  most  common  modulation  schemes. 

The  second  column  of  Table  3.1  displays  the  Eb/N0  required  for  each  modulation  scheme 

to  achieve  a  bit  error  rate  (BER)  of  10-5  .  By  definition  the  BER  ’’gives  the  probability  of 
receiving  and  erroneous  bit”  [15].  In  regard  to  Table  3.1,  a  BER  of  10~5  implies  that  the 
data  transmitted  will  have  only  one  error  in  every  105  bits.  This  is  an  acceptable  BER. 

The  best  modulation  scheme  for  this  mission  is  the  one  that  provides  the  best  BER 
performance.  The  best  performance  is  defined  as  the  one  which  minimizes  the  Eb/Na. 
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Table  3.1:  Modulation  Schemes 


Modulation 

Eb/Nofor  BER=  10'5  (dB) 

BPSK 

9.6 

DPSK 

10.3 

QPSK 

9.6 

FSK 

13.3 

8FSK 

9.2 

BPSK  and  QPSK  Plus 

R-l/2  Viterbi  Decoding 

4.4 

BPSK  and  Plus  RS  Viterbi 
Decoding 

2.7 

8FSK  Plus  R-l/2  Viterbi 
Decoding 

4.0 

MSK 

9.6 

Though  not  listed  in  the  table,  the  Shannon  Limit/minimum  Eb/Na  for  a  BER  of  10-5  is 
-1.5  dB  [15].  This  means  that  the  minimum  available  Eb/N0  of  2.7  dB  is  only  4.2  dB 
greater  than  the  Shannon  Limit. 

The  binary  phase  shift  keying  (BPSK)  and  plus  Reed-Solomon  (RS)  Viterbi 
Decoding  technique  works  by  concatenating  modulation  and  coding.  First,  BPSK 
modulates  the  phase  of  the  transmitted  signal  between  two  settings-0  and  180  degrees 
[15].  These  two  phases  represent  0  and  1  respectively  [15].  Next,  the  signal  is  encoded 
using  forward  error  correction.  Forward  error  correction  is  implemented  to  ’’reduce  the 
Eb/Na  requirement”  [15].  In  this  case  the  method  of  forward  error  correction  is 
convolutional  coding  with  RS  Viterbi  decoding.  The  rate  of  convolutional  coding  is  1/2, 
meaning  that  for  every  data  bit  two  bits  are  transmitted  [15].  Finally,  when  using  BPSK 
modulation  the  bandwidth  is  equal  in  magnitude  to  the  data  rate  [15]. 

Rearranging  Equation  3.5  produces  Equation  3.10,  which  is  an  equation  for  the  data 
rate  as  a  function  of  EIRP  and  distance  since  all  of  the  other  variables  in  the  equation  are 
constants. 

(. EIRP)LsLaGr 

R  =  - - ^  (3.10 

Eb 
kTs  — 
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The  Eb/Na  used  in  calculating  the  data  rate  is  a  unitless  1.862  which  is  produced  after 
converting  2.7  dB  into  its  non-dimensional  variant.  A  graph  of  Equation  3.10  for  a  range 
of  EIRP  and  distances  is  shown  in  Figure  3.6  to  give  a  clearer  picture  of  the  relationship. 
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Figure  3.6:  Mission  Data  Rate 


In  conclusion,  based  on  the  size,  type  and  power  of  the  antennas  flown  the  mission 
data  rate  can  be  found  using  Equation  3.10.  The  size,  type  and  power  of  the  antennas  also 
affect  other  communications  aspects  which  will  be  explored  in  the  following  subsection. 

3.3.3  Antennas.  In  the  Background  and  Theory  the  communications  architecture 
of  previous  spacecraft  were  discussed.  Most  spacecraft  had  multiple  antennas  of  different 
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types.  Spacecraft  antennas  fall  into  three  categories  high  gain  antennas  (HGA),  medium 
gain  antennas  (MGA),  and  low  gain  antennas  (LGA).  High  gain  antennas  possess  high 
data  rates  and  narrow  half-power  beamwidths  (HPBW).  On  the  other  hand,  low  gain 
antennas  have  wide  HPBW  but  low  data  rates.  Finally,  the  medium  gain  antennas  make  up 
the  space  between  the  extremes.  Each  antenna  type  is  suited  for  a  different  role.  The  HGA 
is  suited  to  high  bit  rate  transmissions  of  mission  data.  Low  gain  antennas  are  primarily 
used  for  transmitting  data  at  near  Earth  distances.  On  the  other  hand,  medium  gain 
antennas  are  used  to  acquire  Earth  and  initiate  the  link  between  the  spacecraft  and  the 
ground  station. 

Equation  3.6  displayed  how  the  gain  is  calculated  for  a  parabolic  antenna.  Hence,  to 
manipulate  the  gain  the  engineer  must  manipulate  the  diameter  of  the  antenna.  However, 
changing  the  diameter  also  affects  the  HPBW.  The  HPBW  is  dictated  by  the  expression 

HPBW  =  (3.11) 

JGHz D 

where  /  is  the  signal  frequency,  D  is  the  antenna  diameter  and  the  HPBW  is  measured  in 
degrees.  The  HPBW  is  defined  as  the  ’’angle  across  which  the  gain  is  50  percent  of  the 
peak  gain”  [15].  From  Equation  3.11  and  Equation  3.6  it  is  clear  that  the  HPBW  and  gain 
are  inversely  related.  Therefore,  the  higher  the  gain  the  more  narrow  the  beamwidth  and 
vice  versa. 

The  data  rate,  HPBW,  and  gain  will  be  the  driving  factors  in  the  design  of  the 
spacecraft  antennas,  but  the  antennas  are  just  one  piece  of  the  communications  and  data 
handling  subsystem. 

3.3.4  Autonomy  and  Lifespan.  Radio  frequency  communications  travel  at  the 
speed  of  light  so  latency  becomes  a  problem  with  increasing  distances.  At  550  AU  it 
would  take  approximately  3.16  days  for  a  signal  from  the  spacecraft  to  reach  the  Earth.  As 
a  result,  the  spacecraft  will  require  a  high  level  of  autonomy  and  a  lot  of  storage  space. 
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Autonomy  requires  a  very  capable  processor  .  An  example  processor  is  the  BAE 
RAD750.  Flown  on  the  Juno  and  Curiosity  missions,  it  is  capable  of  operating  at  133 
MHz  and  withstanding  1  Mrad  of  radiation.  What  is  more,  the  RAD750  has  490  year 
mean  time  between  failures  [40].  Finally,  it  has  a  mass  of  only  9  g  and  a  power  dissipation 
of  5  Wat  133  MHz  [40]. 

To  store  the  flight  data  the  spacecraft  will  use  the  same  proven  set-up  as  New 
Horizons,  which  was  outlined  in  the  Background  and  Theory.  New  Horizons  utilized  a  64 
Gbit  solid  state  recorder  with  recording  capabilities  of  13  Mbits/sec.  Unlike  New 
Horizons  the  FSM  will  use  a  radio  telescope  instead  of  high  resolution  optical  telescopes; 
thus,  its  data  requirements  will  be  lower.  That  was  the  reasoning  behind  this  course  of 
action.  Consequently,  the  mission  data  rate  requirement  will  be  500  bits/sec  based  on  the 
payload  data  rates.  This  will  allow  5  Gbits  of  mission  data  to  be  downloaded  in 
approximately  240  days  with  12  hour  passes. 

3.4  Power 

There  are  two  primary  requirements  in  powering  the  FSM  that  will  dictate  what 
power  systems  may  be  used.  The  first  requirement  is  to  provide  sufficient  power  to  operate 
the  mission  at  the  EOF.  The  second  requirement  is  to  provide  enough  power  to  propel  the 
spacecraft  during  the  boost  phase  of  the  mission. 

At  the  EOF  the  spacecraft  will  be  at  least  800  AU  from  the  sun;  thus,  solar  power 
will  not  be  an  option.  That  leaves  the  options  of  nuclear  thermal  power  and  radioisotope 
electric  power.  Nuclear  thermal  can  be  eliminated  from  contention  immediately  because  it 
does  not  meet  the  COTS  requirement. 

The  three  major  RTGs  are  the  MMRTG,  GPHS-RTG,  and  the  ASRG.  The  first  two 
were  briefly  mentioned  in  the  Background  and  Theory  and  have  previously  been  flown. 
The  ASRG  is  currently  in  development  at  NASA  Glenn  Research  Center.  Its  TRF,  or 
Technology  Readiness  Fevel,  is  TRF  6.  A  TRF  of  6  means  that  a  prototype  has  been 
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tested  to  high  fidelity  in  a  laboratory.  Thus  the  ASRG  is  undergoing  lifetime  and 
performance  testing  and  does  not  classify  as  COTS.  Table  2.4  compares  the  two  remaining 
power  options  based  on  mass,  dimensions,  and  power  density. 

Both  of  the  power  sources  utilize  Pu-238  which  has  a  half-life  of  87.7  years.  As  a 
result,  the  mission  duration  should  account  for  this.  Meaning,  if  possible,  the  spacecraft 
should  arrive  and  have  transmitted  a  large  amount  of  data  before  this  time  is  reached.  The 
degradation  in  power  will  also  result  in  lower  data  rates  and  limited  spacecraft  operation. 

Finally,  the  second  requirement,  which  is  to  provide  power  for  the  boost  phase  of  the 
spacecraft,  does  not  preclude  the  use  of  a  solar  panel  array  along  with  the  RTG. 

3.5  Thruster 

3.5.1  Introduction.  The  thruster  analysis  is  split  into  four  parts.  The  first  part  of 
the  analysis  examines  the  spacecraft  AV  and  how  it  can  be  optimized  to  produce  the 
minimum  trip  time  .  The  second  piece  of  the  analysis  of  the  payload  mass  fraction.  After 
that,  the  analysis  addresses  thruster  efficiency  and  finally  gravity  assists. 

3.5.2  AC  The  AV,  or  velocity  change  required  to  reach  a  certain  target  in  a  desired 
amount  of  time,  analysis  is  an  integral  part  of  the  spacecraft  thruster  analysis  because  it 
determines  the  level  of  performance  required  by  the  thruster. 

The  FSM  AV  analysis  is  built  on  the  assumption  that  the  satellite  is  initially  placed  in 
a  270  km  circular  orbit,  and  that  corresponds  to  an  orbital  velocity  of  7.74  km/ s.  This 
assumption  is  based  on  the  fact  that  placing  a  satellite  in  an  orbit  of  270  km  is  well  within 
the  capabilities  of  modem  launch  vehicles.  From  there,  the  satellite  is  impulsively  boosted 
into  an  Earth  escape  trajectory  by  its  upper  stage.  At  270  km,  the  minimum  velocity  to 
achieve  Earth  escape  is  10.95  km/s,  which  can  be  found  using  Equation  3.2.  Thus  a  AV, 
of  3.21  km/s  is  required.  Furthermore,  since  these  numbers  are  based  off  of  an  assumption 
it  is  important  to  bound  its  effect.  If  the  assumed  altitude  was  the  minimum  acceptable 
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low  earth  orbit  altitude  of  200  km  then  the  AV7  would  be  3.22  km/ s  and  if  the  insertion 
orbit  was  400  km  the  AV7  would  be  3.17  km/ sec.  So  these  values  give  credence  to  the 
assumption  by  demonstrating  that  only  a  small  deviation  occurs  due  to  the  assumed 
altitude  of  270  km. 

This  new  velocity,  10.95  km/ s,  corresponds  to  a  specific  mechanical  energy,  e,  of 
119.9  km2 / s2.  Another  one  of  the  assumptions  of  this  analysis,  and  a  general  assumption 
in  AV  calculations,  was  that  the  specific  mechanical  energy,  e,  of  an  orbit  remains 
constant.  As  a  result,  it  is  easy  to  calculate  the  spacecraft’s  velocity  at  the  termination  of 
Earth’s  sphere  of  influence,  whose  location  is  determined  by  Equation  3.12,  by  plugging 
rSoi  in  for  r  in  Equation  3.13  and  solving  for  v. 


rsoi  =  aplanet(r^^\  [35]  (3.12) 

\  mSun  ) 

2 

e=\~-  [15]  (3.13) 

2  r 

If,  like  postulated,  the  spacecraft  left  at  exactly  the  escape  velocity  it  would  arrive  at 
the  sphere  of  influence  with  a  parabolic  trajectory  and  velocity  of  0  km/ s  relative  to  the 
Earth.  However,  it  would  still  possesses  a  velocity  of  29.78  km/ s  relative  to  the  sun.  Not 
surprisingly,  this  is  because  29.78  km/ s  is  the  orbital  velocity  of  the  Earth  relative  to  the 
sun.  Hence,  for  all  intents  and  purposes  the  spacecraft  is  assumed  to  be  in  Earth’s  exact 
orbit.  This  is  a  good  simplifying  assumption  since  the  Earth’s  sphere  of  influence  is 
minuscule  relative  to  1  AU.  More  specifically,  Earth’s  SOI  is  only  .006%  of  1  AU 

As  discussed  in  the  launch  vehicle  section,  it  is  possible  for  the  spacecraft  to  arrive  at 
the  sphere  of  influence  with  a  velocity  higher  than  0  km/ s  relative  to  the  Earth.  This 
velocity  is  known  as  the  hyperbolic  excess  velocity,  14, .  This  requires  a  correspondingly 
higher  AV.  This  new  AV  can  be  found  with  the  following  equation. 
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(3.14) 


AV  =  ^119.05 +  V£  -  J ^ 

In  equation  3.14,  Vo  represents  the  desired  hyperbolic  excess  velocity,  r  represents 
the  initial  circular  orbit  radius  and  /uE  represents  the  gravitational  parameter  of  Earth. 
Figure  3.7  shows  the  AV  required  to  achieve  an  excess  velocity  in  the  range  from  0  km/s 
to  20  km/s  for  an  initial  orbit  radius  of  270  km. 

An  important  excess  velocity  to  look  at  is  11.99  km/s,  or  41.77  km/s  relative  to  the 
Sun.  This  is  the  escape  velocity  for  the  Sun  at  Earth,  which  can  be  found  by  substituting 
the  mass  of  the  Sun  for  the  mass  of  the  Earth  and  the  radius  from  the  the  Sun  to  the  Earth 
in  Equation  3.2.  The  resulting  AV  to  achieve  solar  escape  is  8.47  km/s.  So,  with  this  AV 
the  spacecraft  would  have  enough  velocity  to  escape  the  sun’s  gravity;  thus,  it  could  coast 
all  the  way  to  the  mission  distance  without  a  thruster.  Finally,  with  these  numbers  and 
equations  in  hand  the  author  set  up  a  MATLAB  script  that  would  allow  for  the  analysis  of 
a  AV  trade  space. 


0  5  10  15  20 

Excess  Velocity  (km/sec) 


Figure  3.7:  Delta-V  Required  for  Excess  Velocity 

The  goal  of  the  aforementioned  script  was  to  input  the  excess  velocity,  the  initial 
spacecraft  mass,  the  final  spacecraft  mass,  thrust,  and  the  Isp;  and  output  AV,  payload 
mass  fraction,  and  time  of  flight. 
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For  this  mission  direct  impulsive/chemical  burns  are  not  a  viable  option  because  of 
the  propellent  mass  requirements  inherent  in  yielding  a  AV  in  the  km/ s  range.  So  by 
necessity,  the  low  thrust/high  specific  impulse  trade  space  was  explored. 

Unlike  the  high  thrust/approximately  impulsive  case  where  explicit  analytical 
equations  are  used,  the  low  thrust  case  directly  utilizes  the  equations  of  motion. 


r  =  u  [41] 

(3.15) 

v2  u 

u- - -  +  A  sin0  [41] 

(3.16) 

r  rl 

uv 

v  = - h  Acos0  [41] 

(3.17) 

r 


Equations  3.15,  3.16,  3.17,  are  the  polar  representation  of  the  two  dimensional 
equations  of  motion.  In  these  equations  r  expresses  the  distance  from  the  gravitational 
body,  u  is  the  radial  velocity,  v  is  the  circumferential  velocity,  A  is  the  acceleration  due  to 
thrust,  and  (p  is  the  polar  thrust  angle.  The  AV  is  calculated  by  integrating  the  equations  of 
motion  over  time  until  the  desired  radial  distance  is  achieved  and  using  the  Ideal  Rocket 
Equation,  Equation  3.3  .  Using  the  equations  of  motion,  the  Ideal  Rocket  Equation, 
MATLAB  ode45  (a  numerical  integrator),  and  by  varying  the  aforementioned  inputs 
Figure  3.8  and  Figure  3.9  were  created.  In  both  Figure  3.8  and  Figure  3.9  the  AV  was 
manipulated  by  changing  the  propellant  mass  and  keeping  the  spacecraft  mass  constant. 
Finally,  the  inputs  were  varied  based  on  attainable  and  demonstrated  thrust  levels  and 
specific  impulses. 

Figure  3.8  displays  the  effect  of  varying  the  thrust  level  on  trip  time  to  550  AU.  For 
this  data  set  the  specific  impulse  was  held  at  a  constant  3000  sec.  From  the  figure  it  is 
clear  that  an  increase  in  thrust  will  result  in  a  requisite  decrease  in  trip  time.  Another  fact 
that  can  be  gleaned  from  this  figure  is  that  for  every  thrust  level,  at  a  given  specific 
impulse,  there  exists  an  optimal  AV  which  minimizes  the  trip  time.  This  fact  is  harder  to 
discern  but  equally  valid.  Therefore,  to  optimize  trip  time  one  must  locate  the  optimal 
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Figure  3.8:  Delta-V  for  Various  Thrust  Levels 


point  for  a  given  specific  impulse  and  thrust  level,  and  use  that  point  to  find  optimal  ratio 
between  the  final  and  initial  masses. 

Figure  3.9  displays  the  effect  of  varying  the  specific  impulse  on  trip  time.  For  this 
data  set  the  thrust  was  held  at  a  constant  100  mN.  It  is  evident  that  an  increase  in  specific 
impulse  will  result  in  a  decrease  in  trip  time.  Also,  similar  to  Figure  3.8,  there  exist  an 
optimal  AV.  Notice  that  the  optimal  AV  is  the  same  for  the  3000  sec  series  in  Figure  3.9 
and  the  100  mN  series  in  Figure  3.8,  as  they  should  be. 

Conclusively,  the  MATLAB  script  along  with  the  data  on  available  thrusters  was  used 
to  calculate  the  mission  AV. 
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Figure  3.9:  Delta-V  for  Various  Specific  Impulses 


3.5.3  Payload  Mass  Fraction.  The  payload  mass  fraction  represents  the 
percentage  of  the  spacecraft  that  remains  for  payload  after  the  power  source  and 
propellant.  The  equation  for  the  payload  mass  fraction  is 

(YsiJspti  0  ) 

mf - j - (m‘  ~  m/) 

fpay  =  - — -  [15]  (3.18) 

nij 

where  as  is  the  power  system  specific  mass  ,  //,  is  the  thruster  efficiency,  in,  is  the  initial 
mass  ,  nif  is  the  final  mass  ,  and  r  is  the  burn  time.  The  specific  mass  has  a  direct  linear 
relationship  with  the  payload  mass  fraction.  On  the  other  hand,  the  thruster  efficiency  has 
more  complicated  relationship  with  the  payload  mass  fraction.  An  increase  in  thruster 
efficiency  causes  an  asymptotic  approach  to  a  maximum  payload  mass  fraction  for  the 
given  variables. 
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3.5.4  Power.  For  a  thruster  to  operate  it  needs  power.  The  jet  power  is  the  amount 
of  power  required  for  the  the  thruster  to  operate  at  a  given  thrust  and  specific  impulse.  The 
equation  for  the  jet  power  is  given  by 


where 


_  mdotv2e 

2 


U  ~  IspgO 


(3.19) 


(3.20) 


and  rridot  is  the  mass  flow  in  kg/s,  T  is  thrust  in  newtons,  and  Isp  is  the  specific  impulse. 

However,  since  thrusters  are  not  100  percent  efficient  a  thruster  efficiency  term,  the 
same  one  from  the  payload  mass  fraction,  is  introduced  into  the  jet  power  equation.  This 
produces  the  equation  for  source  power. 


Ps  = 


mdotv : 


TV, 


n  t 


2rjt  2rjt 


[15] 


(3.21) 


As  a  result,  higher  thrust  and  higher  specific  impulses  require  correspondingly  higher 
powers.  Therefore,  the  thrust  and  specific  impulse  of  the  spacecraft  will  be  limited  by  the 
power  available. 


3.5.5  Gravity  Assist.  Looking  back  at  Figure  3.9  and  Figure  3.8,  to  reach  550  AU 
in  approximately  90  years  would  require  a  AV  in  the  60-70  km/ sec  range.  To  decrease 
that  requirement  the  spacecraft  could  make  use  of  a  gravity  assist.  Gravity  assists  utilize 
the  one  or  multiple  gravity  wells  to  substantially  increase  the  velocity  of  the  spacecraft 
relative  to  the  sun.  Thus,  the  options  available  are:  Jupiter  Gravity  Assist  (JGA),  Saturn 
Gravity  Assist  (SGA),  Uranus  Gravity  Assist  (UGA),  Neptune  Gravity  Assist  (NGA), 
Jupiter  and  Saturn  Gravity  Assist  (JSGA),  Jupiter  and  Uranus  Gravity  Assist  (JUGA), 
Jupiter  and  Neptune  Gravity  Assist  (JNGA), Solar  Gravity  Assist  (SoGA),  Venus  and 
Solar  Gravity  Assist  (VSGA).  Furthermore,  if  the  spacecraft  is  to  be  launched  in  a  specific 
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direction  examining  the  sidereal  period  of  the  planets  produces  the  time  between  launch 
opportunities,  which  is  pictured  in  Figure  3.10. 


180 


Gravity  Assist 


Figure  3.10:  Time  Between  Launch  Opportunities  [8]  [9] 

Figure  3.10  shows  that  the  only  plausible  gravity  assists  are  the  SoGA,  VSGA,  JGA, 
and  SGA.  The  others  are  unreasonable  because  the  time  between  launch  opportunities  is 
too  large,  i.e.  greater  than  50  years.  The  SoGA  and  VSGA  would  require  close 
approaches  to  the  sun  necessitating  the  use  of  large,  heavy  radiation  shielding. 
Furthermore,  the  physics  of  gravity  assist  dictate  that  the  more  massive  the  planet  the 
more  velocity  it  provides;  thus,  since  Saturn  is  much  less  massive  than  Jupiter  out  of  the 
four  remaining  trajectories  the  most  promising  is  the  JGA. 
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The  physics  of  gravity  assists  can  be  understood  using  analytic  expressions.  The 
following  approach  to  gravity  assists  assumes  an  approximately  coplanar  flyby. 
Furthermore,  the  calculations  for  a  gravity  assist  at  Jupiter  require  input  values  for  the 
radius  of  perigee,  rp;  insertion  velocity  at  “infinity”,  v^;  mass  of  the  planet,  My  velocity 
of  the  planet,  Vy  and  the  orientation  angle,  (p.  The  goal  of  the  gravity  assist  calculation  is 
to  find  how  much  AV  the  planet  provides.  The  first  required  expression  is  that  of  the 
semi-major  axis.  The  semi-major  axis  is  defined  as  the  largest  radius  of  the  orbit’s  conic 
section.  This  relation  is  stated  as 

a  =  _^Ml  [IQ]  (3.22) 

vz 

oo 

where  the  only  non-input  variable  is  G  the  gravitational  constant.  Following  the 
calculation  of  the  semi-major  axis  is  the  eccentricity  e,  which  measures  the  eccentricity  of 
the  gravity  assist  orbit. 

2 

e=l  +  7^T  [1Q]  (3-23) 

LrMj 

In  Equation  3.23,  all  of  the  variables  have  been  previously  defined.  Next,  the  value  of 
the  parameter  can  be  calculated  with  the  eccentricity  and  the  semi-major  axis.  The 
parameter  is  a  property  of  conic  sections  dictated  by  the  expression 

p  =  a(  1  -  e2)  [10]  (3.24) 

After  the  parameter  the  next  quantity  of  interest  is  the  periapsis  velocity  vp.  The 
periapsis  velocity  is  the  velocity  at  the  minimum  orbital  radius  from  Jupiter. 

[10]  (3.25) 

This  velocity  can  then  be  used  to  calculate  the  angular  momentum  h  of  the  orbit. 

h  =  rpvp  [10]  (3.26) 

The  values  of  Equations  3.22,  3.23,  3.24,  3.25,  and  3.26  are  all  constants  throughout 
the  calculations.  The  first  equation  of  note  with  a  varying  expression  is  the  true  anomaly. 
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The  true  anomaly  measures  the  position  of  the  spacecraft  in  the  orbit  relative  to  periapsis, 
it  is  negative  prior  to  it  and  positive  following  it.  Thus,  periapsis  has  a  true  anomaly  of 
zero.  True  anomaly  is  given  by  the  expression 


/  =  cos 


1  _  i 
-1  I  r 


[10] 


(3.27) 


where  r  is  the  current  radius  of  the  gravity  assist.  Furthermore,  to  find  the  initial  and  final 
true  anomaly  examine  Equation  3.27  and  allow  r  to  go  to  infinity,  because  the  initial  point 
and  the  final  point  of  the  gravity  assist  are  an  arbitrarily  large  distance  away.  As  a  result 
Equation  3.27  becomes 

.  /  1  \ 

(3.28) 


foe  =  cos  1 1-- 1  [10] 

Following  that,  Equation  3.29  produces  the  current  radius  of  the  gravity  assist  at  any 
point  in  the  orbit. 


[10] 


(3.29) 


1  +  ecosf 

Varying  true  anomaly  from  -  frx  to  +/«,  creates  an  entire  map  of  the  radii  from  start  to 
finish.  These  radii  can  then  be  inserted  into  Equation  3.30  to  yield  its  corresponding 
velocity, 

bn  Mr 

(3.30) 


v=J^i+vi  [101 


Next,  the  range  angle  [3  can  be  thought  of  as  a  measure  of  displacement  from  the 
initial  true  anomaly. 

0  =  /«,+/[  10]  (3.31) 


The  flight  path  angle  y  is  given  by 

y  =  +  cos-1  [10]  (3.32) 

where  the  radius  and  velocity  of  the  orbit  are  constantly  changing.  The  minus  plus  sign  is 
present  because  while  the  true  anomaly  is  negative  the  flight  path  angle  is  also  negative, 
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and  while  the  true  anomaly  is  positive  the  flight  path  angle  is  positive.  Then,  combining 
the  flight  path  angle  and  range  angle  into  an  equation  produces  the  turn  angle  6. 

6  =  (3  -  y  -  90°  [10]  (3.33) 

The  turn  angle  ”is  a  measure  of  how  much  the  spacecraft’s  velocity  orientation  has  been 
rotated  from  its  starting  direction  toward  its  final  direction”  [10].  Finally,  Equation  3.34 
builds  on  the  knowledge  provided  by  the  previous  equations  to  produce  the  spacecraft’s 
velocity  relative  to  the  sun  V. 

V  =  yjv2  +  Vj-  2 vVj  cos(0  +  6)  [10]  (3.34) 

Figure  3.11  demonstrates  the  effect  of  gravity  assists  very  clearly  and  concisely.  It  shows 
the  velocity  profile  of  Voyager  I’s  Jupiter  gravity  assist.  Upon  entering  its  Jupiter  gravity 
assist  maneuver  Voyager  I  possessed  a  velocity  of  12.62  km/ sec  relative  to  the  Sun.  At  its 
departure  from  the  maneuver  Voyager  I’s  velocity  relative  to  the  sun  had  increased  by 
10.77  km/ sec  to  23.39  km/ sec.  Nonetheless,  Voyager  I’s  velocity  relative  to  Jupiter  was 
the  same  on  entry  and  departure. 
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Figure  3.1 1:  Voyager  I  Gravity  Assist  Velocity  Relative  to  the  Sun  [10] 
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3.6  Attitude  Control 


The  attitude  control  subsystem  is  driven  by  the  pointing  requirements.  For  the  FSM 
the  pointing  requirements  are  dictated  not  by  the  mission-pointing  at  the  Sun-but  by 
communications-pointing  at  the  Earth.  The  basic  pointing  requirement  equation  can  be 
expressed  as 

#  =  t  [35]  (3.35) 

h 

where  4>  is  the  pointing  accuracy,  D  is  the  target  diameter,  and  It  is  the  distance  from  the 
spacecraft  to  the  target.  Interestingly,  as  the  pointing  requirement  decreases  it  becomes 
harder  to  maintain  and  more  accurate  sensors  are  necessary.  What  is  more,  while  h  is 
explicitly  stated  D  must  be  calculated.  For  the  case  of  the  FSM,  the  target  diameter  is  the 
distance  subtended  by  half  the  HPBW  at  Earth,  because  with  this  tolerance  the  Earth  will 
remain  within  the  HPBW. 


Figure  3.12:  Pointing  Accuracy  Geometry 
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The  geometry  of  the  pointing  requirement  is  simple.  At  the  apex  of  Figure  3.12  is  the 
spacecraft  transmitter.  The  straight  line  distance  between  the  current  location  of  the 
spacecraft  and  Earth  is  represented  by  h.  The  half  of  the  HPBW,  which  was  mentioned 
earlier  in  the  Communications  and  Data  Handling  subsection  on  page  45,  subtends  the 
target  diameter  at  Earth.  Finally,  6  is  a  quarter  of  the  HPBW  and  was  introduced  to 
simplify  solving  for  the  target  diameter. 

Looking  at  Figure  3.12  the  relationship  between  8  and  D  can  be  stated  as 

D  =  2fctan(0)  (3.36) 


Thus  the  target  diameter  increases  with  the  HPBW.  By  inserting  Equation  3.11  in  for  9 
Equation  3.36  becomes 


D  =  2 h  tan 


HPBW 


2 h  tan 


4  fcHzDa 


(3.37) 


Equation  3.37  illustrates  the  inverse  relationship  between  the  target  diameter  and  the 
antenna  diameter.  Furthermore,  plugging  Equation  3.37  back  into  Equation  3.35  produces 


=  2  tan 


4  fcHzDa 


(3.38) 


which  is  the  pointing  requirement  equation  in  its  simplest  form.  Solving  this  equation 
using  the  specifics  for  the  FSM  produces  a  pointing  accuracy  of  approximately  20  degrees; 
however,  by  limiting  the  pointing  accuracy  to  five  degrees  the  data  rate  will  see  a  requisite 
increase  due  to  smaller  pointing  losses.  Finally,  to  achieve  the  specified  pointing  accuracy 
of  five  degrees  the  spacecraft  needs  a  robust  attitude  control  sensor  and  actuator  suite. 

Table  2.2  displays  the  attitude  sensors  and  actuators  used  on  NASA’s  deep  space 
missions.  The  sensors  include  star  trackers,  sun  sensors,  and  inertial  reference  units.  Sun 
sensors  track  the  spacecraft’s  attitude  by  knowing  its  own  position  on  the  spacecraft 
reference  frame  and  comparing  that  with  the  sun  vector  it  measures.  Unfortunately,  sun 
senors  can  only  track  two  attitude  dimensions  whether  that  be  pitch  and  yaw,  pitch  and 
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roll,  or  yaw  and  roll  [35].  Star  trackers,  on  the  other  hand,  can  target  multiple  stars  and 
produce  information  on  all  three  attitude  dimensions.  Finally,  inertial  reference  units  do 
not  require  an  external  reference.  They  are  capable  of  measuring  torques  on  the  spacecraft 
using  well-calibrated  gyroscopes  and  accelerometers.  The  measured  outputs  from  the 
gyroscopes  and  accelerometers  are  then  translated  into  information  about  the  spacecraft’s 
attitude.  The  FSM  will  use  only  two  of  the  sensors  to  ensure  that  it  can  achieve  the 
pointing  accuracy  required  to  complete  the  mission.  A  sun  sensor  will  not  be  used  because 
operating  at  the  mission  distance  the  sun  will  appear  to  be  just  another  star.  Knowing  the 
spacecraft  attitude  is  not  enough.  The  FSM  must  also  be  able  to  control  its  attitude. 

The  primary  active  deep  space  attitude  actuators  are  thrusters.  Thrusters  are  a  reliable 
solution  for  long  term  attitude  actuation.  They  were  chosen  over  the  inclusion  of  CMGs 
and  reaction  wheels  because  while  CMGs  and  reaction  wheels  would  decrease  the  need 
for  propellant  they  have  a  much  higher  propensity  to  fail  on  a  long  duration  mission 
according  to  Larson  and  Wertz’s  text  [15].  The  principal  thruster  propellant  is 
MMH/N204,  because  of  its  desirable  specific  impulse  and  ability  to  be  stored  for  long 
durations.  Thus  thrusters  will  be  used  as  the  active  actuator  on  the  FSM.  Of  course, 
thrusters  require  propellant  which,  in  turn,  require  propellant  tanks. 

The  amount  of  propellant  required  is  determined  from  the  Ideal  Rocket  Equation, 
Equation  3.3.  In  this  case,  the  known  values  are  the  initial  mass,  specific  impulse,  and  AV 
since  the  amount  of  AF  provided  by  the  thruster  is  predetermined  based  on  an  expectation 
of  disturbance  torques.  And  the  desired  value  is  the  final  mass.  Once  the  final  mass  is 
found  the  propellant  mass  is  the  difference  between  the  initial  and  final  mass. 

The  propellant  tank  size  is  dependant  upon  the  tank  material  and  volume  of  the 
propellant.  For  the  FSM  the  tank  material  will  be  Aluminum-2219,  because  of  its  ultimate 
strength  and  density.  Aluminum-2219  has  an  ultimate  strength  of  .413  GPa  and  density  of 
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2800  kg/ m3  [38].  Moreover,  it  is  compatible  with  both  MMH  and  N204  [38].  Nitrogen 
Tetroxide  and  MMH  have  densities  of  1440  kg/m 3  and  878  kg/m 3  respectively  [38]. 


After  using  Equation  3.3  to  calculate  the  amount  of  total  propellant  required  the  total 
propellant  mass  is  partitioned  based  on  the  desired  fuel  to  oxidizer  mass  ratio.  Then,  the 
following  equation  is  used  to  calculate  the  radius  of  the  spherical  tank. 


r,  = 


3(l.l)Vp 


4  n 


[38] 


(3.39) 


where  Vp  is  the  volume  of  the  propellant.  The  constant  value  of  1.1  takes  care  of  the 
ullage  and  trapped  volumes.  The  ullage  and  trapped  volumes  account  for  the  fact  that  the 
tanks  will  not  be  completely  full  with  propellant  but  will  leave  space  so  that  there  will  be 
enough  pressure  to  push  the  final  propellant  through  the  pipes.  In  this  case  a  conservative 
ten  percent  is  left  empty  in  the  tank.  Next,  Equation  3.40  calculates  the  area  of  the  tank 
using  the  radius. 

A,  =  4 nr;  [38]  (3.40) 


The  tank  wall  thickness  can  be  calculated  using  the  radius  along  with  the  burst  pressure 
and  ultimate  strength. 

=  P8]  (3.41) 

2  Fu 

where 

Pb  =  2  P  me  (3.42) 


In  Equation  3.42,  PMe  is  the  maximum  expected  operating  pressure.  For  both  the  fuel  and 
the  oxidizer  the  maximum  operating  pressure  will  be  around  1.7  MPa  so  that  the  pressure 
in  the  tank  is  sufficient  to  meet  the  demands  of  a  bipropellant  thruster.  The  constant  2  in 
the  equation  is  the  factor  of  safety.  Finally,  the  mass  of  the  tank  can  be  calculated  through 


ms  =  Astspmat  [38]  (3.43) 

where  pmat  is  the  density  of  the  tank  material.  Outside  of  the  active  actuation  of  thrusters 
there  are  also  passive  attitude  actuators. 
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For  deep  space  missions  there  are  only  two  viable  passive  attitude  actuators-nutation 
dampers  and  spin  stabilization.  Nutation  dampers  absorb  torques  on  a  given  axis  and 
convert  them  into  heat.  The  second  form  of  passive  actuation,  spin  stabilization,  uses  ’’the 
conservation  of  angular  momentum  to  maintain  a  constant  inertial  orientation  of  one  of  its 
axes”  [35].  Thus  the  spacecraft  can  remain  pointed  toward  Earth  or  the  Sun  and  be  slewed 
by  the  thrusters  when  necessary.  Spin  stabilization  cannot  be  used  on  the  FSM  during  the 
boost  because  of  the  need  to  keep  the  solar  panels  pointed  at  the  sun. 

3.7  Environmental  Control 

The  environment  of  the  FSM  will  not  be  decidedly  different  then  that  of  previous 
deep  space  missions.  There  are  two  principal  environmental  concerns.  They  are  thermal 
balance  and  irradiation. 

Maintaining  thermal  balance  means  having  the  the  amount  of  heat  coming  into  the 
spacecraft  equal  the  amount  of  heat  being  radiated  away  from  the  spacecraft,  which  is 
mathematically  illustrated  by  the  following  equation 

Qin  =  Qout  (3.44) 

where  heat  is  designated  by  the  quantity  Q.  This  is  the  steady  state  thermal  balance 
equation.  The  quantity  on  the  left  side  of  the  expression  includes  heat  generated  by  the 
spacecraft  and  imparted  on  the  spacecraft  by  planets  and  the  sun.  The  right  side  accounts 
for  heat  radiated  away  by  passive  or  active  means.  Active  means  includes  radiators, 
louvers,  heaters,  and  heat  pipes.  Passive  means  include  blankets  and  reflective  materials. 
The  goal  of  thermal  balance  is  to  keep  the  spacecraft  temperature  in  a  desirable  range. 
Table  3.2  provides  typical  thermal  ranges  for  spacecraft  components.  The  outer 

spacecraft  temperatures  can  range  between  ±200°C  [15].  Accordingly,  a  significant 
amount  of  thermal  protection  is  needed. 
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Table  3.2:  Typical  Thermal  Ranges  for  Spacecraft  Components  [15] 


Component 

Typical  Temperature  Range  (°C) 

Operational 

Survival 

Inertial  Reference  Unit 

0-40 

-10-50 

Star  Trackers 

0-30 

-10-40 

C&DH 

-20  -  60 

-40  -  75 

Hydrazine  Tanks 

15-40 

5-50 

Antennas 

-100-100 

-120-120 

Solar  Panels 

-150-110 

-200-  130 

Total  Range 

15-30 

5-40 

Radioisotope  heaters  units  ( RHU )  will  be  used  as  an  active  means  to  maintain  the 
specific  temperature  range  of  each  component.  These  devices,  discussed  in  the 
Background  and  Theory,  radiate  heat  from  the  decay  of  radioisotopes.  Furthermore,  the 
FSM  will  make  use  of  low-emittance  aluminized  mylar  multi-layered  insulation  blankets 
to  prevent  excessive  heating/heat  loss. 

Not  only  must  thermal  balance  be  maintained  but  heat  must  also  distributed  around 
the  spacecraft.  A  high  concentration  of  heat  in  one  area  or  the  lack  heat  in  another  would 
be  detrimental  to  the  spacecraft.  However,  heat  pipes  will  not  be  used  to  distribute  heat 
throughout  the  spacecraft.  Heat  will  be  distributed  by  thermally  tying  the  spacecraft 
together. 

Like  the  thermal  environment  the  radiation  environment  will  closely  match  that  of 
previous  space  missions  therefore  previous  radiation  hardening  schemes  can  be  used.  The 
novel  threat  the  spacecraft  will  face  is  prolonged  exposure  to  the  same  rate  of  irradiation. 
This  will  increase  the  likelihood  of  a  malfunction.  Nonetheless  providing  protection 
against  irradiation  is  a  simple  matter.  Either  dense  materials  must  be  introduced  to  stop 
the  high  energy  particles  from  reaching  the  electronics  or  high  reflectance  materials  to 
reflect  particles  away  from  the  spacecraft.  There  are  also  internal  radiation  threats  to  the 
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spacecraft  due  to  the  use  of  RTGs.  As  a  result,  the  RTGs  must  be  placed  to  minimize  their 
effect  on  the  primary  payload. 

3.8  Conclusion 

In  conclusion  this  chapter  analyzed  and  discussed  how  to  optimize  the  gravitational 
lens  mission.  The  following  chapter  will  utilize  the  findings  of  this  chapter  to  design  the 
FSM. 
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4  Design 


4.1  Design  Flow 

The  spacecraft  requirements  have  been  presented;  thus,  the  purpose  of  this  chapter  is 
to  synthesize  the  outcomes  of  the  previous  chapters  and  the  spacecraft  requirements  into 
one  coherent  spacecraft  design  through  the  iterative  process  shown  in  the  Design  Flow 
Diagram,  Figure  4.1. 

At  the  top  of  the  diagram  is  the  propulsion  system.  As  stated  in  Section  3.2,  the  type 
of  propulsion  system,  high  specific  impulse  or  low  specific  impulse,  decided  which  launch 
vehicle  is  utilized.  What  is  more,  the  propulsion  system  also  has  a  direct  effect  on  the  AV. 
It  is  for  these  two  reasons  that  the  propulsion  system  was  chosen  as  the  starting  point  of 
the  iterative  process. 

Under  the  propulsion  system  is  the  launch  vehicle.  For  this  analysis,  the  launch 
vehicle  was  directly  dependent  on  only  the  propulsion  subsytem.  The  launch  vehicle  then 
puts  limits  on  the  mass  and  volume  of  the  entire  spacecraft.  Progressing  to  the  next  level, 
every  spacecraft  subsystem  was  interdependent. 

In  conclusion,  there  were  several  iterations  through  the  design  flow  due  to  the 
hierarchy  and  coupling  of  dependencies.  The  following  sections  present  the  final  result  of 
these  iterations. 
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Figure  4.1:  Design  Flow  Diagram 
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4.2  Launch  Vehicle 


The  analysis  in  Chapter  III  stipulated  that  when  a  low  thrust  propulsion  system  is 
present  the  optimal  launch  vehicle  choice  is  the  one  that  provides  the  largest  delivered 
mass  to  the  parabolic  excess  energy.  Therefore,  from  the  chart  of  options  on  pg.  35,  the 
obvious  choice  is  the  Delta  IV-H/Star  48/Star  37  configuration.  This  provides  a  delivered 
mass  of  8608  kg;  thus,  limiting  the  total  spacecraft  mass  to  8608  kg. 


4.292  m 


Payload 

Module 


12.192  m 


5.080  m 


Figure  4.2:  Payload  Fairing  [11] 


With  that,  the  launch  vehicle’s  payload  fairing  constrains  the  dimensions  of  the 
spacecraft.  Figure  4.2  displays  the  standard  Delta  IV  metallic  payload  fairing  used  by  the 
United  States  government.  Clearly  from  this  illustration  the  majority  of  spacecraft  must  be 
less  than  5.08  meters  in  diameter  and  12.19  meters  in  height  in  its  stowed  configuration. 
Or  if  it  is  rectangular  it  must  have  sides  that  measure  less  than  3.59  meters  so  that  it  will  fit 
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inside  the  circular  envelope.  Finally,  The  spacecraft  may  have  appendages  that  extend  up 
to  4.29  meters  from  the  apex  of  the  main  spacecraft. 

4.3  Spacecraft 

4.3.1  Payload.  The  primary  payload  is  a  12-meter  diameter  parabolic  center-feed 
antenna/radio  telescope.  It  is  this  size  so  that  it  can  accurately  measure  variations  in  the 
frequency  magnification,  detect  possible  communications  sent  through  the  gravitational 
lens,  and  take  images  images  using  radio  telescopy. 

The  dual  purpose  primary  payload  will  be  a  deployable  metallic  mesh  similar  to 
Galileo  and  the  Tracking  and  Data  Relay  Satellites  [42].  The  metallic  mesh  has  an 
approximate  specific  mass  of  .12  kg/m2.  With  an  area  of  113.10  m2  the  antenna  mass  will 
be  13.5  kg.  The  antenna  mass  must  also  account  for  the  mass  of  the  feed  array  structure, 
which  is  approximately  10  kg  [15].  Therefore,  the  total  mass  of  the  primary  payload  is 
23.5  kg.  The  power  seen  by  the  primary  payload  will  vary  with  the  degradation  of  the 
power  source,  but  the  goal  is  to  maintain  15  W.  The  secondary  payloads  were  outlined  in 
the  Background  and  Theory.  Figure  4.1  consolidates  the  mass  and  power  requirements  of 
the  payloads. 

4.3.2  Communications  and  Data  Handling.  The  primary  payload  will  also  serve 
as  the  high  gain  antenna.  Accompanying  it  will  be  a  medium  gain  parabolic  antenna 
whose  mission  it  will  be  to  acquire  Earth  and  to  maintain  beacon  monitoring  between  the 
spacecraft  and  the  DSN,  located  in  California,  USA;  Madrid,  Spain;  and  Canberra, 
Australia. 

Using  Equation  3.6  the  gain  of  the  primary  HGA  is  46.59  dB.  Since  the  goal  of  the 
MGA  is  to  acquire  Earth  it  is  desirable  for  it  to  have  a  high  HPBW.  Furthermore, 
according  to  Equation  3.11  to  do  that  the  diameter  of  the  antenna  needs  to  be  minimized. 
In  contrast,  the  gain  equation  suggests  a  large  diameter.  As  a  result,  the  chosen  MGA 
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Table  4.1:  Payload  Performance  [15] 


Payload 

Mission 

Mass 

(kg) 

Data 

Rate 

(Mbps) 

Power 

(W) 

Attitude 

Requirements 

Section 

VBSDC 

Dust  Detection 

1.6 

0.02 

5.1 

- 

2.3 

FGM 

Magnetic  Field 
Measurements 

15.25 

0.03 

11.3 

2.3 

SHM 

Magnetic  Field 
Measurements 

9.08 

0.01 

6.5 

" 

2.3 

PEPSSI 

Detect  Neutral 

and  Ionized 

Particles 

1.5 

0.01 

2.5 

2.3 

Radio 

Telescope 

Characterize 

Gravitational 

Lens 

23.5 

2 

15 

5° 

4.3.1 

TOTAL 

50.93 

2.07 

40.4 

diameter  is  2-meters  in  diameter.  This  choice  balances  the  need  for  a  respectable  gain  in 
deep  space  and  the  desire  for  a  higher  HPBW.  The  MGA,  according  to  its  area  and  the 
specific  mass  of  metallic  mesh,  weighs  .38  kg.  Concurrently,  the  array  structure  for  the 
MGA  weighs  5  kg.  Its  area  also  corresponds  to  a  HPBW  of  4.58°,  which  does  not  sound 
large  but,  in  fact,  it  is  six  times  larger  than  the  HGA  HPBW  of  0.7636°.  The  gain  of  the 
MGA  is  31.027  dB.  The  corresponding  data  rates  will  be  addressed  in  Section  4.5. 

Finally,  in  regard  to  data  handling,  the  processor  will  be  the  RAD750  and  the  FSM 
will  use  a  64  Gbit  solid  state  recorder  as  stated  in  Chapter  III.  The  data  handling  and 
processing  functions  of  the  spacecraft  are  expected  to  draw  approximately  15  W  [15]. 
Figure  4.3  shows  the  payload/communications  architecture  for  the  FSM. 


4.3.3  Power.  The  ASRG  has  only  attained  TRL  6  thus  it  cannot  be  considered 
COTS.  As  a  result,  the  flight  proven  GPHS  was  chosen  based  on  its  superior  power 
density  compared  to  the  other  options  in  Table  2.4.  The  GPHS  is  pictured  in  Figure  4.4. 
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Figure  4.3:  Payload  and  Communications  Architecture.  Note:  The  primary  payload  serves 
as  a  radio  telescope  and  HGA. 


GPHS-RTG 


Cooling  Tube 
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RIG  Mounting  Silicon-Germanium  Source  Support 

Range  ^uhi-foil  (Si-C.e)  Unicouple 


Figure  4.4:  General  Purpose  Heat  Source  RTG  [12] 
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The  number  of  generators  required  is  usually  determined  by  the  EOL  requirements, 
which  were  taken  into  account,  but  for  the  FSM  the  power  requirement  dictated  by  the 
thruster  superseded  that  of  the  EOL.  Thus  the  beginning  of  life  generator  power  level  is 
4914  W  so  that  thrust  can  be  maintained  after  jettisoning  the  boost  phase  at  Jupiter.  This 
number  was  found  by  trading  the  mass  of  the  power  system,  propellant  mass,  trip  time,  the 
power  provided  by  the  power  system,  and  the  propulsion  levels.  Figure  4.5  displays  the 
result.  Therefore,  the  optimum  number  of  GPHS  RTGs  for  the  FSM  is  20,  because  it 
corresponds  to  the  minimum  trip  time. 
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Figure  4.5:  Power  Source  Trade 
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This  corresponds  to  an  EOL  power  of  approximately  1 184  W  at  1000  AU,  which  is 
more  than  sufficient  for  mission  operations  based  on  the  total  power  draw  of  the  spacecraft 
presented  in  Table  4.1 1. 

At  57.8  kg  per  GPHS  the  total  weight  of  the  GPHS  power  system  is  1156  kg,  which 
requires  156  kg  of  plutonium-238  .  This  presents  an  intractable  problem  that  will  be 
discussed  in  Chapter  V.  The  stellar  performance  demonstrated  by  its  forbearers  on 
Voyager  and  Pioneer  imply  that  these  generators  can  operate  for  extremely  long  durations. 
Finally,  Figure  4.14  displays  the  expected  degradation  of  the  power  source  as  well  as  the 
power  draw. 

The  second  challenge  for  the  power  system  of  the  FSM  is  to  provide  sufficient  power 
for  the  boost  phase  prior  to  the  JGA.  Analysis  of  this  problem  showed  that  the  optimal 
power  strategy  would  be  to  utilize  a  solar  electric  phase.  Furthermore,  the  solar  array 
should  be  able  to  provide  sufficient  power  for  the  thruster  throttling  profile.  A  good 
candidate  is  the  Spectrolab  Ultra  Triple  Junction  GaInP2/GaAs/Ge  solar  cells.  A  60 -nr, 
900  kg  optimized  array  of  these  cells  would  theoretically  produce  approximately  22  kW 
of  power  at  air  mass  zero  (AMO)  [38]  [31].  That  may  seem  heavy  but  it  would  take  more 
than  5000  kg  of  this  mission’s  RTGs  to  produce  a  similar  power. 

The  cells  have  a  28%  BOF  efficiency,  and  a  24.3%  EOF  efficiency,  where  EOF  is 
defined  as  after  1  MeV  electron  fluence  of  1  x  1015^  [31].  Also,  since  solar  arrays  are 
sensitive  to  the  sun  angle  the  attitude  control  system  will  keep  the  FSM’s  solar  array 
pointed  toward  the  sun.  Finally,  as  the  efficiency  degrades  and  the  solar  panels  travel  away 
from  their  optimal  range  the  power  will  drop  off  steeply;  hence,  when  the  spacecraft 
reaches  Jupiter  the  solar  stage  will  be  jettisoned.  Table  4.2  consolidates  the  information 
presented  in  this  section. 
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Table  4.2:  Power  Performance 


BOL  Power  (kW) 

Total  Dimensions  (m) 

Total  Mass  (kg) 

ASRG 

4.91 

.422  x  1.14 

1156 

Solar  Panels 

22 

3x  20 

900 

Total 

26.91 

- 

2056 

4.3.4  Propulsion.  The  NASA  Evolutionary  Xenon  Thruster  NEXT  is  by  far  the 
best  thruster  choice.  It  has  records  for  the  most  hours  of  operation,  highest  propellant 
throughput,  greatest  total  impulse  demonstrated,  and  longest  hollow  cathode  operation 
[43].  More  specifically,  it  has  demonstrated  the  ability  to  process  greater  632  kg  of  Xenon 
over  4.2  years. 

The  NEXT  thruster  can  be  throttled  through  its  forty  throttle  levels  from  6.9  kW 
down  to  .5  kW  [43].  This  corresponds  to  a  specific  impulse  range  from  approximately 
4300  sec  to  1480  sec  and  thrust  range  of  25.4  mN  to  234.6  mN  [43][44],  Table  4.3  shows 
the  sizing  and  power  requirements  of  the  thruster  while  Figure  4.6  displays  the  entire 
thruster  system. 


Table  4.3:  NEXT  Sizing  and  Power  [13] 


Mass  (kg) 

Dimensions  (cm) 

Power  (W) 

Thruster 

12.7 

Diameter  -  55 
Length  -  44 

6900-500 

Power  Processing 
Unit 

33.9 

42  x  53  x  14 

N/A 

Propellant 

Management 

System 

5 

High  Pressure 
Assembly 

33x15x6.4 

Low  Pressure 
Assembly 
38x30x6.4 

High  Pressure 
Assembly 

4.3 

Low  Pressure 
Assembly 

15.9 

Gimbal 

6 

61x72 

N/A 

Total 

57.6 

Thus  far  NEXT  has  processed  632  kg  of  Xenon  in  the  course  of  the  long  duration 
testing.  Yet  it  has  demonstrated  only  slightly  decreasing  efficiency.  Whereas  NSTAR, 
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Figure  4.6:  NEXT  Thruster  System  [13] 

NEXT’S  predecessor,  lost  9%  efficiency  over  the  same  period  NEXT  has  lost  only  2.5% 
engine  efficiency  [43].  Furthermore,  NEXT  is  forecasted  to  be  capable  of  processing  800 
kg  at  the  maximum  power  level,  failing  primarily  due  to  accelerator  grid  erosion,  which  is 
the  result  of  charge-exchange  erosion  [43].  Charge-exchange  erosion  is  based  on 
propellant  throughput  and  input  powers  meaning  the  thruster  is  life  limited  by  the  amount 
of  propellant  processed. 

The  analysis  which  was  used  to  chose  the  number  of  power  sources  also  determined 
the  propellant  load  of  4885  kg  of  xenon.  This  is  what  remains  of  the  delivered  mass  after 
the  spacecraft  dry  mass-2514  kg-and  the  boost  phase  dry  mass-1200  kg-are  taken  into 
account.  Then,  the  required  number  of  thrusters  was  derived  from  the  propellant  mass  by 
taking  the  total  propellant  mass  and  dividing  it  by  640  kg,  which  is  the  predicted 
maximum  throughput  minus  twenty  percent  for  margin.  The  resulting  number  of  thrusters 
is  7.63  thrusters.  Hence,  a  conservative  eight  thrusters  are  required  to  process  that  amount 
of  propellant.  Meaning  that  each  thruster  operating  one  at  a  time  will  process,  on  average, 
610  kg  of  Xenon.  Eight  thrusters  require  the  components  listed  in  Table  4.4.  Finally, 
Equation  4.1  displays  the  mass  of  the  xenon  tank  as  a  function  of  the  mass  of  xenon 
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propellant. 


mXe  tank  =  52  +  0.075 ntprop  +  .154 mlrop  [38]  (4.1) 

The  arrangement  of  the  spacecraft  and  the  reasons  for  it  will  be  discussed  in  more  detail  in 
Section  4.5.3  and  Section  4.4,  but  at  this  moment  it  is  important  to  note  that  the  spacecraft 
will  be  split  into  two  segments.  The  first  segment  will  be  a  solar  electric  boost  stage  and 
the  second  segment  will  be  the  main  spacecraft.  It  is  important  to  note  at  this  point 
because  the  xenon  will  be  split  into  two  tanks-one  in  the  boost  phase  and  one  in  the  solar 
electric  phase-to  produce  a  higher  total  AV.  Therefore,  using  Equation  4.1,  the  mass  of 
the  first  tank  is  308.7  kg  and  the  second  is  217.13  kg,  which  correspond  to  propellant 
throughputs  of  2996  kg  and  1888  kg  respectively. 


Table  4.4:  NEXT  Component  List  [13] 


Number  of  Components  Mass  Power  (W) 

HPA 

1 

1.9  1.9 

LPA 

4 

12.4  3.1 

PPU 

4 

135.6  610-7220 

Gimbal 

4 

24 

Thruster 

4 

50.8  500-6900 

Total 

17 

224.7  615-7225 

4.3.5  Attitude. 

4.3.5. 1  Attitude  Determination.  Attitude  determination  will  be  managed  by 
star  trackers  and  an  inertial  reference  unit.  This  section  will  present  a  COTS  component 
that  meets  the  mission  requirement. 

The  Ball  Aerospace  CT-633  will  serve  as  the  COTS  star  tracker.  It  is  capable  of 
tracking  up  to  five  stars  and  providing  an  angular  accuracy  of  25  arc  seconds,  temporal 
accuracy  of  6.5  arc  seconds,  and  attitude  accuracy  of  15  arc  seconds  at  the  end  of  life  [45]. 
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The  CT-633  has  a  power  draw  of  9  W,  mass  of  2.49  kg  and  dimensions  measuring  .135  m 
in  diameter  by  .142  m  in  length  [45]. 

The  inertial  reference  unit  chosen  for  the  FSM  is  the  Honeywell  Miniature  Inertial 
Measurement  Unit  ( MIMU ).  It  can  sense  attitude  changes  in  all  three  axes  in  the  range  of 
±375  deg  I  sec  [46].  It  has  a  mass  of  4.44  kg,  power  consumption  of  22  W,  and  measures 
.233  in  in  diameter  by  .169  m  in  height  [46]. 


Table  4.5:  Attitude  Determination  Mass  and  Power 


Mass  (kg) 

Power  (W) 

Star 

Tracker 

2.49 

9 

MIMU 

4.4 

22 

Total 

6.89 

31 

43.5.2  Attitude  Control.  The  active  attitude  actuation  of  the  spacecraft  will 
be  handled  by  twelve  .8  N  N2O4/MMH  thrusters.  At  a  1-to-l  mixture  ratio  the  specific 
impulse  of  the  thruster  is  318  sec  [38].  Furthermore,  similar  thrusters  were  used  on  both 
of  the  Voyager  missions  and  New  Horizons.  The  New  Horizons  qualification  testing  found 
that  thrusters  are  reliable  up  to  at  least  400,000  cycles  [5].  On  the  Voyager  mission  these 
thrusters  lasted  for  more  than  500,000  cycles  [5].  Moreover,  the  FSM  attitude  control 
profile  requires  5  corrections  a  month.  Over  two  hundred  years  that  corresponds  to 
240,000  cycles. 

Similar  thrusters  weigh  in  a  range  of  .2  kg  to  .4  kg  each  [47].  Thus,  for  the  FSM  each 
thruster  will  be  assumed  to  have  a  mass  of  .4  kg  and  dimensions  of  .4  m  in  diameter  and 
.15  m  in  height.  Also,  the  catalyst  bed  for  each  thruster  will  require  2.2  W  of  power  and 
the  entire  thruster  system  can  draw  up  to  50  W  [15]. 


76 


Thruster  performance  was  examined  by  looking  at  the  slew  rate  of  the  spacecraft  in 
two  different  configurations.  The  first  configuration  was  with  the  boost  phase  attached  and 
the  second  configuration  was  with  the  boost  phase  removed.  The  first  step  in  this  analysis 
was  to  calculate  the  mass  moments  of  inertia  for  the  spacecraft.  This  was  done  by  using 
the  principles  of  superposition  and  the  parallel  axis  theorem.  Some  basic  assumptions 
were  necessary. 

Two  assumptions  were  made  for  the  mass  moment  of  inertia  analysis.  The  first 
assumption  of  the  analysis  was  that  the  spacecraft  and  its  appendages  could  accurately  be 
modeled  as  simple  three-dimensional  shapes,  e.g.  rectangular  prisms  and  cylinders.  The 
second  assumption  of  the  analysis  was  that  the  mass  was  distributed  evenly  throughout 
these  sections  so  that  the  center  of  mass  of  each  of  these  shapes  was  also  the  centroid. 

Table  shows  the  calculated  moments  of  inertia  for  the  spacecraft.  The  axes  of  the 
spacecraft  are  arranged  such  that  the  x-axis  is  parallel  to  the  face  of  the  solar  panels,  the 
z-axis  goes  traverses  the  spacecraft  from  bottom  to  top  and  the  y-axis  completes  the 
right-hand  rule.  The  moments  of  inertia  for  the  spacecraft  with  the  boost  phase  are  much 
larger  than  those  for  the  spacecraft  without  the  boost  phase,  which  is  to  be  expected  due  to 
the  size  of  the  solar  panels,  the  body  of  the  boost  phase,  and  amount  of  propellant  used 
between  the  two  analysis  points. 


Table  4.6:  FSM  Mass  Moments  of  Inertia 


Axis 

With  Boost  Phase 

Without  Boost 

(kg-m2) 

Phase  (kg-m2) 

46190  1782 
14011  5555 
16174  3414 


Then  the  slewing  rates,  Equation  4.2,  were  calculated  using  the  moments  of  inertia  /, 
the  force  of  the  thrusters  F,  the  thruster  burn  time  t,  and  the  thruster’s  distance  from  the 


77 


center  of  mass  L. 


©  =  —  [15]  (4.2) 

The  results  are  shown  in  Table  4.7.  For  the  spacecraft  with  its  boost  phase  the  thrust 
time  was  assumed  to  be  5  seconds  and  for  the  spacecraft  sans  the  boost  phase  the  thrust 
time  was  assumed  to  be  1  second.  This  difference  is  to  aid  the  thrusters  with  the  large 
variations  in  the  size  of  the  mass  moments  of  inertia.  The  distance  from  the  thruster  to  the 
center  of  mass  for  both  cases  is  approximately  1.5  meters  since  the  diameter  and  height  of 
the  main  spacecraft  body  is  3  meters. 


Table  4.7:  FSM  Slew  Rate 


Axis  With  Boost  Phase  Without  Boost 

(deg/sec)  Phase  (deg/sec) 


.0055 

.009 


Based  on  the  number  of  thruster  cycles  the  FSM  disturbance  environment  is 
comparable  to  previous  missions,  e.g.  New  Horizons,  so  the  thrusters  will  provide  a 
similar  130  m/s  of  AF  for  margin  (unexpected  and  disturbance)  torques  and  200  m/s  of 
AF  for  alignment  and  slewing  to  meet  the  5  degree  pointing  requirement  of  the 
communications  system.  This  corresponds  to  a  total  propellant  mass  of  150.72  kg  with 
75.36  kg  of  N204  and  75.36  kg  of  MMH.  Table  4.8  was  created  using  the  propellant 
masses  and  equations  from  Section  3.6. 

Attitude  thrusters  also  require  a  propellant  pressure  system  to  pump  the  propellant 
out  of  the  chamber  and  into  the  thrusters.  The  propellant  pressure  system  on  the  FSM  will 
have  a  mass  of  10  kg.  In  conclusion,  Figure  4.7  diagrams  the  attitude  control  subsystem 
on  the  FSM. 
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Table  4.8:  Propellant  Tank  Design 


N2O4  MMH 


Propellant  Mass  (kg) 

75.36 

75.36 

Propellant  Volume  (  m3  ) 

.0523 

.0858 

Burst  Pressure  (Pa) 

1.72  x  106 

1.72  x  106 

Tank  Radius  (  m  ) 

.2395 

.2824 

Tank  Area  ( m2 ) 

.7208 

1.003 

Tank  Thickness  (m) 

4.99  x  10-4 

5.88x10-4 

Tank  Mass  (kg) 

1.01 

1.65 

4.3.6  Environmental  Control.  Environmental  control  is  one  of  the  simplest 
subsystems.  Furthermore,  since  the  environment  of  the  FSM  is  not  decidedly  different 
than  that  of  previous  space  missions  the  majority  of  this  subsystem  will  be  carried  over 
from  those  aforementioned  missions.  Nonetheless,  the  environmental  control  subsystem 
still  serves  as  a  multifaceted  defense  against  the  space  environment . 

Its  first  line  of  defense  is  on  the  exterior  of  the  spacecraft  in  the  form  of  MFI 
blankets.  The  number  of  layers  of  insulation  was  determined  to  be  forty.  This  number  was 
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derived  from  the  steady  state  heat  transfer  analysis  based  on  the  following  equation  which 
was  itself  derived  from  the  thermal  balance  equation, 


6  = 


Qb 


<rMTl/c  -  Ti) 


[15] 


(4.3) 


In  Equation  4.3  ,  e  is  the  emissivity  of  the  spacecraft,  Qin  is  the  internally  generated  heat 
plus  the  solar  flux,  cr  is  the  Stefan-Boltzmann  constant,  A  is  the  exterior  area  of  the 
spacecraft,  Ts/c  is  the  temperature  of  the  spacecraft,  and  7'OT  is  the  temperature  of  the 
surroundings.  Furthermore,  the  spacecraft  temperature  was  assumed  to  be  295  K  since 
that  is  the  mean  acceptable  operating  temperature  for  electronics  as  defined  by  Table  3.2 
and  the  temperature  of  the  surroundings  was  assumed  to  be  2.7  K. 

The  analysis  examined  four  different  points  in  the  spacecraft’s  trajectory:  spacecraft 
at  Earth  with  thruster  on,  spacecraft  at  Jupiter  with  thrusters  off,  spacecraft  post-Jupiter 
with  thrusters  on,  and  spacecraft  after  thruster  burnout.  The  internal  heat  generation  was 
varied  based  on  these  four  points  show  in  Figure  4.9  .  The  spacecraft  generates  the  least 
heat  in  the  fourth  case  when  the  spacecraft  thrusters  have  finished  the  thrust  profile  and  the 
solar  flux  has  declined  dramatically.  As  a  result,  the  emissivity  required  to  maintain 
thermal  balance  came  out  to  be  the  lowest  at  this  point.  More  specifically,  the  emissivity 
required  was  .0020.  That  correlates  to  approximately  forty  layers  of  MLI  [15].  At  that 
thickness  the  multi-layered  insulation  blankets  to  be  used  on  this  mission  have  a  specific 
mass  of  .73  kg/m 2  [15].  Accordingly,  the  total  mass  of  the  MLI  blankets  is  41.61  kg  since 
the  spacecraft  main  body  is  57  m2. 


Since  the  number  of  layers  of  MLI  is  defined  based  on  the  coldest  case  radiators  are 
needed  to  raise  the  emissivity  when  the  spacecraft  heat  generation  is  higher.  When  the 
spacecraft  is  close  to  the  Sun  and  firing  thrusters  it  will  be  the  warmest,  and  at  that 
temperature  the  emissivity  required  to  maintain  thermal  balance  is  .26.  This  is  within  the 
capability  of  COTS  radiators.  Thus  the  flight  qualified  radiator  louver  assembly  will  have 
a  maximum  emittance  of  .88  to  bring  the  total  emittance  of  the  spacecraft  up  to  .26  from 
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Table  4.9:  Spacecraft  Heat  Balance 


Condition 

Qin  (kW) 

Emissivity 

Thrusting  at  Earth 

6.5 

.2655 

Jupiter  w/  Engine  Off 

.839 

.0343 

Post  Jupiter  Thrusting 

1.928 

.0787 

Post  Burnout 

.051 

.0020 

the  radiator  off/louver  closed  emittance  of  .20.  The  radiator  louver  assembly  has 
dimensions  of  45.7  cm  by  58.2  cm,  an  approximate  specific  mass  of  4.5  kg/m 2  ,  and  will 
draw  up  to  20  W[15]. 

While  MLI  blankets  and  radiator  louver  assemblies  will  maintain  the  temperature  in 
the  main  spacecraft  the  FSM  will  use  seven  of  the  RHUs  described  in  the  Background  & 
Theory  to  maintain  the  thermal  balance  of  external  sensors. 

The  final  and  most  important  line  of  defense  against  the  space  environment  is 
radiation  hardening.  As  a  result,  the  spacecraft  bus  will  be  encapsulated  in  a  tantalum 
radiation  vault.  Tantalum  has  a  density  of  16.69  -fjp,  which  means  it  is  almost  fifty  percent 
more  dense  than  lead.  This  will  keep  the  radiation  level  the  bus  sees  under  the  100  krad 
mission  requirement.  The  vault  will  be  0.35  cm  in  thickness  and  have  dimensions  of 
.75  x  .75  x  .75  meters.  Thus,  the  vault  will  be  189  kg. 

4.4  Model 

This  section  displays  models  of  the  spacecraft  design. 
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3  m 


NEXT 

-  Two  are  behind  the  two  shown 


High  Pressure  Assembly 

-  0.33  x  0.15  x  .064  m 


Electronics  Vault 
-ADCS 

-  Processor 

-  Power  Regulation 

-  Communications 

-  Environmental  Control 

-  Thruster  Control 

-  Payload  Processing 

-  SSR 


Low  Pressure  Assembly 

-  0.38  x  .305  x.064  m 

-  Two  are  behind  the  two  shown 


Power  Processing  Unit 

-  0.42  x  0.53  x  .14  m 

-  Two  are  behind  the  two  shown 


Figure  4.8:  Internal  Configuration 


•  Helium  pressurant  present  in  Xenon,  N204  and  MMH  tanks 


•  All  items  fit  within  the  3  x  3  meter  spacecraft  body 


•  Width  was  necessitated  by  propellant  tanks 
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Figure  4.9:  Stowed  Configuration 


•  Shows  the  stowed  configuration 

•  Fits  within  the  payload  and  base  module  of  the  payload  fairing 

Payload  fairing  shown  in  pink 
10.56  meters  tall 
4.76  meters  at  widest 

•  HGA/Radio  Antenna  wrapped  around  antenna  stem 

•  Solar  panel  folded  in  thirds  and  retracted 
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Figure  4.10:  Zoomed  Out  View 

•  HGA/Radio  Antenna  in  deployed  configuration 

•  HGA/Radio  Antenna  is  12  meters  in  diameter 

•  Antenna  Stem  is  6  meters  in  height 

•  MGA  is  2  meters  in  diameter 

•  Four  next  thrusters  in  a  cluster 

•  Magnetometers  on  booms 
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Figure  4.11:  RTG  View 


.  GPHS 

Four  per  boom 

Five  booms  on  the  spacecraft,  for  a  total  of  twenty 
•  MLI 

Represented  by  the  bronze  exterior  coloring 
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StarTrac 


-i 


VBSDC 


Attitude 

Thruster 


Radiator/ 

Louver 


Magnetometers 


Figure  4.12:  Zoomed  In  View 


•  Shows  all  secondary  payloads 

Venetia  Burney  Student  Dust  Counter 
Flux  Gate  Magnetometer 
Scalar  Helium  Magnetometer 

Pluto  Energetic  Particle  Spectrometer  Science  Investigation 
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Figure  4.13:  Spacecraft  with  Boost  Phase 

•  Boost  Phase  Configuration 

•  Solar  Panels 

3  x  10  m  each 

•  Thrusters 

These  four  thrusters  bring  the  total  to  eight 

4.5  Performance  Analysis 

4.5.1  Mass.  Table  4.10  accounts  for  the  masses  present  on  the  spacecraft.  If  an 
exact  number  was  not  known  a  25%  margin  was  added  to  a  best  guess.  Thus,  the  total 
mass  of  the  spacecraft  minus  the  weight  of  NEXT  propellant  is  2345.93  kg.  The  total  mass 
of  the  boost  phase  minus  the  NEXT  propellant  is  approximately  1350  kg  and  the  mass  of 
the  NEXT  xenon  propellant  is  4885  kg  bringing  the  total  launch  mass  to  8580.93  kg. 
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Table  4.10:  Spacecraft  Dry  Mass  Budget 


Mass  (kg)  Margin 

Section 

Attitude  Control  Propellant 
Tanks 

2.66 

.665 

4.3. 5. 2 

Attitude  Control  Propellant 

150.72 

37.68 

4.3. 5. 2 

Attitude  Determination 

6.89 

- 

Xenon  Propellant  Tank 

217.13 

54.28 

4.3.4 

Cabling 

20 

5 

- 

MLI  Blankets 

41.61 

10.40 

4.3.6 

RHU  (10) 

.397 

- 

2.4.4 

Primary  Payload  and  MGA 

23.73 

5.93 

4.3.1 

GPHS  (20) 

1156 

- 

4.3.3 

Data  Processing/Handling 

15 

3.75 

- 

Thrusters  (4) 

224.7 

- 

4.3.4 

Louvers/Radiators  (10) 

11.96 

- 

4.3.6 

Secondary  Payloads 

27.43 

- 

4.3.1 

Radiation  Vault 

189 

47.25 

4.3.6 

Structure 

75 

18.75 

- 

Total  (with  Margin) 

2162.23(2345.93) 

4.5.2  Power.  Table  4. 1 1  presents  the  power  draw  of  all  of  the  spacecraft 
subsystems  along  with  a  25%  margin.  The  total  power  draw  is  well  under  the  power 
provided  by  the  RTGs  displayed  in  Figure  4.14.  In  the  figure  the  spike  represents  the 
post-Jupiter  thrust  profile  which  will  be  discussed  in  the  following  section.  Furthermore,  it 
is  also  clear  from  the  figure  that  the  spacecraft  should  provide  enough  power  to  send  data 
for  over  200  years. 


Table  4.11:  Spacecraft  Power  Consumption 


Power  fW) 

Margin  fW) 

Section 

Secondary  Payloads 

25.4 

- 

4.3.1 

Communications 

15 

3.75 

4.3.2 

Attitude  Thrusters 

50 

12.5 

4. 3. 5. 2 

Attitude  Sensors 

31 

7.75 

4.3. 5.1 

Thermal 

20 

5 

4.3.6 

Data 

Processing/Handling 

15 

3.75 

4.3.2 

Total  (w/  Margin) 

156  (188.75) 
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5000 


....  Power  Draw 
—  GPHS  Power 


Time  (Years) 


Figure  4.14:  RTG  Power  Level  and  Spacecraft  Power  Draw 


4.5.3  Trajectory  Analysis.  The  trajectory  analysis  was  split  into  three  subsections 
to  facilitate  the  analysis.  The  first  subsection  is  the  Pre-Jupiter  Phase.  This  phase 
commences  at  the  end  of  the  launch  vehicle  insertion.  Thus,  at  the  start  of  this  phase  the 
spacecraft  has  the  same  velocity  relative  to  the  Sun  as  Earth-29.78  km/ sec.  Furthermore, 
this  phase  lasts  until  the  Jupiter  gravity  assist  at  which  point  the  Jupiter  Phase  begins  using 
the  final  values  from  the  first  phase  as  its  initial  conditions.  The  Jupiter  Phase  consists 
solely  of  the  gravity  assist.  Finally,  the  last  phase  is  the  Post-Jupiter  Phase.  This  phase 
lasts  from  the  termination  of  the  gravity  assist  until  the  spacecraft  reaches  550  AU  where 
it  will  begin  to  examine  the  gravitational  lens. 


89 


4.5.3. 1  Pre-Jupiter  Phase.  The  spacecraft  trajectory  starts  at  Earth.  Further,  it 
is  assumed  that  Earth  and  Jupiter  are  properly  aligned  for  rendezvous.  The  initial  angular 
position  of  Earth  relative  to  the  Sun  is  taken  as  the  initial  angle  of  the  trajectory.  Meaning 
the  trajectory  will  begin  at  zero  degrees.  The  complete  trajectory  from  the  Earth  to  Jupiter 
is  shown  in  Figure  4.15.  Since  the  trajectory  is  using  low  thrust  propulsion  the  spacecraft 
naturally  spirals  around  the  body  it  is  orbiting,  and  in  this  case  that  body  is  the  sun. 


Figure  4.15:  Pre -Jupiter  Trajectory  (Distance  -  AU)(Angle  -  degrees) 

Travelling  from  the  Earth  to  Jupiter  on  this  trajectory  requires  16.9  years,  and  during 
those  years  the  spacecraft  is  thrusting  the  entire  time.  It  is  important  to  note  that  while  the 
spacecraft  is  thrusting  the  entire  time  it  is  also  being  throttled  to  lower  powers  to 
compensate  for  the  dwindling  solar  flux.  The  code  in  Appendix  A.l  shows  how  this  was 
simulated. 

This  trajectory  requires  17.17  km/ sec  of  AV,  which  corresponds  to  a  propulsive  mass 
of  2996.2  kg.  Hence,  the  final  spacecraft  mass  for  this  phase  is  5603.8  kg.  Finally,  the 
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spacecraft  velocity  at  Jupiter  is  13.68  km/ sec.  The  entire  velocity  profile  from  Earth  to 
Jupiter  is  shown  in  Figure  4.16.  In  the  figure  the  velocity  of  the  spacecraft  is  on  a  constant 
decline  due  to  the  pull  of  the  sun’s  gravitation.  However,  the  thrust  of  the  spacecraft 
decreases  the  deceleration  effect  of  the  sun  allowing  the  spacecraft  to  push  itself  all  the 
way  to  Jupiter.  Finally,  the  jagged  portion  of  the  profile  is  the  result  of  the  spacecraft’s 
varying  thrust  levels. 


Figure  4.16:  Pre-Jupiter  Velocity  Profile 


4.53.2  Jupiter  Phase.  At  Jupiter,  the  spacecraft  enters  into  its  gravity  assist 
maneuver  and  drops  its  1200  kg  boost  stage  at  Jupiter’s  periapsis.  As  a  result,  at  the  end  of 
the  gravity  assist  maneuver  the  spacecraft  mass  is  down  to  4403.8  kg.  All  of  the  gravity 
assist  data  was  derived  from  the  equations  presented  in  Section  3.5.2.  As  expected,  the 
velocities  of  the  spacecraft  relative  to  Jupiter  entering  and  exiting  the  maneuver  are  equal, 
but  the  velocities  relative  to  the  Sun  are  distinctly  different.  In  fact,  relative  to  the  Sun  the 
spacecraft  experiences  a  velocity  increase  of  13.905  km/ sec  bringing  its  exit  velocity  to 
27.585  km/ sec. 
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4.53.3  Post-Jupiter  Phase.  The  final  phase  characterizes  the  trajectory  from 
Jupiter  to  550  AU.  The  initial  velocity  at  Jupiter  is  27.585  km/ sec.  By  the  time  the 
spacecraft  reaches  550  AU  its  velocity  has  increased  to  29.6  km/ sec,  or  6.24  AU  /year. 
Figure  4.17  presents  the  complete  velocity  profile. 


Figure  4.17:  Post-Jupiter  Velocity  Profile 


Looking  at  the  figure,  the  velocity  of  the  spacecraft  continues  to  decline  until  a 
minimum  velocity  is  reached.  This  is  due  to  the  waning  effect  of  the  sun’s  gravitation  and 
the  constant,  long  duration  thrust  of  the  propulsion  system.  Prior  to  the  minimum  the 
sun’s  gravitational  force  is  dominant  and  decelerates  the  spacecraft.  After  the  minimum 
the  thrust  of  the  propulsion  system  dominates,  accelerating  the  spacecraft  back  to 
approximately  29.6  km/ sec.  Moreover,  the  velocity  profile  corresponds  to  a  AV  of  20.72 
km/ sec  and  a  thrust  period  of  16.79  years. 


4.53.4  Total.  These  three  phases  result  in  a  total  AV  of  5 1 .79  km/ sec,  total 
thrust  period  of  33.69  years,  and  total  trip  time  of  108.72  years.  This  may  appear  to  be  a 
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Figure  4.18:  Post-Jupiter  Trajectory  (Distance  -  AU)(Angle  -  degrees) 

daunting  trip  time  but  compare  that  to  Voyager  I.  It  was  launched  35  years  ago  in  1977 
and  it  would  still  take  Voyager  I  another  1 19  years  to  reach  550  AU  at  3.599  AU/year  [48]. 

4.5.4  Communications.  Figure  4.19  depicts  the  data  rate  for  the  mission  assuming 
the  communications  architecture  is  receiving  a  full  15  W  for  transmit.  That  transmitter 
power  corresponds  to  an  EIRP  of  15225.19  W  (41.82  dB)  for  the  MGA  and  548106.7  W 
(57.39  dB)  for  the  HGA.  Moreover,  the  data  rate  in  Figure  4.19  was  calculated  using 
Equation  3.10  with  a  spacecraft  BER  of  10  5,  BPSK,  and  Plus  RS  Viterbi  Decoding. 
Finally,  Figure  4.19  show  the  link  requirement.  If  the  HGA  and  MGA  are  above  the 
minimum  there  is  link  margin.  On  the  other  hand,  if  it  dips  below  the  minimum  it  is  no 
longer  meeting  the  data  transfer  requirements.  Therefore,  at  approximately  250  AU  the 
MGA  no  longer  meets  the  data  transfer  rate.  This  was  expected  and  by  this  time  the  HGA 
will  be  handling  data  transfer. 
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Figure  4.19:  Mission  Data  Rate 


4.5.5  Cost.  Cost  can  be  estimated  a  number  of  ways.  One  approach  is  to  utilize 
cost  models  and  that  is  what  was  done  in  this  thesis.  Two  models  were  examined.  The  first 
model  investigated  was  the  NASA  Advanced  Missions  Cost  Model  (AMCM).  Its 
self-described  purpose  is  to  ’’provide  a  useful  method  for  quick  turnaround, 
rough-order-of-magnitude  [cost  estimation]”  [49].  The  inputs  to  the  model  are  quantity, 
dry  weight,  mission  type,  launch  year,  block  number,  and  difficulty.  For  the  FSM,  the 
quantity  is  one,  the  dry  weight  is  2345  kg,  the  mission  type  is  ’’spacecraft  -  planetary”,  the 
launch  year  was  chosen  to  be  2020,  the  block  number  (which  is  described  as  the  level  of 
inheritance  from  previous  spacecraft)  is  one  since  the  FSM  constitutes  a  new  system,  and 
the  difficulty  is  very  high  [49].  The  cost  estimation  equation  used  by  the  AMCM  model  is 

Cost  =  a  *  Qb  *  Wc  *  ds  *  e'oc-isoo  *  g  f  *  gD  [50]  (4.4) 
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where  Q  is  the  number  of  spacecraft,  W  is  the  weight,  s  is  a  value  derived  from  the 
mission  type,  IOC  is  the  launch  year,  B  is  is  block  number,  D  is  a  value  derived  from  the 
difficulty,  and  a-g  are  model  parameters  [50].  The  resulting  cost  is  $3,927  billion  (in 
FY11  US$). 

The  total  cost  is  expected  to  be  in  the  range  of  $3-5  billion.  An  exact  cost  cannot  be 
tied  down  due  to  cost  overruns,  inflation,  and  the  nebulous  costs  of  labor.  Nonetheless,  the 
pay  of  a  single  GS-14  employee  was  analyzed  to  get  a  general  idea  of  these  costs.  In  2012, 
a  GS-14  working  in  Denver,  CO  will  make  $103,771  [51].  In  the  108  years  that  it  will  take 
the  spacecraft  to  reach  550  AU  and  given  a  conservative  2%  raise  per  year  that  will 
increase  to  $880,835.  That  means,  assuming  spacecraft  were  launched  today  with  a 
complement  of  20  GS-14  engineers,  labor  costs  in  2120,  when  the  spacecraft  arrives, 
would  be  approximately  $17  million  per  year.  What  is  more,  integrating  the  cost  from  the 
launch  date  to  108  years  shows  that  the  total  cost  over  that  time  span  of  all  twenty 
positions  is  approximately  $800  million.  Outside  of  labor  costs  there  are  also  operating 
costs  like  electricity,  building  maintenance,  and  equipment  replacement  costs  which  are 
outside  of  the  scope  of  this  analysis. 
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5  Conclusions  and  Recommendations 


5.1  Introduction 

Outside  of  the  deliverables  this  thesis  had  three  goals.  The  first  goal  was  to  produce  a 
piece  of  work  with  academic  merit.  The  second  goal  was  to  present  a  novel  design.  The 
final  goal  was  to  present  actionable  recommendations  to  the  scientific  community.  The 
first  four  chapters  of  this  thesis  accomplished  the  first  two  goals  and  this  final  chapter 
accomplishes  the  third  goal  by  presenting  the  conclusions  of  this  research  and  analyzing 
them  to  come  up  with  actionable  recommendations. 

5.2  Conclusions 

In  simple  terms,  this  thesis  demonstrated  that  it  is  possible  to  create  an  interstellar 
precursor  spacecraft  solely  from  commercial  of  the  shelf  parts.  Unfortunately,  it  would 
take  more  than  a  century  for  the  spacecraft  to  reach  the  desired  distances  meaning  there 
would  be  significant  risk  and  astronomical  costs  involved.  Also,  that  century  would  eclipse 
the  span  of  several  mission  directors  and  engineers  lifetimes,  which  would  introduce  even 
more  risk.  Furthermore,  if  this  mission  were  to  fly  successfully  in  spite  of  the  obstacles 
the  total  gain  from  the  sun  and  the  12-m  antenna  would  only  be  1 15  dB.  A  terrestrial 
based  8.4  km  radio  telescope  would  have  the  same  gain,  and  there  are  already  telescopes 
of  that  size.  The  gravitational  lens  still  boasts  superior  angular  resolution,  staring  time, 
and  a  clear  view  of  its  target  without  atmospheric  interference.  Nonetheless,  there  is  no 
need  to  explore  this  mission  further  at  this  point  because  of  the  costs  and  trip  times 
involved.  Still,  this  thesis  provided  insight  into  what  the  current  deep  space  capabilities 
are,  and  anyone  interested  in  deep  space  flight  could  glean  helpful  information  from  it. 
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5.3  Recommendations 


•  The  first  recommendation  is  to  spend  more  effort  deveioping  spacecraft  power 
instead  of  new  thrusters.  At  this  moment,  power  is  the  iimiting  factor  in  deep  space 
expioration;  thus,  it  wouid  most  significantiy  improve  deep  space  flight.  For 
exampie,  if  the  power  density  was  increased  so  that  the  NEXT  thruster  couid 
maintain  maximum  thrust  for  the  entire  FSM  mission  the  trip  time  wouid  decrease 
by  neariy  10%  to  98  years.  If  two  thrusters  could  be  powered  there  would  be  a  25% 
decrease  in  trip  time  to  80  years.  This  is  because  at  higher  powers  the  thruster  can 
maintain  higher  thrust  and  use  the  propellant  mass  more  efficiently  with  higher 
specific  impulses.  Spacecraft  power  can  be  developed  through  the  exploration  of 
new  radioisotopes  and  higher  efficiency  solar  panels  or  the  introduction  of  new 
power  methods,  i.e.  nuclear  reactors. 

•  Again  this  recommendation  relates  to  power.  Though  the  exact  amount  of 
plutonium-238  the  government  possesses  is  a  secret,  it  is  estimated  that  there  may 
only  be  enough  plutonium-238  for  one  more  mission.  [52]  Hence,  the  growing 
difficulty  in  procuring  radioisotopes  necessitates  developing  a  novel  power  method 
such  as  nuclear  reactors  or  restarting  radioisotope  production  if  deep  space 
exploration  is  to  continue. 

•  In  this  thesis,  the  long  term  reliability  of  all  instruments  was  taken  on  the  fact  that 
the  instruments  chosen  are  manufactured  to  withstand  the  radiation  that  the  FSM 
will  experience.  Moreover,  it  would  not  be  feasible  to  do  a  long  duration  test  that 
lasts  more  than  a  decade  because  by  that  time  the  device  will  be  out  of  date.  Hence, 
a  novel  long  duration  test  method  must  be  developed. 

•  In  regard  to  communications,  new  modulation  schemes  must  continue  to  be 
explored  to  push  data  rates  closer  to  their  Shannon  Limit. 
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•  Increase  emphasis  on  low  thrust  propulsion,  investigate  replacement  propellants 
because  of  the  difficulty  in  procuring  xenon.  Xenon  is  a  rare  gas  and  it  is  only 
produced  at  air  separation  plants.  There  are  only  about  75  of  these  plants  in  the 
entire  world  [53].  These  75  plants  produce  just  over  9  million  liters  or  53,046  kg  of 
xenon  a  year  which  means  the  FSM,  and  any  similar  mission,  would  require  nearly 
10%  of  the  world’s  yearly  supply  of  xenon  [53].  Thus  xenon  is  not  a  sustainable 
propellant  option. 

•  The  continued  development  of  a  heavy  launch  vehicle  on  the  scale  of  NASA’s  Space 
Launch  System  would  enable  deep  space  probes  to  carry  significantly  more 
propellant  and  power  which  would  in  turn  decrease  the  trip  time  as  shown  in  the 
first  bullet. 
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Appendix:  MATLAB 


This  appendix  contains  all  of  the  MATLAB  code  used  during  my  analysis,  which  was 
mentioned  on  pg.  51  in  Chapter  III. 

A.l  Pre- Jupiter 

This  portion  of  the  code  takes  in  the  initial  mass,  velocity,  and  orbital  radius  relative 
to  the  sun  as  inputs.  It  then  uses  loops  to  calculate  the  AV,  trip  time,  trajectory,  and  final 
mass  of  the  spacecraft  at  the  Jovian  orbit. 


1  clc;clear  all 

2 

3  %%  Initialize  Variables 

4  %  In  this  section  the  initial  coniditions  and  global  variables  . . . 

are  set 

5 

6  global  I_sp  m_f  thrust  g_0 


7 


8 

g_0  =  9.80665; 

% [m/s"2] 

9 

excess  =  0; 

% [km/ sec 

10 

11 

m_f  =  2514.23; 

%  [kg] 

12 

m_i  =  8600; 

%  [kg] 

13 

14  mass.count  =  0; 

15 

16  %%  Evaluation  Loop 

17  %The  primary  calculations  are  made  within  these  loops. 

18  for  m_i  =  m_i 

19  x_0  =  [29 . 78+excess  0  149598261  0  m_i];%ds  %dr  %r  %th  %m 
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20 


21  % 


Time  Step  Parameters 


22  t  _  0  =  0 ; 


o_ 

o 


[sec] 


23  t_f  =  31536000/10; 


o_ 

o 


[  sec] 


24 

25  %Counter/condition  for  "while"  loop 

26  count  =  0 ; 

27  r  =  x_0  (3)  ; 

28 

29  %Counter  for  "for"  loop 

30  massxount  =  massxount  +  1; 

31 

32  %Inital  Thrust  Level 

33  TL  =  1;  %  [counter] 

34  I_sp  =  4190;  %  [sec] 

35  thrust  =  236e-3;  %[N] 

36  input_power  =  7.22;  %[kW] 

37 

38  %Initial  Conditions  for  trajectory  plotting 

39  r_prev{mass_count}  =  149598261;  %[km] 

40  th_prev{mass_count}  =  0;  %[rad] 

41  v_prev{mass_count}  =  sqrt  (x_0  (1)  ~2+x_0  (2)  ~2)  ;  %  [km/sec] 

42 

43  %While  loop  that  executes  up  to  the  orbit  of  Jupiter 

44  while  r  <  5.204*149598261  % 

45 

46  %ODE45  Executor.  Integrates  Equations  of  Motion. 

47  [t , x] =ode4  5 ( @ derivative, linspace (t_0 , t_f, 100)  ,  x_0)  ; 

48 

49  %Saves  the  last  entry  and  stores  it  as  new  intial  condition 

50  [last_entry  -i]  =  size  (x)  ; 

51  x_0  =  x  (last-entry,  : )  ; 


100 


%radius  [km} 


52  r  =  x_0  (3)  ; 


53 

m_c  =  x_0  (5)  ;  %current  mass  [kg] 

54 

55 

%Throttling 

56 

power  =  (  (1/ (r/149598261) ) *1.37) *60*. 28; 

%  [kW] 

57 

if  power  <  input_power 

58 

switch  TL 

59 

casejl} 

60 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

61 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

62 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

63 

64 

%New  TL  Properties 

65 

I  ~sp  =  4150;  % [ sec] 

66 

thrust  =  208e-3;  % [N] 

67 

input.power  =  6.39;  %[kW] 

68 

69 

%  Sets  current  mass  and  the  new 

inital  mass 

70 

m_i  =  m_c;  %#ok<FXSET> 

71 

display (' level  change') 

72 

73 

74 

case{2 } 

75 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

76 

tau (TL)  =  . . . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

77 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

78 

79 

%New  TL  Properties 

80 

I_sp  =  3910;  %[sec] 

81 

thrust  =  221e-3;  % [N] 
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82 

input_power  =  6.385;  %[kW] 

83 

84 

%  Sets  current  mass  and  the  new 

inital  mass 

85 

m_i  =  m_c;  %#ok<FXSET> 

86 

display (' level  change') 

87 

88 

case{3} 

89 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

90 

tau (TL)  =  . . . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

91 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

92 

93 

%New  TL  Properties 

94 

I -Sp  =  4150;  %  [  sec ] 

95 

thrust  =  208e-3;  % [N] 

96 

input-power  =  6.39;  %[kW] 

97 

98 

%  Sets  current  mass  and  the  new 

inital  mass 

99 

m_i  =  m_c ;  %#ok<FXSET> 

100 

display (' level  change') 

101 

102 

case{ 4 } 

103 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

104 

tau (TL)  =  . . . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

105 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

106 

107 

%New  TL  Properties 

108 

I_sp  =  3690;  %  [sec] 

109 

thrust  =  208e-3;  % [N] 

110 

input-power  =  5.78;  %[kW] 

102 


%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 


116 

case{5} 

117 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

118 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

119 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

120 

121 

%New  TL  Properties 

122 

I_sp  =  3875;  %  [sec] 

123 

thrust  =  194e-3;  % [N] 

124 

input_power  =  5.66;  %[kW] 

125 

126 

%  Sets  current  mass  and  the  new 

inital  mass 

127 

m_i  =  m_c;  %#ok<FXSET> 

128 

display (' level  change') 

129 

130 

case{ 6} 

131 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

132 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

133 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

134 

135 

%New  TL  Properties 

136 

I  ~sp  =  4150;  %  [  sec ] 

137 

thrust  =  181e— 3;  % [N] 

138 

input-power  =  5.6;  %[kW] 

139 

140 

%  Sets  current  mass  and  the  new 

inital  mass 

m_i  =  m_c; 


5#ok<FXSET> 
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142 

display (' level  change') 

143 

144 

case{ 7 } 

145 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

146 

tau (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

147 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

148 

149 

%New  TL  Properties 

150 

I_sp  =  3660;  % [sec] 

151 

thrust  =  184e-3;  % [N] 

152 

input_power  =  5.085;  %[kW] 

153 

154 

%  Sets  current  mass  and  the  new 

inital  mass 

155 

m_i  =  m_c;  %#ok<FXSET> 

156 

display (' level  change') 

157 

158 

case{8} 

159 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

160 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

161 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

162 

163 

%New  TL  Properties 

164 

I_sp  =  3395;  %  [sec] 

165 

thrust  =  192e-3;  % [N] 

166 

input-power  =  4.965;  %[kW] 

167 

168 

%  Sets  current  mass  and  the  new 

inital  mass 

169 

m_i  =  m_c;  %#ok<FXSET> 

170 

display (' level  change') 

171 
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172 

case{ 9} 

173 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

174 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

175 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c) ) /1000; % [km/ sec] 

176 

177 

%New  TL  Properties 

178 

I_sp  =  3875;  %  [sec] 

179 

thrust  =  169e-3;  % [N] 

180 

input_power  =  4.92;  %[kW] 

181 

182 

%  Sets  current  mass  and  the  new 

inital  mass 

183 

m_i  =  m_c;  %#ok<FXSET> 

184 

display (' level  change') 

185 

186 

case{ 1 0} 

187 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

188 

tau  (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

189 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

190 

191 

%New  TL  Properties 

192 

I -Sp  =  4100;  % [ sec ] 

193 

thrust  =  158e-3;  % [N] 

194 

input_power  =  4.87;  %[kW] 

195 

196 

%  Sets  current  mass  and  the  new 

inital  mass 

197 

m_i  =  m_c ;  %#ok<FXSET> 

198 

display (' level  change') 

200 

201 


case{ 1 1} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 


105 


202  tau  (TL)  =  ... 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) 731536000;% [years] 


203 

Delta_V (TL)  =  (I_sp 

*g_0*log (m_i/m_c) ) /1000; % [km/ sec] 

204 

205 

%New  TL  Properties 

206 

I_sp  =  3360; 

% [ sec] 

207 

thrust  =  169e-3; 

%  [N] 

208 

input_power  =  4.49; 

%  [kW] 

209 

210 

%  Sets  current  mass 

and  the  new  inital  mass 

211  m_i  =  m_c;  %#ok<FXSET> 

212  display  (' level  change') 

213 

214  case{l2} 

215  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

216  tau  (TL)  =  ... 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 

217  Delta_V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c ))/ 1000 ;% [km/sec] 

218 

219  %New  TL  Properties 

220  I_sp  =  3660;  %  [sec] 

221  thrust  =  160e-3;  %  [N] 

222  input  _power  =  4.4  55;  %[kW] 

223 

224  %  Sets  current  mass  and  the  new  inital  mass 

225  m_i  =  m_c;  %#ok<FXSET> 

226  display  (' level  change') 

227 

228  case{l3} 

229  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

230  tau  (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) 731536000;% [years] 
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(I_sp*g_0*log (m_i /m_c ) ) /1000 ;  % [km/ sec] 


23i  Delta_V  (TL)  = 


232 

233  %New  TL  Properties 


234 

I_sp  =  3830;  %  [sec] 

235 

thrust  =  147e— 3;  % [N] 

236 

input_power  =  4.315;  %[kW] 

237 

238 

%  Sets  current  mass  and  the  new 

inital  mass 

239 

m_i  =  m_c;  %#ok<FXSET> 

240 

display (' level  change') 

241 

242 

case{ 1 4} 

243 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

244 

tau  (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

245 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

246 

247 

%New  TL  Properties 

248 

I_sp  =  4310;  %[sec] 

249 

thrust  =  134e-3;  % [N] 

250 

input_power  =  4.235;  %[kW] 

251 

252 

%  Sets  current  mass  and  the  new 

inital  mass 

253 

m_i  =  m_c ;  %#ok<FXSET> 

254 

display (' level  change') 

255 

256 

case{ 1 5} 

257 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

258 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

259 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

260 
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261 

%New  TL  Properties 

262 

I  ~sp  =  3615;  % [ sec ] 

263 

thrust  =  139e-3;  % [N] 

264 

input-power  =  3.91;  %[kW] 

265 

266 

%  Sets  current  mass  and  the  new 

inital  mass 

267 

m_i  =  m_c;  %#ok<FXSET> 

268 

display (' level  change') 

269 

270 

case{ 1 6} 

271 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

272 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

273 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

274 

275 

%New  TL  Properties 

276 

I_sp  =  3360;  %  [sec] 

277 

thrust  =  147e— 3;  % [N] 

278 

input.power  =  3.86;  %[kW] 

279 

280 

%  Sets  current  mass  and  the  new 

inital  mass 

281 

m_i  =  m_c ;  %#ok<FXSET> 

282 

display (' level  change') 

283 

284 

case{ 1 7 } 

285 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

286 

tau (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

287 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

288 

289 

%New  TL  Properties 

290 

I-sp  =  4025;  % [ sec ] 
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291 

thrust  =  125e-3;  % [N] 

292 

input-power  =  3.76;  %[kW] 

293 

294 

%  Sets  current  mass  and  the  new 

inital  mass 

295 

m_i  =  m_c ;  %#ok<FXSET> 

296 

display (' level  change') 

297 

298 

case{ 1 8} 

299 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

300 

tau (TL)  =  . . . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

301 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

302 

303 

%New  TL  Properties 

304 

I ~sp  =  4190;  %  [ sec ] 

305 

thrust  =  107e-3;  % [N] 

306 

input_power  =  3.46;  %[kW] 

307 

308 

%  Sets  current  mass  and  the  new 

inital  mass 

309 

m_i  =  m_c;  %#ok<FXSET> 

310 

display (' level  change') 

311 

312 

case{ 1 9} 

313 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

314 

tau  (TL)  =  .  .  . 

(  (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

315 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

316 

317 

%New  TL  Properties 

318 

I ~sp  =  3125;  %  [ sec ] 

319 

thrust  =  137e-3;  % [N] 

320 

input_power  =  3.425;  %[kW] 
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321 


322 

%  Sets  current  mass  and  the  new  inital  mass 

323 

m_i  =  m_c;  %#ok<FXSET> 

324 

display (' level  change') 

325 

326 

case{2  0} 

327 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

328 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

329 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % 

[km/ sec] 

330 

331 

%New  TL  Properties 

332 

I_sp  =  3800;  % [sec] 

333 

thrust  =  118e— 3;  % [N] 

334 

input_power  =  3.415;  %[kW] 

335 

336 

%  Sets  current  mass  and  the  new  inital  mass 

337 

m_i  =  m_c;  %#ok<FXSET> 

338 

display (' level  change') 

339 

340 

case{2 1} 

341 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

342 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

343 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

344 

345 

%New  TL  Properties 

346 

I_sp  =  3325;  %  [sec] 

347 

thrust  =  128e-3;  % [N] 

348 

input.power  =  3.39;  %[kW] 

349 

350 

%  Sets  current  mass  and  the  new  inital  mass 
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354 

355 

356 


m_i  =  m_c ;  %#ok<FXSET> 
display (' level  change') 

case{2  2} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 


357 

DeltaJV (TL)  =  (I_sp* 

g_0*log (m_i/m_c) ) /1000;  % 

[km/ sec] 

358 

359 

%New  TL  Properties 

360 

I_sp  =  3910; 

% [ sec] 

361 

thrust  =  99.9e-3; 

%  [N] 

362 

input-power  =  3.08; 

%  [kW] 

363 

364 

%  Sets  current  mass 

and  the  new  inital  mass 

365 

m_i  =  m_c;  %#ok<FXSET> 

366 

display (' level  change') 

367 

368 

case{2  3} 

369 

%Finds  the  Burn  Time 

and  Delta  V  from  previous  TL 

370 

tau  (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

371 

Delta_V (TL)  =  (I_sp* 

g_0*log (m_i/m_c) ) /1000;  % 

[km/ sec] 

372 

373 

%New  TL  Properties 

374 

I_sp  =  3090; 

% [ sec] 

375 

thrust  =  119e— 3; 

%[N] 

376 

input_power  =  3.01; 

%  [kW] 

377 

378 

%  Sets  current  mass 

and  the  new  inital  mass 

379 

m_i  =  m_c;  %#ok<FXSET> 

380 

display (' level  change') 
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381 


382 

case{2  4} 

383 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

384 

tau (TL)  =  . . . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

385 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

386 

387 

%New  TL  Properties 

388 

I_sp  =  3490;  %  [sec] 

389 

thrust  =  108e-3;  % [N] 

390 

input-power  =  2.97;  %[kW] 

391 

392 

%  Sets  current  mass  and  the  new 

inital  mass 

393 

m_i  =  m_c ;  %#ok<FXSET> 

394 

display (' level  change') 

395 

396 

case{2  5} 

397 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

398 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years] 

399 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

400 

401 

%New  TL  Properties 

402 

I_sp  =  3690;  %  [sec] 

403 

thrust  =  94.3e-3;  % [N] 

404 

input_power  =  2.765;  %[kW] 

405 

406 

%  Sets  current  mass  and  the  new 

inital  mass 

407 

m_i  =  m_c;  %#ok<FXSET> 

408 

display (' level  change') 

409 

410 

case{2  6} 
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411 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

412 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

413 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

414 

415 

%New  TL  Properties 

416 

I_sp  =  3250;  %  [sec] 

417 

thrust  =  101e-3;  %[N] 

418 

input-power  =  2.62;  %[kW] 

419 

420 

%  Sets  current  mass  and  the  new  inital  mass 

421 

m_i  =  m_c;  %#ok<FXSET> 

422 

display (' level  change') 

423 

424 

case{2  7} 

425 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

426 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

427 

Delta_V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

428 

429 

%New  TL  Properties 

430 

I_sp  =  4000;  %  [sec] 

431 

thrust  =  80.2e-3;  % [N] 

432 

input_power  =  2.585;  %[kW] 

433 

434 

%  Sets  current  mass  and  the  new  inital  mass 

435 

m_i  =  m_c ;  %#ok<FXSET> 

436 

display (' level  change') 

437 

438 

case{2  8} 

439 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
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442 

443 

444 


447 

448 


464 

465 


tau (TL)  =  . . . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c )) /1000 ;% [km/sec] 

%New  TL  Properties 
I_sp  =  3395;  %  [sec] 

thrust  =  86.7e-3;  % [N] 

input_power  =  2.415;  %[kW] 

%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 


452 

case{2  9} 

453 

%Finds  the  Burn  Time 

and  Delta  V 

from  previous  TL 

454 

tau (TL)  =  . . . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

455 

Delta_V(TL)  =  (I_sp* 

g_0*log (m_i/ 

m_c  )  )  / 1 0 0 0 ; % [km/ sec ] 

456 

457 

%New  TL  Properties 

458 

I_sp  =  3735; 

% [ sec] 

459 

thrust  =  74.9e-3; 

%[N] 

460 

input_power  =  2.3; 

%  [kW] 

%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 

case{3  0} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
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(I_sp*g_0*log (m_i /m_c ) ) /1000 ;  % [km/ sec] 


469  Delta_V  (TL)  = 

470 

471  %New  TL  Properties 


472 

I ~sp  =  3155;  %  [ sec ] 

473 

thrust  =  80.6e-3;  % [N] 

474 

input_power  =  2.16;  %[kW] 

475 

476 

%  Sets  current  mass  and  the  new 

inital  mass 

477 

m_i  =  m_c;  %#ok<FXSET> 

478 

display (' level  change') 

479 

480 

case{3l} 

481 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

482 

tau  (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

483 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

484 

485 

%New  TL  Properties 

486 

I_sp  =  3525;  %  [sec] 

487 

thrust  =  70.7e-3;  % [N] 

488 

input_power  =  2.09;  %[kW] 

489 

490 

%  Sets  current  mass  and  the  new 

inital  mass 

491 

m_i  =  m_c ;  %#ok<FXSET> 

492 

display (' level  change') 

493 

494 

case{32} 

495 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

496 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

497 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

498 
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%New  TL  Properties 


499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 


I_sp  =  3240;  %[sec] 

thrust  =  65e-3;  % [N] 

input-power  =  1.825;  %[kW] 

%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 

case{3  3} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
Delta_V  (TL)  =  ( l_sp*g_0 *log (m_i /m_c ) ) /1000; % [km/sec] 


%New  TL  Properties 

I_sp  =  3015; 

% [sec] 

thrust  =  60.4e-3; 

%[N] 

input.power  =  1.635; 

%  [kW] 

%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 

case{3  4} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau  (TL)  =  . . . 

( (m_i-m_c ) /  (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
Delta_V (TL)  =  ( l_sp*g_0 *log (m_i /m_c )) /1000; % [km/sec] 

%New  TL  Properties 

I-sp  =  2885;  %  [sec] 
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529  thrust  =  57.8e-3;  %  [N] 

530  input-power  =  1.52;  %[kW] 

531 

532  %  Sets  current  mass  and  the  new  inital  mass 

533  m_i  =  m_c;  %#ok<FXSET> 

534  display  (' level  change') 

535 

536  case{35} 

537  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

538  tau  (TL)  =  ... 


( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) 731536000;% [years] 


539 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

540 

541 

%New  TL  Properties 

542 

I_sp  =  2745;  %[sec] 

543 

thrust  =  55.1e— 3;  % [N] 

544 

input_power  =  1.415;  %[kW] 

545 

546 

%  Sets  current  mass  and  the  new 

inital  mass 

547 

m_i  =  m_c;  %#ok<FXSET> 

548 

display (' level  change') 

549 

550 

case{3  6} 

551 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

552 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 

553  Delta_V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c )) /1000 ;% [km/sec] 

554 

555  %New  TL  Properties 

556  I_sp  =  2450;  %  [sec] 

557  thrust  =  49.2e-3;  %  [N] 

558  input_power  =  1.210;  %[kW] 
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559 


560 

%  Sets  current  mass  and  the  new  inital  mass 

561 

m_i  =  m_c;  %#ok<FXSET> 

562 

display (' level  change') 

563 

564 

case{3  7 } 

565 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

566 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years [ 

567 

Delta_V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % 

[km/ sec] 

568 

569 

%New  TL  Properties 

570 

I  ~sp  =  2400;  % [ sec ] 

571 

thrust  =  48.1e— 3;  % [N] 

572 

input_power  =  1.175;  %[kW] 

573 

574 

%  Sets  current  mass  and  the  new  inital  mass 

575 

m_i  =  m_c;  %#ok<FXSET> 

576 

display (' level  change') 

577 

578 

case{3  8} 

579 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

580 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

581 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

582 

583 

%New  TL  Properties 

584 

I_sp  =  1855;  %  [sec] 

585 

thrust  =  37.2e-3;  % [N] 

586 

input_power  =  .865;  %[kW] 

587 

588 

%  Sets  current  mass  and  the  new  inital  mass 
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589 

m_i  =  m_c ;  %#ok<FXSET> 

590 

display (' level  change') 

591 

592 

case{3  9} 

593 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

594 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

595 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

596 

597 

%New  TL  Properties 

598 

I ~sp  =  1585;  %  [ sec ] 

599 

thrust  =  31.8e— 3;  % [N] 

600 

input-power  =  .74;  %[kW] 

601 

602 

%  Sets  current  mass  and  the  new 

inital  mass 

603 

m_i  =  m_c;  %#ok<FXSET> 

604 

display (' level  change') 

605 

606 

case{ 4  0} 

607 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

608 

tau  (TL)  =  ... 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

609 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

610 

611 

%New  TL  Properties 

612 

I ~sp  =  1585;  %  [ sec ] 

613 

thrust  =  31.8e— 3;  % [N] 

614 

input_power  =  .74;  %[kW] 

615 

616 

%  Sets  current  mass  and  the  new 

inital  mass 

617 

m_i  =  m_c;  %#ok<FXSET> 

618 

display (' level  change') 
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619 


620 

621 

622 


case{4 1}  %TL  01 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau (TL)  =  . . . 


( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
Delta_V (TL)  =  ( l_sp*g_0*log (m_i/m_c) ) /1000; % [km/sec] 

%New  TL  Properties 


626 

I_sp  = 

1400; 

%  [ sec] 

627 

thrust 

=  25. 5e- 

-3; 

%  [N] 

628 

input_power  = 

.  610; 

%  [kW] 

629 

630 

%  Sets 

current 

mass 

and  the 

631 

632 


m_i  =  m_c;  %#ok<FXSET> 
display (' level  change1) 

case{42}  %off 
tau  (TL)  =  .  .  . 


636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 


( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years] 


Delta_V(TL)  = 

( I-sp 

*g_0*log (m_i/m_c) ) /1000; % [km/ sec] 

%New  TL  Properties 

I_sp  =  1480; 

o. 

o 

[sec] 

thrust  =  0; 

o, 

o 

[N] 

input_power  = 

0;  % 

[kW] 

end 

%Go  to  Next  thrust  level (TL)  at  next  power  change 
TL  =  TL+1 ; 

end 
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649  %Trajectory 

650  r_save{mass_count}  =  vertcat  (  r_prev{mass_count}  ,  x  (  :  ,  3)  )  ; 

651  r_prev{mass_count}  =  r_save{mass_count} ; 

652 

653  th_save{mass_count}  =  vertcat (th_prev{mass_count}  ,  x (:, 4 ))  ; 

654  th_prev{mass_count}  =  th_save{mass_count}; 

655 

656  v_save  {mas  s_count }  =  ... 

vertcat ( v_prev{mas s .count }  ,  sqrt (x ( : ,  1)  . ~2+x ( : , 2 )  . "2 ) )  ; 

657  v_prev{mass_count}  =  v_save{mass_count} ; 

658 

659  count  =  count  +  1; 

660  end 

661 

662  %%  Final  Calculations 

663 

664  %Calcualtes  Burn  Time  and  Delta  V  for  final  Thrust  Level  if 

thruster  does 

665  %not  throttle  to  the  off  position 

666  if  TL<42 

667  tau  (TL)  =  (  (m_i-m_c)  /  (thrust/ ( l_sp*g_0  ))) /31536000;  %  [years] 

668  Delta_V(TL)  =  ( I_sp*  g_0  *log  (m_i /m_c  ))/ 1 000  ;%  [km/sec] 

669  end 

670 

671  %Totals 

672  Total.Time  (mass.count )  =  ( count *t_f ) /31536000 ; 

673  Total-DV  (mass-count)  =  sum  (Delta.V)  ;  %  [km/sec] 

674  Total.tau  (mass.count )  =  sum(tau);  %  [years] 

675 

676  %Store 

677  Tau{mass_count}  =  tau; 

678  DV{mass_count}  =  Delta.V; 


% [years ] 
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679  %Clear 


680  clear  tau  Delta_V 

681 

682  end 


A.2  Post- Jupiter 

This  block  of  code  picks  after  the  Jovian  gravity  assist.  It  uses  the  final  values  from 
the  gravity  assist  as  initial  values  and  calculates  the  trajectory,  trip  time,  and  AV  from 
Jupiter  to  the  mission  distance.  Finally,  to  ensure  that  the  code  skips  the  power  levels  the 
spacecraft  is  not  capable  of  achieving  the  thrust  levels  up  to  level  16  were  blacked  out. 


1  clc;clear  all 

2 

3  %%  Initialize  Variables 

4  %  In  this  section  the  initial  coniditions  and  global  variables  . . . 

are  set 

5 

6  global  I_sp  m_f  thrust  g_0 

7 

8  g_0  =  9.80665;  %[m/s"2] 

9  excess  =  0;  % [km/sec] 


10 

11 

m_f  =  2514.23; 

%  [kg] 

12 

m_i  =  4403.8; 

%  [kg] 

13 

14 

RTG_power  =  4.914; 

%  [kw] 

15 

initial-time  =  17.12; 

% [years 

16 

17 

%Counter 
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is  mass_count 


0; 


19 

20  %%  Evaluation  Loop 

21  %The  primary  calculations  are  made  within  these  loops. 

22  for  m_i  =  m_i 

23  x_0  =  [27.585  0  5.2*149598261  47.6497  m_i];%ds  %dr  %r  %th 

24 

25  %Time  Step  Parameters 

26  t_0  =  0;  %  [sec] 

27  t_f  =  31536000/10;  %[sec] 

28 

29  %Counter/condition  for  "while"  loop 

30  count  =  0; 

31  r  =  x_0  (3)  ; 

32 

33  %Counter  for  "for"  loop 

34  mass.count  =  mass.count  +  1; 

35 

36  %Inital  Thrust  Level 

37  TL  =  15;  %  [counter] 

38  I -Sp  =  3615;  %  [ sec] 

39  thrust  =  139e-3;  %[N] 


40 

input_power  =  3.91; 

%  [kW] 

41 

42 

%Initial  Conditions 

for  trajectory  plotting 

43 

r_prev{mass_count}  = 

:  5.2*149598261;  % 

[km] 

44 

th_prev{mass  .count } 

=  47.6497; 

% [ rad] 

45 

v_prev{mass_count}  = 

:  sqrt (x_0 (1) "2+x_0 (2) 

"2) ;  % [km/sec 

46 

47 

%While  loop  that  executes  up  to  the  orbit 

of  Jupiter 

48 

while  r  <  550*149598261 

49 
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50  %ODE45  Executor.  Integrates  Equations  of  Motion. 

51  [  t ,  x]  =ode4  5  ( @  derivative,  linspace  (t_0,t_f,100),x_0)  ; 

52 

53  %Saves  the  last  entry  and  stores  it  as  new  intial  condition 

54  [  last_entry  ->]  =  size(x); 

55  x_0  =  x  ( last_entry ,  :  )  ; 

56  r  =  x_0  (3)  ;  %radius  [km} 

57  m_c  =  x_0  (5)  ;  %current  mass  [kg] 

58 

59  %Throttling 

60  power  =  ... 

RTG_power* . 5 "  (  (  (  ( count *t_f) /31536000) +initial_time) /87 . 7)  ; 

%  [kW] 

61  if  power  <  input_power  +  .1 

62  switch  TL 

63  %  casejl} 

64  %  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

65  %  tau  (TL)  =  ... 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years] 

66  %  Delta_V  (TL)  =  ( I_sp  *  g_0  *log  (m_i /m_c  ))/ 1000 ;%  [km/sec] 

67  % 

68  %  %New  TL  Properties 

69%  I_sp  =  4150;  %[sec] 

70  %  thrust  =  208e—  3;  %  [N] 

71  %  input_power  =  6.39;  %[kW] 

72  % 

73  %  %  Sets  current  mass  and  the  new  inital  mass 

74  %  m_i  =  m_c ;  %#ok<FXSET> 

75  %  display  (' level  change') 

76  % 

77  % 

78  %  case{2} 
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79  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

80  % 

tau (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) 731536000;% [years 


81  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec] 

82  % 


83  % 

%New  TL  Properties 

84  % 

I_sp  =  3910;  %[sec] 

85  % 

thrust  =  221e-3;  % [N] 

86  % 

input-power  =  6.385;  %[kW] 

87  % 


88  % 

%  Sets  current  mass  and  the  new  inital  mass 

89  % 

m_i  =  m_c ;  %#ok<FXSET> 

90  % 

display (' level  change') 

91  % 

92  %  case{3} 


93  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

94  % 

tau (TL)  =  ... 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


95  % 

Delta.V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % [km/sec] 

96  % 


97  % 

%New  TL  Properties 

98  % 

I  -S  p  =  4150;  % [sec] 

99  % 

thrust  =  208e-3;  % [ N ] 

100  % 

input_power  =  6.39;  %[kW] 

101  % 


102  % 

%  Sets  current  mass  and  the  new  inital  mass 

103  % 

m_i  =  m_c;  %#ok<FXSET> 

104  % 

display (' level  change') 

105  % 

106  %  case{4} 

107  %  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
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108  %  tau  (TL)  =  ... 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


109  % 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec] 

no  % 


111  % 

%New  TL  Properties 

112  % 

I_sp  =  3690;  %[sec] 

113  % 

thrust  =  208e-3;  % [N] 

114  % 

input_power  =  5.78;  %[kW] 

115  % 


116  % 

%  Sets  current  mass  and  the  new  inital  mass 

117  % 

m_i  =  m_c;  %#ok<FXSET> 

118  % 

display (' level  change') 

119  % 

120  %  case{5} 


121  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

122  % 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


123  % 

Delta_V (TL)  =  ( l_sp*g_0*log (m_i/m_c) ) /1000;  % [km/sec] 

124  % 


125  % 

%New  TL  Properties 

126  % 

I_sp  =  3875;  %[sec] 

127  % 

thrust  =  194e-3;  % [N] 

128  % 

input_power  =  5.66;  %[kW] 

129  % 


130  % 

%  Sets  current  mass  and  the  new  inital  mass 

131  % 

m_i  =  m_c ;  %#ok<FXSET> 

132  % 

display (' level  change') 

133  % 

134  %  case{6} 


135  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

136  % 

tau (TL)  =  . . . 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 
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137  % 

Delta.V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec] 

138  % 


139  % 

%New  TL  Properties 

140  % 

I  -S  p  =  4150;  % [sec] 

141  % 

thrust  =  181e— 3;  % [N] 

142  % 

input_power  =  5.6;  %[kW] 

143  % 


144  % 

%  Sets  current  mass  and  the  new  inital  mass 

145  % 

m_i  =  m_c;  %#ok<FXSET> 

146  % 

display (' level  change') 

147  % 

ms  %  case{7} 


149  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

150  % 

tau (TL)  =  . . . 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


151  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec] 

152  % 


153  % 

%New  TL  Properties 

154  % 

I_sp  =  3660;  %[sec] 

155  % 

thrust  =  184e-3;  % [N] 

156  % 

input_power  =  5.085;  %[kW] 

157  % 


158  % 

%  Sets  current  mass  and  the  new  inital  mass 

159  % 

m_i  =  m_c;  %#ok<FXSET> 

160  % 

display (' level  change') 

161  % 

i62  %  case{8} 


163  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

164  % 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years 


165  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec ] 

166  % 
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167  % 

%New  TL  Properties 

168  % 

I_sp  =  3395;  %[sec] 

169  % 

thrust  =  192e-3;  % [N] 

170  % 

input-power  =  4.965;  %[kW] 

171  % 


172  % 

%  Sets  current  mass  and  the  new  inital  mass 

173  % 

m_i  =  m_c;  %#ok<FXSET> 

174  % 

display (' level  change') 

175  % 

i76  %  case{9} 


177  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

178  % 

tau (TL)  =  ... 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


179  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ; % [km/sec] 

180  % 


181  % 

%New  TL  Properties 

182  % 

I_sp  =  3875;  %[sec] 

183  % 

thrust  =  169e-3;  % [N] 

184  % 

input_power  =  4.92;  %[kW] 

185  % 


186  % 

%  Sets  current  mass  and  the  new  inital  mass 

187  % 

m_i  =  m_c;  %#ok<FXSET> 

188  % 

display (' level  change') 

189  % 

i9o  %  casejlO} 


191  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

192  % 

tau (TL)  =  . . . 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


193  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec] 

194  % 


195  % 

%New  TL  Properties 

196  % 

I  -S  p  =  4100;  % [sec] 
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197  % 

thrust  =  158e-3;  %[N] 

198  % 

input-power  =  4.87;  %[kW] 

199  % 


200  % 

%  Sets  current  mass  and  the  new  inital  mass 

201  % 

m_i  =  m_c ;  %#ok<FXSET> 

202  % 

display (' level  change') 

203  % 

204%  casejll} 


205  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

206  % 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


207  % 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ; % [km/sec] 

208  % 


209  % 

%New  TL  Properties 

210  % 

I_sp  =  3360;  %[sec] 

211  % 

thrust  =  169e-3;  % [N] 

212  % 

input_power  =  4.49;  %[kW] 

213  % 


214  % 

%  Sets  current  mass  and  the  new  inital  mass 

215  % 

m_i  =  m_c;  %#ok<FXSET> 

216  % 

display (' level  change') 

217  % 

218  %  case{l2} 


219  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

220  % 

tau (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


221  % 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % [km/sec] 

222  % 


223  % 

%New  TL  Properties 

224  % 

I_sp  =  3660;  %[sec] 

225  % 

thrust  =  160e-3;  % [N] 

226  % 

input_power  =  4.455;  %[kW] 
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227  % 


228  % 

%  Sets  current  mass  and  the  new  inital  mass 

229  % 

m_i  =  m_c;  %#ok<FXSET> 

230  % 

display (' level  change') 

231  % 

232  %  case{l3} 


233  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

234  % 

tau (TL)  =  . . . 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years 


235  % 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % [km/sec ] 

236  % 


237  % 

%New  TL  Properties 

238  % 

I_sp  =  3830;  %[sec] 

239  % 

thrust  =  147e-3;  % [N] 

240  % 

input_power  =  4.315;  %[kW] 

241  % 


242  % 

%  Sets  current  mass  and  the  new  inital  mass 

243  % 

m_i  =  m_c;  %#ok<FXSET> 

244  % 

display (' level  change') 

245  % 

246  %  case{l4} 


247  % 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

248  % 

tau (TL)  =  . . . 

( (m_i-m_c) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years 


249  % 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ; % [km/sec] 

250  % 


251  % 

%New  TL  Properties 

252  % 

I_sp  =  4310;  %  [sec] 

253  % 

thrust  =  134e-3;  % [N] 

254  % 

input_power  =  4.235;  %[kW] 

255  % 


256  % 

%  Sets  current  mass  and  the  new  inital  mass 
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257  %  m_i  =  m_c ;  %#ok<FXSET> 

258  %  display  (' level  change') 

259  % 

260  %  case{l5} 

261  %  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

262  %  tau  (TL)  =  ... 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years] 

263  %  Delta.V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c ))/ 1000 ;% [km/sec ] 

264  % 

265  %  %New  TL  Properties 

266  %  I_sp  =  3615; 

267  %  thrust  =  139e-3; 

268  %  input-power  =  3.91; 

269  % 

270  %  %  Sets  current  mass  and  the  new  inital  mass 

271  %  m_i  =  m_c ;  %#ok<FXSET> 

272  %  display  (' level  change') 

273 


274 

case{ 1 6} 

275 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

276 

tau  (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

277 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) / 1000 ;  % 

[km/ sec] 

278 

279 

%New  TL  Properties 

280 

I_sp  =  3360;  %  [sec] 

281 

thrust  =  147e— 3;  % [N] 

282 

input.power  =  3.86;  %[kW] 

283 

284 

%  Sets  current  mass  and  the  new  inital  mass 

285 

m_i  =  m_c ;  %#ok<FXSET> 

286 

display (' level  change') 

%  [sec] 
%[N] 

%  [kW] 
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287 

288 

case{ 1 7 } 

289 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

290 

tau (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

291 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

292 

293 

%New  TL  Properties 

294 

I ~sp  =  4025;  % [ sec ] 

295 

thrust  =  125e-3;  %[N] 

296 

input-power  =  3.76;  %[kW] 

297 

298 

%  Sets  current  mass  and  the  new 

inital  mass 

299 

m_i  =  m_c ;  %#ok<FXSET> 

300 

display (' level  change') 

301 

302 

case{l 8} 

303 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

304 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

305 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

306 

307 

%New  TL  Properties 

308 

I ~sp  =  4190;  %  [ sec ] 

309 

thrust  =  107e-3;  %[N] 

310 

input_power  =  3.46;  %[kW] 

311 

312 

%  Sets  current  mass  and  the  new 

inital  mass 

313 

m_i  =  m_c;  %#ok<FXSET> 

314 

display (' level  change') 

315 

316 

case{ 1 9} 
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317 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

318 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

319 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

320 

321 

%New  TL  Properties 

322 

I  -Sp  =  3125;  % [ sec ] 

323 

thrust  =  137e-3;  %[N] 

324 

input-power  =  3.425;  %[kW] 

325 

326 

%  Sets  current  mass  and  the  new  inital  mass 

327 

m_i  =  m_c;  %#ok<FXSET> 

328 

display (' level  change') 

329 

330 

case{2  0} 

331 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

332 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

333 

Delta_V  (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

334 

335 

%New  TL  Properties 

336 

I_sp  =  3800;  %  [sec] 

337 

thrust  =  118e-3;  %[N] 

338 

input_power  =  3.415;  %[kW] 

339 

340 

%  Sets  current  mass  and  the  new  inital  mass 

341 

m_i  =  m_c ;  %#ok<FXSET> 

342 

display (' level  change') 

343 

344 

case{2 1} 

345 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
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346 


tau (TL) 


( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) 731536000;% [years] 


347 

Delta_V (TL)  =  (I_sp 

*g_0*log (m_i/m_c) ) /1000;  % 

348 

349 

%New  TL  Properties 

350 

I_sp  =  3325; 

% [ sec] 

351 

thrust  =  128e-3; 

%[N] 

352 

input.power  =  3.39; 

%  [kW] 

353 

354 

%  Sets  current  mass 

and  the  new  inital  mass 

357 

358 

359 


m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 

case{2  2} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau (TL)  =  . . . 


361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 


( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
Delta_V (TL)  =  ( l_sp*g_0 *log (m_i /m_c ) ) /1000; % [km/sec] 

%New  TL  Properties 


I_sp  =  3910; 

%  [  sec 

thrust  =  99.9e-3; 

%[N] 

input_power  =  3.08; 

%  [kW] 

%  Sets  current  mass  and  the  new  inital  mass 
m_i  =  m_c;  %#ok<FXSET> 
display (' level  change') 

case{2  3} 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
tau  (TL)  =  .  . . 


( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;%[ years] 
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375 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % 

[km/ sec] 

376 

377 

%New  TL  Properties 

378 

I_sp  =  3090;  %  [sec] 

379 

thrust  =  119e— 3;  % [N] 

380 

input_power  =  3.01;  %[kW] 

381 

382 

%  Sets  current  mass  and  the  new 

inital  mass 

383 

m_i  =  m_c ;  %#ok<FXSET> 

384 

display (' level  change') 

385 

386 

case{2  4} 

387 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

388 

tau  (TL)  =  .  .  . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000 ; 

% [years ] 

389 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c) ) /1000; % 

[km/ sec] 

390 

391 

%New  TL  Properties 

392 

I_sp  =  3490;  %  [sec] 

393 

thrust  =  108e-3;  % [N] 

394 

input_power  =  2.97;  %[kW] 

395 

396 

%  Sets  current  mass  and  the  new 

inital  mass 

397 

m_i  =  m_c ;  %#ok<FXSET> 

398 

display (' level  change') 

399 

400 

case{2  5} 

401 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

402 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000 ; 

% [years ] 

403 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % 

[km/ sec] 

404 
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405 

%New  TL  Properties 

406 

I_sp  =  3690;  % [sec] 

407 

thrust  =  94.3e-3;  % [N] 

408 

input-power  =  2.765;  %[kW] 

409 

410 

%  Sets  current  mass  and  the  new 

inital  mass 

411 

m_i  =  m_c;  %#ok<FXSET> 

412 

display (' level  change') 

413 

414 

case{2  6} 

415 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

416 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

417 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

418 

419 

%New  TL  Properties 

420 

I_sp  =  3250;  %  [sec] 

421 

thrust  =  101e-3;  % [N] 

422 

input_power  =  2.62;  %[kW] 

423 

424 

%  Sets  current  mass  and  the  new 

inital  mass 

425 

m_i  =  m_c;  %#ok<FXSET> 

426 

display (' level  change') 

427 

428 

case{2  7} 

429 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

430 

tau (TL)  =  . . . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

431 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

432 

433 

%New  TL  Properties 

434 

I_sp  =  4000;  %  [sec] 
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435 

thrust  =  80.2e-3;  % [N] 

436 

input-power  =  2.585;  %[kW] 

437 

438 

%  Sets  current  mass  and  the  new 

inital  mass 

439 

m_i  =  m_c ;  %#ok<FXSET> 

440 

display (' level  change') 

441 

442 

case{2  8} 

443 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

444 

tau (TL)  =  . . . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

445 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec) 

446 

447 

%New  TL  Properties 

448 

I_sp  =  3395;  %  [sec] 

449 

thrust  =  86.7e-3;  % [N] 

450 

input_power  =  2.415;  %[kW] 

451 

452 

%  Sets  current  mass  and  the  new 

inital  mass 

453 

m_i  =  m_c;  %#ok<FXSET> 

454 

display (' level  change') 

455 

456 

case{2  9} 

457 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

458 

tau  (TL)  =  .  .  . 

(  (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

459 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec) 

460 

461 

%New  TL  Properties 

462 

I_sp  =  3735;  %  [sec] 

463 

thrust  =  74.9e-3;  % [N] 

464 

input_power  =  2.3;  %[kW] 
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465 


%  Sets  current  mass  and  the  new  inital  mass 


467 

m_i  =  m_c;  %#ok<FXSET> 

468 

display (' level  change') 

469 

470 

case{3  0} 

471 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

472 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

473 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

474 

475 

%New  TL  Properties 

476 

I ~sp  =  3155;  %  [ sec ] 

477 

thrust  =  80.6e-3;  % [N] 

478 

input_power  =  2.16;  %[kW] 

479 

480 

%  Sets  current  mass  and  the  new 

inital  mass 

481 

m_i  =  m_c;  %#ok<FXSET> 

482 

display (' level  change') 

483 

484 

case{3l} 

485 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

486 

tau (TL)  =  . . . 

( (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years 

487 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

488 

489 

%New  TL  Properties 

490 

I_sp  =  3525;  %  [sec] 

491 

thrust  =  70.7e-3;  % [N] 

492 

input_power  =  2.09;  %[kW] 

493 

494 

%  Sets  current  mass  and  the  new 

inital  mass 
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495 

m_i  =  m_c ;  %#ok<FXSET> 

496 

display (' level  change') 

497 

498 

case{32} 

499 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

500 

tau  (TL)  =  .  .  . 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

501 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

502 

503 

%New  TL  Properties 

504 

I_sp  =  3240;  %  [sec] 

505 

thrust  =  65e-3;  % [N] 

506 

input-power  =  1.825;  %[kW] 

507 

508 

%  Sets  current  mass  and  the  new 

inital  mass 

509 

m_i  =  m_c;  %#ok<FXSET> 

510 

display (' level  change') 

511 

512 

case{3  3} 

513 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

514 

tau  (TL)  =  ... 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

515 

DeltaiV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  ) / 1000; % [km/ sec] 

516 

517 

%New  TL  Properties 

518 

I ~sp  =  3015;  %  [ sec ] 

519 

thrust  =  60.4e-3;  % [N] 

520 

input_power  =  1.635;  %[kW] 

521 

522 

%  Sets  current  mass  and  the  new 

inital  mass 

523 

m_i  =  m_c;  %#ok<FXSET> 

524 

display (' level  change') 
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525 


526 

case{34} 

527 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

528 

tau  (TL)  =  .  .  . 

(  (m_i-m_c )/ (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

529 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c )  )  /1000; % [km/ sec] 

530 

531 

%New  TL  Properties 

532 

I  ~sp  =  2885;  %[sec] 

533 

thrust  =  57.8e-3;  % [N] 

534 

input-power  =  1.52;  %[kW] 

535 

536 

%  Sets  current  mass  and  the  new 

inital  mass 

537 

m_i  =  m_c ;  %#ok<FXSET> 

538 

display (' level  change') 

539 

540 

case{3  5} 

541 

%Finds  the  Burn  Time  and  Delta  V 

from  previous  TL 

542 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) 

) ) /31536000; % [years) 

543 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i / 

m_c ) ) /1000; % [km/ sec] 

544 

545 

%New  TL  Properties 

546 

I_sp  =  2745;  %[sec] 

547 

thrust  =  55.1e— 3;  % [N] 

548 

input_power  =  1.415;  %[kW] 

549 

550 

%  Sets  current  mass  and  the  new 

inital  mass 

551 

m_i  =  m_c;  %#ok<FXSET> 

552 

display (' level  change') 

553 

554 

case{3  6} 
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555 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

556 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

557 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

558 

559 

%New  TL  Properties 

560 

I -Sp  =  2450;  %  [  sec ] 

561 

thrust  =  49.2e-3;  % [ N ] 

562 

input-power  =  1.210;  %[kW] 

563 

564 

%  Sets  current  mass  and  the  new  inital  mass 

565 

m_i  =  m_c;  %#ok<FXSET> 

566 

display (' level  change') 

567 

568 

case{3  7 } 

569 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

570 

tau  (TL)  =  .  .  . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

571 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

572 

573 

%New  TL  Properties 

574 

I  ~sp  =  2400;  % [ sec ] 

575 

thrust  =  48.1e— 3;  % [N] 

576 

input_power  =  1.175;  %[kW] 

577 

578 

%  Sets  current  mass  and  the  new  inital  mass 

579 

m_i  =  m_c;  %#ok<FXSET> 

580 

display (' level  change') 

581 

582 

case{3  8} 

583 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 
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584 


tau (TL) 


(  (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 

585  Delta_V  (TL)  =  ( I_sp*  g_0  *log  (m_i /m_c  )) /1000 ;%  [km/sec] 

586 

587  %New  TL  Properties 


588 

I_sp  = 

1855; 

%  [ sec] 

589 

thrust 

=  37. 2e- 

-3; 

%  [N] 

590 

input_power  = 

.865; 

%  [kW] 

591 

592 

%  Sets 

current 

mass 

and  the  new  inital  mass 

593  m_i  =  m_c;  %#ok<FXSET> 

594  display  (' level  change') 

595 

596  case{3  9} 

597  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

598  tau  (TL)  =  ... 

( (m_i -m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 

599  DeltaJV  (TL)  =  ( I_sp*  g_0  *log  (m_i /m_c  ))/ 1000 ;%  [km/sec] 

600 

6oi  %New  TL  Properties 


602 

I_sp  = 

1585; 

%  [  sec] 

603 

thrust 

=  31. 8e- 

-3; 

%[N] 

604 

input_power  = 

•  74; 

%  [kW] 

605 

606 

%  Sets 

current 

mass 

and  the  new  inital  mass 

607  m_i  =  m_c;  %#ok<FXSET> 

608  display  (' level  change') 

609 

610  case{40} 

611  %Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

612  tau  (TL)  =  ... 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000;% [years] 
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(I_sp*g_0*log (m_i /m_c ) ) /1000 ;  % [km/ sec] 


613  Delta_V  (TL)  = 

614 

615  %New  TL  Properties 


616 

I  ~sp  =  1585;  % [ sec ] 

617 

thrust  =  31.8e— 3;  % [N] 

618 

input_power  =  .74;  %[kW] 

619 

620 

%  Sets  current  mass  and  the  new  inital  mass 

621 

m_i  =  m_c;  %#ok<FXSET> 

622 

display (' level  change') 

623 

624 

case{4 1}  %TL  01 

625 

%Finds  the  Burn  Time  and  Delta  V  from  previous  TL 

626 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

627 

DeltaJV (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ; % 

[km/ sec] 

628 

629 

%New  TL  Properties 

630 

I ~sp  =  1400;  %  [  sec ] 

631 

thrust  =  25.5e-3;  % [N] 

632 

input_power  =  .610;  %[kW] 

633 

634 

%  Sets  current  mass  and  the  new  inital  mass 

635 

m_i  =  m_c;  %#ok<FXSET> 

636 

display (' level  change') 

637 

638 

case{42}  %off 

639 

tau (TL)  =  . . . 

( (m_i-m_c ) / (thrust/ (l_sp*g_0) ) ) /31536000; 

% [years ] 

640 

Delta_V (TL)  =  ( I_sp* g_0 *log (m_i /m_c ) ) /1000 ;  % 

[km/ sec] 

641 

642  %New  TL  Properties 
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643 

I_sp  =  1480; 

%  [sec 

644 

thrust  =  0; 

%[N] 

645 

input.power  =  0; 

%  [kW] 

646 

647 

648  end 

649  %Go  to  Next  thrust  level  (TL)  at  next  power  change 

650  TL  =  TL+1; 

651  end 

652 

653  %Trajectory 

654  r_save{mass_count}  =  vertcat  ( r_prev{mass_count} ,  x  (  :  ,  3)  )  ; 

655  r_prev{mass_count}  =  r_save{mass_count} ; 

656 

657  th_save{mass_count}  =  vertcat  (th_prev{mass_count}  ,  x  (:,  4  ))  ; 

658  th_prev{mass_count}  =  th_save{mass_count}  ; 

659 

660  v_save{mass_count}  =  ... 

vertcat ( v_prev{mas  s .count }  ,  sqrt (x(:,l)  . ~2+x ( : , 2 )  . "2 ) )  ; 

661  v_prev{mass_count}  =  v_save{mass_count} ; 

662 

663  count  =  count  +  1; 

664  end 

665 

666  %%  Final  Calculations 

667 

668  %Calcualtes  Burn  Time  and  Delta  V  for  final  Thrust  Level  if  . 

thruster  does 

669  %not  throttle  to  the  off  position  due  to  lack  of  power 

670  %  if  TL<3 

671  %  tau(TL)  =  (  (m_i -m_c  )/ (thrust /( l_sp*g_0  ))) /3153  6000  ;%  [years  ] 

672  %  Delta.V  (TL)  =  ( I_sp*  g_0  *log  (m_i /m_c  )) /1000  ;%  [km/sec] 
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673 


end 


674 

675  %Totals 

676  Total-Time  (mass_count )  =  (count*t_f ) /31536000  +  initial_time;  %  [years] 

677  Total_DV  (mass_count )  =  sum  (Delta_V)  ;  %  [km/sec] 

678  Total_tau  (mass_count )  =  sum(tau);  %  [years] 

679 

680  %Store 

681  Tau{mass_count}  =  tau; 

682  DV{mass_count}  =  DeltalV; 

683  %Clear 

684  clear  tau  Delta_v 

685 

686  end 


A.3  ODE  Derivative 

This  final  block  of  code  displays  the  derivative  function  necessary  for  the  operation 
of  the  MATLAB  ODE  45  code. 


1  function  xdot=derivative (t ,  x) 

2 

3  global  I_sp  m_f  thrust  g_0 

4  mu_sun  =  132712440018; 

5  phic  =  0 ;  % [ rad] 

6 

7  ds  =  x  ( 1 ) ; 

8  dr  =  x  (2 ) ; 

9  r  =  x  (3)  ; 

10  m  =  x  ( 5 )  ; 

11 
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12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 


m_dot =  -thrust/ (l_sp*g_0); 


if  m+m_dot*t  >  m_f 

m_dot_norm  =  m_dot/m; 
a  =  thrust/x (5) /1000; 

else 

a  =  0  ; 
m_dot  =  0; 
m_dot_norm  =  0; 

end 

xdot  =  [— dr*ds/r+a*cos (phic) /  (l+m_dot_norm*t) ;  .  .  . 

ds~2/r— mu.sun/ r~2+a*sin  (phic)  /  (l+m_dot_norm*t)  ;  .  .  . 
x  (2 )  ; 
ds/r; . . . 
m_dot ;  .  .  . 

]  ; 

end 
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